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suffers a deformation. 


Le présent article représente un effort que nous 
avions fait 1l y a quelques années pour comprendre 
et expliquer aux spécialistes des piles atomiques les 
phénomènes de transport des neutrons dans les 
solides. Cet article est resté longtemps sous forme 
d’un rapport interne, cependant nous croyons faire 
œuvre utile en le publiant aujourd’hui. 


I. Introduction. — La longueur de diffusion, Z, 
est une des grandeurs fondamentales permettant de 
classer les matériaux de pile atomique utilisés 
comme ralentisseurs. En effet, elle représente, en 
symétrie plane, la longueur moyenne parcourue par 
un neutron à partir d’un plan x — 0 avant de dis- 
paraître. Les ralentisseurs seront d’autant meilleurs 
que la longueur de diffusion sera plus grande. Dans 
la théorie classique du ralentissement ZL est reliée 
à deux autres grandeurs fondamentales : le libre 
parcours moyen de transport À et le libre parcours 
moyen de capture A 


z= (ss) | 1 


1/2 
: =) (1-1) 
o 3N°? 6tr Oa 


où ox eb © sont respectivement les sections efli- 
caces microscopiques moyennes de transport et de 
capture et AV le nombre d’atomes par unité de 
volume. 
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Résumé. — On développe un calcul de la longueur de diffusion des poudres cristallines tenant 
compte de la différence de libre parcours entre les neutrons de longueur d’onde supérieure à la 
première longueur d’onde de Bragg et les neutrons de longueur d’onde inférieure. Ce calcul permet 
indirectement de déterminer la valeur de la section efficace de capture. On montre également 
l’existence d’une déformation du spectre de Maxwell. H 

Le calcul est appliqué au cas du béryllium et de la glucine. 


Abstract. — A calculation of the neutron diffusion length in a crystal is developed. We take 
into account the difference of mean free path between neutrons of wave length greater than the 
first Bragg wave length and the neutrons ofsmaller wawelength. This calculation allows an indirect 
determination of the absorption cross-section. 


We show also that the Maxwell distribution 


The calculation is then applied to the case of beryllium and beryllium oxide. 


La longueur de diffusion est une grandeur direc- 
tement mesurable [1](*) ; la section efficace de cap- 
ture dans les milieux faiblement capturants est 
difficile à déterminer expérimentalement. Pour la 
calculer il faut connaître la section efficace de 
transport. Celle-ci peut être soit mesurée direc- 
tement [2], soit calculée théoriquement 


Str(0) = Ototai(v) — 5a(v) < COS 0 (y) (1-2) 


où < cos 0 > représente la valeur moyenne de 
l’angle de déviation pour un neutron de vitesse v, 

Dans le cas d’une poudre cristalline on peut cal- 
culer la section efficace élastique [3] ainsi que la 
valeur de < cos 0 >, (cf. IVe partie de cet 
article), on peut donc, en théorie, déterminer la 
section efficace de transport, et connaissant L dé- 
duire ca. | 

Dans les corps cristallins, les réflexions élastiques 
sont commandées par la loi de Bragg. Au-dessus 
d’une certaine longueur d’onde, ces réflexions 
n'existent plus, et seules se produisent des chocs 
inélastiques. La section efficace tombe à une valeur 
très faible, car elle ne contient plus que les termes 
inélastiques et de capture. Il était possible de se 


(*) Les chiffres entre crochets renvoient à la bibliogra- 
phie en fin d’article. 
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demander quel était le comportement des neutrons 
de longueur d’onde supérieure à la longueur d’onde 
de coupure, et éventuellement s’il était encore 
possible de définir une longueur de diffusion. Ce 
problème a été résolu, et nous avons pu déterminer 
la déformation du spectre de Maxwell qui en 
résulte. 

Dans une première partie, qui reprend l’essentiel 
d’un article de J. Horowitz et A. Messiah [5], et le 
complète, on traite le cas d’un milieu homogène. 
Le cas particulier des poudres cristallines est envi- 
sagé dans la troisième partie. Pour terminer le 
calcul est appliqué au cas du beryllium et de la 
glucine. 


II. La longueur de diffusion d’un milieu homo- 
gène diffusant faiblement capturant. — Dans ce 
qui suit nous allons supposer une distribution de 
neutrons dans un milieu homogène, diffusant et 
faiblement capturant en agitation thermique. Nous 
prendrons le cas de la symétrie plane et nous 
essaierons de trouver la forme asymptotique de 
la diffusion des neutrons. | 

Comme nous cherchons à déterminer essentiel- 
lement une longueur de diffusion nous supposons 
que la densité de neutrons a pour forme asympto- 
tique 


e—0# f(v) (2-1) 


où le spectre f(v) est indépendant de la position 
et w est une constante indépendante de la vitesse. 
La longueur de diffusion est l’inverse de «w. 

Une telle forme asymptotique n’est pas du tout 
obligatoire. Dans la limite du modérateur de masse 
infinie la vitesse des neutrons reste constante au 
cours de la collision et le problème se réduit à une 
superposition de problèmes à une vitesse ; à chaque 
vitesse du neutron correspond une longueur de 
diffusion L(v), donnée par w(v) = V3YX(v) Y,(v) et 
l'expression asymptotique prend la forme e** f(v) ; 
f(v) étant la densité au point x = 0. 

(Z représente la section efficace 
scopique : Z — No). 

Dans ce cas extrême il n’est pas possible de 
parler d’une longueur de diffusion, cela tient au 
fait qu'il n’y à pas d'échange d’énergie entre les 
neutrons et le milieu diffuseur de sorte que le 
spectre initial n’est déformé que par la capture. 

Dans la limite d’une capture nulle, au contraire, 
il existe, à une distance suffisante de la source, un 
spectre de vitesse Maxwellien H{(v) correspondant à 
l'équilibre thermique. La densité cherchée est en 
fait M(v) ; elle correspond à une longueur de dif- 
fusion infinie. La densité à la forme (2-1), avec 
== 0. 

Les deux cas considérés ci-dessus sont les limites 
de deux situations opposées. La première est celle 
où le pouvoir de ralentissement ËX,, est beaucoup 
plus petit que la section: efficace d'absorption Z, de 


macro- 
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sorte que le neutron n’a pratiquement pas le temps +. 
de modifier sa vitesse avant d’être capturé. La 
seconde est la limite du cas 


Za < 6 Yse (2-2) 


Ze est la section efficace de diffusion. 

Dans ce cas, le neutron échange suffisamment 
d’énergie avec le diffuseur avant sa capture et un 
certain équilibre thermique peut se réaliser. C’est » 
dans ce cas que l’on s’attend à une densité du 
type (1), avec une longueur de diffusion 1/« que 
nous allons essayer de déterminer. Dans le cas d’un 
cristal, la masse des atomes diffusant peut 
être considérée comme infinie mais il existe égale- 
ment des collisions inélastiques qui se produisent 
avec échange d’un ou plusieurs phonons, si bien 
que les vitesses des neutrons sont modifiables. 

Loin des sources nous pouvons écrire une équa- 
tion donnant le bilan des neutrons de vitesse v : 
c’est l’équation intégro-différentielle de Boltzmann- 
Lorentz 


v.grad D(r,v) + v Z (v) D(æ, v) 


— fr, &) D(æ, v’) de — 0. (2-3) 
En symétrie plane, nous aurons : 
vu D Da, ©) + v E (o) D(x, v) 
— fre, +) D(x, v’) dv’ = 0 


, +) : est la densité en phase, 
v : est la vitesse du neutron, 
: le cosinus de l’angle de v avec l’axe des x. 

Dans le deuxième terme qui représente la dispa- 
rition des neutrons par chocs, 2 — XX, + X,. et 
Ze — Ze + Zine Somme de la section efficace 
élastique et inélastique macroscopique, 2 repré- 
sente donc la section efficace macroscopique totale, 
et nous la supposerons corrigée de l’effet Doppler 
(négligeable dans le cas d’un cristal car la masse 
apparente est pratiquement infinie). 

Le troisième terme représente l’apparition des 
neutrons de vitesse v par choc. P(v', v) représente 
le noyau de la loi de choc. 

P(v', v) possède deux propriétés remarquables : 
la première exprime que ÂM(v) est solution de 
l’équation (2-4) si la capture est nulle : 

v Ze (0) M(v) — ie P(v’, v) M{v’) dv’ (2-5) 
la seconde, dont la première est une conséquence, 
est l’expression du principe de micro-réversibilité : 


P(v’,v) M(v') = P(v, v’) M(v). (2-6) 
De plus, nous pouvons écrire que : 
D Xe (V) = (re, v’) dv’. (2-7) 
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Les intégrales sont étendues à toutes les valeurs 
de ou de v’. 


A. RÉSOLUTION LE L'ÉQUATION DE BOLTZMANN- 
LORENTZ SANS TRANSPORT. — Nous faisons l’hypo- 
thèse simplificatrice que P(v', v) ne dépend que 
des grandeurs de ® et v’ et non de leurs directions. 
Physiquement cela suppose la diffusion isotrope 
dans le système du laboratoire, ce qui aboutit à 
négliger tout effet de transport. 

L’équation s’écrit : 


Du S D(a, &) + D(x) o[o D + [rt v') dv’| 
— fre", D) D(æ, v’) dv’ = 0. (2-8) 


On a posé 
P(v0’,0} = 47 P{v’, nv). 


Pour résoudre cette équation, nous allons sup- 


poser la densité développable en polynômes de 
_ Legendre. Soit 


| D(x, v) —— D,(x, v) + UD (x, v) 


+ ... + Pi Difs, 0) + ... (2-9) 


a) Résolution de l’équation (2-8) par une méthode 
à un seul harmonique. — Supposons maintenant la 
distribution D(x, v) suffisamment isotrope pour que 
nous puissions ne retenir que les deux premiers 
termes. 

On a 


D(x, v) PS Dilz, v) GE uD;(x, v). (2-10) 


L’équation de Boltzmann-Lorentz sans transport 
devient 


vus (Dolz, 0) + Di(x, 0] 
+ [D,{z, 0) + uDi(x, v)] [o re f P(v, v’) dv'| 


— fre, v) Di(x, v’) dv’ = 0. (2-11) 

Par une manipulation analogue à la théorie élé- 
mentaire classique (intégration en uw de (2-11), 
multiplication de (2-11) par y et intégration), on 
élimine D (x, v) et u. Il vient : 


UD? 
T3 jai Dolr, D) + 0 Za Dof, v) + Difx, v) 


fe, v’) dv” — fre, v) D,(x, v’) du’ = 0. (2-12) 


Si nous cherchons pour D,(x, v) une solution de 
la forme e°* D(v) nous sommes ramenés au pro- 
blème de valeur propre : 

«2? 


3S vD(v) + v Za D(v) + D(v) [ Pt, v’) dv” 


— [rt p) D(v’) do’ — 0. (2-13) 
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x) Calcul direct de ©. — Intégrons l’équation 
(2-13) en v, sur toute l’étendue des vitesses : 


me Su 6 + fo >, (v) D(v) dv 
+ j. D(v) dv [ P{v, v’) dv’ — Î Jk P(v’, v) D{v') dv’ = 0. 


(2-14) 


Les deux intégrales doubles étant égales il vient : 


Re Re 
182 du = se > 
2 (v) Z 
où ®(v) — vD(v) représente le flux des particules 


de vitesse v. Cette équation équivaut à la rela- 
tion (1-1). 

On voit donc que pour calculer la longueur de 
diffusion on devra faire la moyenne de la section 
efficace de capture pondéré? par le flux, ainsi que 
celle de l’inverse de la section efficace totale pon- 
dérée par le flux. La longueur de diffusion corres- 
pondant à la valeur de © la plus petite possible, on 
voit que l’expression (2-15) sera valable à une dis- 
tance infiniment grande pour que dans ®(v) on 
puisse négliger l'influence de racines supérieures 
à ©. (2-15) suppose que l’on connaisse ®(v) ce qui 
n’est généralement pas le cas ; nous allons donc 
donner une autre manière de calculer ©. 

B) Calcul de w par une méthode variationnelle. — 
Dans l’équation (2-13) l’opérateur intégral agissant 
sur D(v') n’est pas symétrique. Effectuons le chan- 
gement de fonction 


D(v) = 4 ER ®(v) (2-16) 
il vient : 
2 
ot = fKke,e)Dw)aw (247) 
et le noyau, 
K(v’,v) — 2? 5(v  — +) 
S) 72 #. 
VÉCETAP NEUE CNT) 


est symétrique comme conséquence du principe de 
micro-réversibilité (équation (2-6)). 

w?/3 est la valeur propre correspondant à la 
fonction propre ®(v) de l’opérateur intégral de 
noyau symétrique Æ, on a donc l’expression sta- 
tionnaire. 

D(v) K(v’,v) D(v') dv dv’ 
& nr J (2-19) 
7 D?{v) du 


En prenant pour o(v) une fonction assez voisine 
de la solution rigoureuse, on aura w?/3 sans grande 
erreur. Une bonne fonction d’essai est celle qui 
correspond à la solution rigoureuse à capture nulle 
c’est-à-dire à la distribution M(v) 
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Soit 
Œ(v) — (1/4x) VoM(v)[E (v). 


L'intégration du numérateur de (2-19) est immé- 
diate (cf. 2,5) et il vient : 


Ÿ v Za (v) M{v) do 


Da (2-20) 
2 
L vM{{v) 4 


2 (v) 


(2-20) est équivalente à (2-15) ou l’on a fait 

_vM(v) = ®(v) mais ici nous voyons que dans le 
cas d’une capture faible nous pouvons prendre le 
spectre de Maxwell comme représentant la distri- 
bution neutronique et ne commettre sur «w qu’une 
erreur du second ordre. 

b) Résolution de l'équation de Boltzmann dans le 
cas ou D ne peut être considéré comme presque 1s0- 
trope. — Revenons maintenant à l’équation de 
Boltzmann sous sa forme (2-8). 

_«) Résolution de l’équation (2-8) par une méthode 
à n harmoniques. — Développons D(x, v) selon 
l’équation (2-9) et bornons nous au /ieme terme. 
Il vient : 


u) Di, v) + ... + P; Dilx, v)] 
+ DE [D,(x, v) + ... + Pi Difx, v)] 
—f[ ?w v) Di(x, v’) dv’ = 0. 


En appliquant la relation de récurrence des poly- 
nômes de Legendre : 


» n +1 n er 
UP —_ On Li Pen + on +1 Pier il vient 
d n A n 
de 3 == L3 Dy+ilx, d) + ne Dyr—ilt, 2] Pau) 


> .S P{u) Dir, b) — fre, D) Dilx, v') du’ = 0. 
GR (2-22) 


En multipliant par P,, n — 0,1 ... et en inté- 


grant 1l vient / + 1 équations 
Dit: = 0 


ù l 
Ee Der Di-lx, vi] + v 2 Difx, v) = 0 


? 
3x Æ = 8 Dn+ilx, v) + er Dh—1(x, vi] 


+ 0 Z Dylæ, v) = 0 (2-23) 


ol 
= 5 Dilx, v) Ce v 2 Di(x, v) 


— fre, b) Dofa, v') dv’ = 0. 
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Si on cherche une solution de la forme : 
D, (x, v) = D,(p}es0 11 

d [ox sera remplacé par — « et on pourra éliminer 
tous les D, jusqu’à D, et D, exclus. Il vient 


ee 
D» 
Di 2 
: 55 * 35 ? 
D ps 
sine 12 


RE ne 
Eat + 1) (2 —1) 
Si nous ne nous étions pas bornés au /ieme terme 
nous aurions eu une fraction continue dont le terme 
2 2 
n 


Z?2(2n + 1) (2n — 1) 


de rang n est : . L’équation 


finale s’écrit 
U w? 1 
aol sera) 
DST 


+ vE Di(v  P(v’ v) D,{v’} du’ = 0. (2-28 


On peut montrer que : 
_® 1 un 2 © 
AM ES ss 


Nous allons retrouver ces résultats par un calcul 
plus direct de ©. 

B) Résolution de l’équation de Boltzmann par une 
méthode directe. — L’équation (2-8) n’est qu’un cas 
particulier de l’équation (2-3). On peut l’écrire sans 
tenir compte de la symétrie plane : 


v.grad D(r,®) + vE (v) D(r,v) 


= af 2e, v) D(r,v’) dv’ dQ’ (2-25) 
avec : 


[pee à = pr, v); v— Qui OPEN 


dQ est l'élément d’angle solide. On sait que l’équa- 
tion : 
@.grad D + ED = f{r,v) 


équivaut à l’équation : 
(e,e] 
Dir v)= [L f(r — s Q, v) e—Etns ds 


l'intégration se faisant sur la demi-droite issue de r 
en sens opposé à Q. 


N°4 
L’équation (2-25) est donc équivalente à 


D(r, v) = :1 F ds e—Etus 


fre s &, v’) P(v', v) aie (2-26) 
0 j s LT ° 
Soit par intégration SUP 
Dir, 0) — Sun pt e— Eos ds 
(É D(r— sQ,v’) P(v’, v) dv. (2-27) 


On cherche une solution de (2-27) sous la forme : 
Dr, v) = D(v) ec 
(c’est-à-dire que nous rétablissons ici la symétrie 
plane) : 


Désignons par uw la composante de Q sur x il 
vient : 


Dir — sQ,v) — D(v) exp [— ox + osul]. 


En reportant dans (2-27) 
obtient l'équation sur « : 


e—®% s’élimine et on 


1 co co 
== fan f ds e Eau [ D(v') P(v’, v) dv 
D 0 0 
(2-28) 
dQ — 2x du 
15 eus de = = | si Z2— œu > 0 
0 Z — ou ; 
Or il faut [u| < 1 ce qui entraîne Z — © > 0. 
Dans ce cas : 
eh Li es Le ’ 
) =; fs 4 D{(v’) P(v’,v) du’. (2-29) 
Soit 
AD 
Dioj— 22 Lee e [?p v) du’. (2-30) 


Si l’on pose x = (w/XZ) < 1 l’équation (2-30) 
prend la forme 
AŒEErT 


—— — ax OÙ a = 
Ce" 


Z (v) v D(v) 
fre, v) D(v’) du 


5 Lo (2-31) 


cette équation a une racine avec 0 << x < 1 si 
eds 

Dans le cas de la distribution de Maxwell 
a — 2/2 est bien supérieur à 1. Remarquons 
cependant que D(v) ne peut être la distribution 
maxWellienne, en effet « doit être indépendant 
de v, ce qui n’est pas réalisé. Dans le cas de la 
capture nulle on retrouve la valeur «© — 0 en 
prenant le spectre de Maxwell pour D(v). 

La condition È — « > 0 est très importante car 
elle conditionne la validité du développement en 
harmoniques. Si & Z X, la méthode n’est plus 
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justifiée et en particulier il n’est plus possible d’uti-- 
liser le premier harmonique comme il a été fait à 
l’alinéa a). 


B. RÉSOLUTION DE L'ÉQUATION DE BOLTZMANN- 
LORENTZ AVEC TRANSPORT. — Examinons mainte- 
nant les corrections de transport. Nous supposons 
évidemment que la condition de convergence est 
remplie. La loi de choc P(v', ») ne peut être consi- 
dérée comme isotrope. Donnons lui la forme 
P(v', v ; 15), dépendant à la fois des grandeurs des 
vitesses initiale et finale et de l’angle qu’elles font 
entre elles (*) uw, désigne le cosinus de cet angle. 


On peut développer P en polynômes de Legendre : 


Pb’, 0 3 eo) = Polv’, v) + bo PAlv’, 0) + 


Les propriétés de micro-réversibilité (2-5) et 
l’équation (2-7) restent toujours valables. Le trai- 
tement exposé ci-dessus fondé sur la quasi isotropie 
de la densité en phase, et non l’isotropie de la loi 
de choc peut être utilisé. Il conduira donc à une 
équation analogue à l’équation (2-13) pour D,;. 


1 2° 
JOLE 


(2-32) 


Da = Q:1Q0D (2-33) 


où @, et ©, sont des opérateurs intégraux de 


noyaux respectifs : 


et 


Supposons que P(v', v) présente un maximum 
aigu si v' — vet que l’on puisse remplacer : 


l'opération : 
2 
[£ P,(v’, v) dv 


It) glo) = eto) [É PAtv’, v) à 


Le sens physique de l’intégrale Z(v) est clair : 


= [EP v) dv 
— L ; 
v 
{ 
= fu Pt 


où << u, > est le cosinus moyen de l’angle de dé- 
flection dans le choc. L'opération Q, se réduit alors 
à une simple multiplication par la section efficace 
de transport : 


Zr (0) = Z (v) — Es (D) < uo > 


par : 


US Uo) dv’ — DE < Ho > 


(*) Ce n’est-pas la loi de choc la plus générale. En milieu 
cristallin la loi de choc peut aussi dépendre des directions 
des vitesses relativement aux directions principales du 
cristal. Ce cas est traité au chapitre III $ B. 
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et l’équation (2-23) prend la forme : 


1 Dy(z, 0) = Et (0) E (v) Dole, v) 
33 o(Z; = £tr 017 
= Life, v) Dix, v’) ne (2-24) 


Le traitement se poursuit de manière identique. 
L'expression (2-20) reste vraie si on remplace Z 
par Dre 


; fr (v) M{(v) dv 


o — 
2 7 vM{v) 
Sir (0) 


(2-25) 


Équation identique à (1-1). 


III. La longueur de diffusion d’une poudre eris- 
talline. — Nous avons vu précédemment que la 
condition de convergence du développement en 
polynômes de Legendre peut s'exprimer par la 
relation 


w < E (v). 


On sait que dans le cas des poudres cristallines, 
il existe une vitesse des neutrons au-dessous de 
laquelle il ne peut y avoir de chocs élastiques, car 
la condition de Bragg ne peut être remplie. La 
section efficace possède une coupure et tombe à une 
valeur très petite. Il sera donc possible d’avoir des 
valeurs de «, qui, calculées sans précautions, 
seraient égales ou supérieures à Z(v). 

Dans le cas du beryllium par exemple si on 
calcule «w par la formule (2-20) on trouve une 
valeur de 0,457 barns. La section efficace au-dessous 
de la coupure est de l’ordre de 0,46, la condition 
de convergence est mal remplie. Ce phénomène 
résulte du fait que les neutrons «lents », c’est-à-dire 
ceux dont la vitesse est inférieure à la vitesse de 
coupure (v < v.), se comportent différemment des 
neutrons « rapides » (v > v.), et, en particulier, 
ont un libre parcours incomparablement plus 
grand. Il est donc nécessaire d’envisager une loi de 
chocs différente pour les deux groupes de neutrons. 


A. ÉTABLISSEMENT ET RÉSOLUTION D'UNE ÉQUA- 
TION DANS LE CAS DE DEUX GROUPES DE NEUTRONS. 
— Comme il a été suggéré ci-dessus, il est naturel 
de séparer les neutrons en deux groupes. Nous repè- 
rerons les neutrons «lents » (v < v,) par l’indice 
«1 » les neutrons «rapides » (v > v.) par l’indice 
(PT 

Le libre parcours des neutrons lents est déter- 
miné par leur section efficace de capture et par leur 
section efficace inélastique. Nous supposerons, ce 
qui est assez bien vérifié expérimentalement pour 
les neutrons de faibles vitesses (*), que les chocs 
inélastiques font passer les neutrons du groupe 


(*) D. J. Hughes, Pile neutrons research, p. 249. 


lent au groupe rapide. De plus, nous ne restreignons 
pas la généralité en supposant que les neutrons n | 
lents sont issus de neutrons rapides par chocs 
inélastiques. “| 

Les chocs inélastiques peuvent être considérés | 
comme isotropes. Il est donc légitime d’envisager 
la distribution des neutrons « lents » comme 1s0- 
trope. Nous ne particulariserons donc pas la direc- 
tion de la vitesse des neutrons lents, et nous défi- ” 
nirons une section efficace macroscopique de pas- 
sage de la vitesse v, à la vitesse *, 


X(u, Tr) 
telle que : 


Le Xe) der = [7 Xe, ve) dx: = 3R (5-1) 
Ve 


où XF représente la section efficace macroscopique, . 
inélastique totale des neutrons de vitesse m1. 

Ces neutrons ont donc un libre parcours À tel 
que : 


= a = Zi + ZT. (3-2) 


1 
À 
Écrivons maintenant le bilan des neutrons 
« rapides » de vitesse », (en symétrie plane). 


Ur #Ÿ D(a, Tr) _. Ur D D(x, Tr) 


— [ Pie, ex) Dia, et) dot — [ s àn = 0. (3-3) 


Les deux premiers termes sont identiques à ceux 
de l’équation de Boltzmann-Lorentz discutée plus 
haut dans le cas général. Les deux derniers termes 
remplacent le troisième terme de cette équation. 


FIGE 


Les neutrons « rapides » sont séparés des neu- 
trons «lents ». $ indique que l’intégrale doit être 


étendue au spectre des neutrons rapides et le troi- 
sième terme représente les neutrons », nés par 
chocs élastiques et inélastiques de neutrons rapides. 
Le dernier terme en fin donne le nombre de neutrons 
rapides créés par chocs inélastiques de neutrons 


N° & 


lents. S, est un terme de source que nous allons 
calculer ci-dessous. 

Les neutrons lents de vitesse », ont été créés 
eux-mêmes par chocs inélastiques de neutrons 
rapides. En un plan d’abscisse 2’ (fig. 1) il en 
existe un nombre : 


Quiz’) — fe, v) P(v, w) dos. (3-4) 
Dans un angle solide il y a Q(x') à 5 : d0 neu- 


trons de vitesse »,. Ceux qui atteignent le plan x, 
issus de x’ et le traversent sont au nombre de : 


Qi(x’) = d6 * ep do. 


(3-5) 
Au total, le nombre de neutrons issus du plan x’ 

et traversant le plan x est 

aË 


ay fF2sin 0 MT 


Qa&)s J cos 0 


d6 dëe 


I 
| 


9 Æil— œlËl) Qu(a') dE 


co 6 et Élu] = [fete ar 
Le nombre de ces neutrons qui passent à la 
vitesse +, est de : 
X(v1, tr) dé 
ME Je ZX, Cr) dèx 


—7 Æf-alél 


LOS 5: dr) Ei(— o«lË|) Q(x’) dé 


Le nombre total de neutrons « lents » passant 
dans le groupe rapide au plan x sera donc : 


—"4 +0co ; 
LS ho de [ dt 
L De ASS R 
5 Xlntr) El nl) Pir£, &) Da’, v{) de. (3-6) 


L’équation du bilan (3-3) devient : 


Ur LU 2 D(z, Tr) + Ur Z D(x, Tr) 


= ' r) der ed d [e d£ 
. n . = } 
P(Xr, Cr) D, Lr) der : J V] À 


ER J—00 
J Xlo, tx) El olËl) P(t?, w) D(x’, v{) def — 0. (3-7) 


Nous allons supposer maintenant que tous les 
chocs sont isotropes, et nous allons développer D 
en série de polynômes de Legendre. Il est tout à 
fait légitime, dans le groupe «rapide », de se borner 
aux deux premiers termes, car Z(v) CS 

D neTmens que du représente l’élément d’angle 
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solide autour de la direction v, (en symétrie plane) 
[ac Ur) du — Avi, Dr). | 


Un traitement analogue à celui déjà fait donne : 


Ur Ù 
3 dr 


— [, Pivé, v) Dal \D;(æ, Ur) dur + à Fan [4 


A 


-— D 11%, Ur) - Ur 2 D(x, Ur) 


1 


d 1 
3 Ur 5x Dolr, Ur) + 3 Ur > D,(x, Ur) — 0, 
Soit : 
2 
— 55 Dole, 0) + dr E Dale, x) 


— f Pi, w) D o(Z, Ur) dur + +3 fan [Ta 


Î Xe ve) Fit alël) Po, wi) Dole’, vf) dut = 0, 


Nous cherchons pour D o(X, v) une solution de la 
forme 
DVD (x, v) = D(v) ee, 


En remarquant que l’on peut écrire 
P(v', v) = vX(v', »), 
que er" — ee e— 0% et en séparant dans le deu- 


xième terme de (3-10) les sections efficaces de cap- 
ture élastique et inélastique, il vient : 


OR 


Dore (Ur) 


D{vr) + Za (Ur) P(Ur) + D(Ur) 


[x X (Ur, Ur) dur + ®D(uxr) [x Ur, VU] 

—[ x (vr, Ur) D(vr) dur +  J. du 

jee à 
CA 


fa X{vr, 1) Dvr) dur = 0. 


X (01, Ur) Ei(— |Ël) e 


(3-11) 


Nous allons utiliser une méthode analogue à celle 
du paragraphe a) de la deuxième partie. 

Intégrons sur l’ensemble du spectre des neutrons 
rapides. Il vient : 


+ fs 
3 (| 


+ [, ot) Lee Dr) [1— Alu) ZË] dm —0 (3-12) 


2 (Ur) dOr + f. Zalvr) P(Vr) dor 


où 
“| +09 à. e 
AO —; f e—0E El— ulË)) d 
A ÈNE 


et où l’on a tenu compte de la relation (3-1). 
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Calculons 


1 D(05) dr j X(0r, ) (1 — A (o) SË] dm. (3-13) 


Le principe de micro-réversibilité s’écrit : 
ÆX(r, D) Vr M{vr) = X(0, vr) v M(vi) 


où M(v) est la répartition de Maxwell à la vitesse v 
(3-13) devient 


fe du f D{vr) X(v1, Ur) SL D 


On peut supposer que le flux de neutrons rapides 
reste peu différent d’une distribution de Maxwell: 


D{vr) # vrM{(tr). 


[1 — Ai(w) XF] der. 


(3-14) 


(3.14) devient : 


[ mn M{u) [EF — Alu) (EP) ?] dm. (3-15) 


Un rapide calcul montre que 


d Log * + © 


20 Œ] — & 


Alo) = 


2 


Remarquons que l'intégrale donnant Ai(w) n’a 
de sens que si 4 > «. 

Or, 4 — Xi, Cette condition est équivalente 
à la tr déjà trouvée au paragraphe b) de 
la deuxième partie. 

On pourra donc écrire l'équation finale ; en 
tenant compte de 


Mo) = 9 e—t'int 


2 co Dr CO 
3 J Sa) “ie Za (x) D(0r) dur 


(2. sit + o 
+ JL C8 e—2/* | 27 — Log du = 0. 
A $ DS pote Se 1 


(3-16) 


Remarquons que le flux de neutrons lents n’est 
pas nécessairement maxWellien. 
En effet, pour les neutrons de vitesse v, on a : 


D) = 7 Ait alé) dE 


+00 LA 4 + ‘ 
d D(v£) e—0# X(v!, wi) du 
Ÿ 


(9 


=? 
= eos f Fe St ne lé) [7 
— CO ve 


D{vr) X{vr, ) dur 


= e—0% A(w) À Mn). (3-17) 
D'où l’écart à la loi de Maxwell : 
Po) LEE 2 + Zu + © 16 


ni Mie Do SSP 
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Si la capture devient nulle, le deuxième membre … À 


tend vers 1, et 
Du) = à Mn). 


B. CAS D’UNE LOI DE CHOCS NON ISOTROPES, —#] 


Dans ce qui précède, nous n’avons pas envisagé le . 


cas d’une loi de choc non isotrope. Ceci a peu 


d'importance dans le cas des chocs Imélastiques que 


nous supposerons toujours isotropes, mais nous 


avons négligé le cas des chocs élastiques qui dans 
le cas d’une poudre cristalline obéissent à la loi de 
Bragg. 

Remarquons que cela ne signifie pas que la dis- 
tribution soit anisotrope. En effet, une première 
direction donne lieu en un point à des diffusions 
obéissant à la loi de Bragg, donc à des angles bien 
déterminés. Les directions diffusées donnent, à leur 
tour, lieu à des diffusions dont les directions ne 
coïncident pas nécessairement avec une des direc- 
tions incidentes. De plus les chocs inélastiques sont 
isotropes et la distribution des neutrons peut donc 
être considérée comme presque isotrope. Le déve- 
loppement en polynômes de Legendre reste justifié 
si la condition © << Z(v) est remplie. Dans ce qui 
suit nous allons séparer les sections efficaces élas- 
tiques et inélastiques, et ne pas tenir compte de la 
source de neutrons lents qui ne subissent que des 
chocs inélastiques, donc isotropes. 

La loi des chocs élastiques doit obéir à la relation 
de Bragg et on pourra écrire : 


_2f Pate) 


(cos 0 — cos 0x) sin 0 40 do 


: sin 6 dôdo 
(3-19) 

Le facteur 2 est un facteur de normalisation 
nécessaire à l’obtention d’une bonne valeur de la 
section efficace. 0 est l’angle entre v et v', 0w l’angle 
de déviation de Bragg correspondant à la 
réflexion A7. 

L’équation de Boltzmann s’écrit 


0) eŸ De, d) + 0 Zne D(x, d) 


2 [ Palo) S (cos 0 — cos Ow) sin 0 dû de 
LT 


fre 2 Pm(v) S(cos 0/—cos 0x) sin 0’d0/de 
S 


M LT û 


(3-20) 
Dans cette équation X,, désigne les termes non 
élastiques de la section efficace 
Pv',v) si je] £ fe 


S(v’ ,v) = 
eve) 0 si [el = fe 


Nu 


Nous allons développer D en série de polynômes 
de Legendre 


D = D, + 1. 
Il vient 
Du (Dale, v) + Di(e, D] 
+ 0 Zn [Dolz, 0) + Di(x, vi] 
+ S {Due v) + wDilx, v) 


fear) $ (cos 0— cos O) sin 0 d0 do 
&T 


— [540 0) Dole, v') av —$ [Di(x, v') + mDA(e, v°] 


120 à (cos 0— cos 0m) sin 0 d0 do 
AT 


Dans cette dernière somme 


= v” d’où D,(x, v') ES D,{x, v) ; D,(x, v') ex Dfx, v). 


Free 29: 


On voit sur la figure 2 où l’on a pris & comme 
axe des z que : 


u’ = cos à — cos 0 cos à + sin à sin 0 cos © 


portant dans (3-21) on a : 
op À [Do(a, v) + wDi(e, v] 
+ 0 Ze (0) [Dolx, v) + uD;(x, v)] 
+ vDil, 0) Zu (v) — f P(v', v) Die, v') du’ 
+ uD,(x, v) : Py(v) [1 — cos Om] = 0. (3-22) 


En intégrant en y cette"équation telle quelle, et 
multipliée par u, il vient 
Die 0) — —? Data, v){ Eté — SE (M) cos Our ji 
d Dilx, v) 


se dx Ztransport £ 


(3-22) 
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On voit que tous les raisonnements faits ci- 


2 
dessus restent valables en remplaçant Fe D par 


Div = pe S Zu (M, v) cos Om. (3-23) 


En particulier l’équation (2-12) devient : 


2 
=. D ss Dale, v) + v Xa (v) Dil, v) + Dolx, v) 
fre, v') dv’ — [76° v) Dix, v') dv’ = 0 (3-24) 


et l’équation (3-16) 


o? CO ®D(vr) 


: dr + ju Za (Ur) D(vr) dur 
C 


3 ve De (Dr) 
R\2 R $ 
Æ 1 Di° e7 jy /r0° E: te Log _. sa Res e| 
p 20 Z +2—0 
‘du = 0. (3-25) 


IV. Application au cas du béryllium et de la 
glueine. — Nous avons vu que dans le cas d’une 
poudre eristalline, l'équation donnant la longueur 
de diffusion prend la forme (3-25). 

Il est possible à l’aide de cette équation de 
calculer par exemple Z,, supposée en 1/v, si on 
connaît ©, Z4r(0), et XF. 


a) CALCUL DE LA SECTION EFFICACE DE TRANS- 
PORT. — Dans le cas d’une poudre cristalline la 
section efficace élastique est donnée par : 


EN e— M" 


(o} = —  ————— 
TL px À 87 M 


(fig. 3) (4-1) 


— 


GS barns 


15 


5 1e ee _ 
2 3 4 À À 
Fic. 3. — Section efficace calculée de la glucine. 


202 JOURNAL DE PHYSIQUE 


e 2% représente le facteur de Debye-Waller : 


GOT TROIE 
D = ol: + (9) dE 210 

© est la température de Debye du corps consi- 
déré, T la température d’expérience, k, la cons- 
tante de Boltzmann, m la masse de l’atome. 

Fest le facteur de forme par maille élémentaire, 
M une translation du réseau réciproque, # le nom- 
bre d’onde du neutron incident, 4 — 1/x. N'est le 
nombre de mailles par unité de volume. o,, est donc 
la section efficace par maille. 

L’angle de Bragg 0, est tel que : 


sin (0x /2) — M2k. 
Soit 
cos 0m — 1 — (M?]2k?) 
La section efficace de transport prend donc la 
forme : | 


FNM 
Fe Qu 
Str ) mor 2k1 8 ° 
x) Cas du beryllium. — Le beryllium a une 


structure hexagonale compacte ; il s’ensuit que : 
M Apr de paeT Cr” 
A=9/aV3  C—Alfe 


P, 4, r entiers 


où a et c sont les paramètres de la maille. 


2PÉP rq% #] 


F = 20m [1 + cos | 3 


où cp est la. section efficace du beryllium soit 
7,54 barns. La température de Debye est de 
1 000 ©K. On calculera 


1 CO k? el y dk 
ne ZE FN exp — OM) 
ke M<2x 2k48T FÈE 


Après estimation du reste de la série (au 800 
terme) par une intégrale on détermine : 


ER Von 
< 1 JStr > 


(4-2) 


— 5,63 barns 


en bon accord avec la valeur de 5,65 barns calculée 
par R. G. Sachs (*). 

Cette valeur correspond à une longueur de trans- 
port de 1,41 cm. 

B) Cas de la glucine. — La formule (4-1) reste 
valable. Nous avons pris pour simplifier les calculs 
le facteur de Debye-Waller égal à l’unité. 


3 A 
F — [ve L d6 2 Von cos T°] 


[1 + cos ES 


b) 
M? = A'{p? + q? + pq) + C'?r? A’=92/a V3 C'= Ac 
(*) Cette valeur de 5,65 barns est citée par D. J. Hughes 


(Piles Neutrons Research, p. 223). Le calcul de Sachs n’a 
pas été publié. Voir référence [4]. 


No£b| 


avec a! et c' paramètre de la maille de BeO. Le | 
calcul est analogue à celui fait pour le beryllium et ; | 
4 | 


——— — 11,55 barns. La valeur réelle, comptes 
1 JOtr = : | 


tenu du facteur de Debye doit être sensiblement | 


inférieure. La longueur de transport correspon- 
dante est de 1,18 cm. 


b) CALCUL DE LA SECTION EFFICACE DE CAPTURE 


DU BERYLLIUM. — &) La section efficace inélastique 


du beryllium au-dessous de la coupure est-donnée. 
par Palevsky et R. R. Smith et est reproduite dans 
le livre de D. J. Hughes (p. 252). 

On peut admettre que cette section efficace est w 
constante et égale à 0,46 barns pour des longueurs 
d’onde comprises entre 6 À et la longueur d’onde 
de coupure. La section efficace inélastique varie 
ensuite en « 1 /v » pour des longueurs d’onde allant « 
de 6 À à l'infini et vaut ZE — 0,0766 À, ZE en" 
barns, à en À. | 

La dernière intégrale de l’équation (2-25) se « 
calcule aisément si l’on suppose que e—”/* est peu 
différent de l’unité. Si l’on prend pour la longueur 
de diffusion une valeur de 21,8 cm, la résolution 
de (2-25) conduit à une valeur de 


Sy — 14,5 mbarns 


pour la section efficace de capture à 2 200 cm/sec. 

B) Dans le cas de la glucine, L vaut 31,8 cm ce 
qui conduit pour w à une valeur de 0,43 barn. La « 
section efficace inélastique a été mesurée par « 
B. Jacrot et G. Sarma (*). On peut la considérer 
comme variant en « 1 /v »et donnée par l’équation : 


x DT (SF en barns, À en À) 
La résolution de l’équation (2-25) conduit à : 
y — 10,5 mbarns à 2 200 m’s. 


La différence entre les deux valeurs reste impor- 
tante. Elle provient vraisemblablement de deux 
causes complémentaires. 

1) Dans le calcul de Af, il n’a pas été tenu compte 
du facteur de Debye, difficile à exprimer dans le 
cas de mailles formées d’atomes différents, ce qui 
conduit à une valeur trop faible de 4 et par consé- 
quent de X*. 

2) Les corps utilisés pour les mesures ne doivent 
pas avoir la même pureté ; à la plus grande pureté 
correspond la plus faible section efficace de capture. 


€) DÉFORMATION DU SPECTRE DE MAXWELL. — 
L'écart à la loi de Maxwell est donné par la for- 
mule (3-18). Cette formule donne la valeur limite 
du rapport du spectre déformé au spectre de 
Maxwell. Cette limite existe comme conséquence du 


(*) Communication privée. 
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F1G. 4. — Déformation du spectre de Maxwell dans le 
béryllium. 


flux des neutrons, 
— densité des neutrons, 
— .— spectre de Maxwell. 


sait qu’il a été possible de définir une longueur de 
diffusion. La partie déformée reste solidaire du 
reste du spectre. D. J. Hughes nous avait signalé 
oralement que cette déformation avait été observée 
expérimentalement dans le cas du beryllium. Nous 


x) CAS DU BÉ 
x en À 3,95 5 6 7 8 
Dr[Dm 1,42 1,42 1,42 1,25 4,17 
B) Cas 
à en À 4,68 5 6 7 8 
DD 2,42 1,73 1,32 1,20 4,14 
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F1G. 5. — Déformation du spectre de Maxwell dans la 
glucine. 


flux des neutrons, 

— densité des neutrons, 

— .— spectre de Maxwell. 
D 


avons fait le calcul de la formule (3-18) pour un 
certain nombre de longueurs d’onde supérieures à 
celles de la coupure. Les résultats pour le béryllium 
et la glucine sont reportés sur les courbes ci-jointes 
(fig. 4 et 5) et dans le tableau ci-dessous. 


RYLLIUM 
9 10 41 12 13 1% 19 
1:12 1,10 1,08 1,06 1,05 1,0% 1,0% 
DE LA GLUCINE 
9 10 ddl 12 13 14 15 
1,10 1,08 1,07 1,053 1,045 4,04 1,03 
Manuscrit reçu le 22 octobre 1960. 

GRAPHIE 


[8] Fermi, Strurm et Sacns, Phys. Rev., 1947, 71, 589. 
[4] Sacs (R. G.), Rapport déclassé : A. E. CG. D. 2269. 
[5] HorowiTz (J.) et Messrax (A.), Rapport SPM n° 5. 


LE JOURNAL DE PHYSIQUE ET LE RADIUM TOME 22, AVRIL 1961, PAGE 204. 


) 


VALIDITÉ DE LA REPRÉSENTATION PAR MODÈLE ÉLECTRIQUE 
DES EFFETS DE RÉSONANCE MAGNÉTIQUE NUCLÉAIRE 


Par GEorGes BONNET, 
Centre d'Études Nucléaires de Grenoble. 


Résumé. — L'étude du comportement d’un transducteur à résonance magnétique formé de 
l’association d’un réseau électrique et d’un ensemble de noyaux dotés de spin, montre qu’il est 
possible d’en déduire une représentation analytiquement équivalente au moyen d’un modèle 
entièrement électrique utilisable pour un régime transitoire aussi bien que pour un régime perma- 
nent. Une étude détaillée des conditions de validité permet d’en justifier l'emploi dans la majorité 
des cas. 

On propose enfin un critère de linéarité des équations de Bloch en régime transitoire, qui cons- 
titue un prolongement de la condition connue de non-saturation en régime stationnaire. 


Abstract. — When studying the behaviour of a magnetic resonance transductor formed by 
the association of an electrical network and of a set of nuclear spins, it is possible to bring about a 
representation that is analytically equivalent by means of an entirely electrical model, availa- 
ble for transients as well as steady-state. À detailed study of the validity conditions justifies 
its use in most cases. : 

Also proposed is a linearity criterion of Bloch’s equations in transient-state that is simply 


the prolongation of the well-known condition of non-saturation in the steady-state. 


I. INTRODUCTION. 


L'emploi d’un réseau électrique analytiquement 
équivalent à un système physique complexe a 
trouvé sa source dans le domaine de lélectro- 
mécanique et les travaux de P. Le Corbeiller, 
Y. Rocard, W. P. Mason — parmi tant d’autres — 
lui ont assuré un grand développement. 

Il semble que l’intérêt principal de cette méthode 
réside dans l’unicité axiomatique qu’elle confère 
aux problèmes d'adaptation réciproque entre sys- 
tèmes d’essence différente. En outre, la théorie des 
réseaux électriques a fait l’objet d’une somme de 
recherches considérable et l'application de ses résul- 
tats à d’autres domaines qui lui sont intégrés par 
cette"voie est toujours bénéfique. 

N. Bloembergen et R. V. Pound [1] avaient fait 
allusion à la similitude entre un montage à réso- 
nance gyromagnétique et un réseau électrique à 
double accord, et tout récemment P. D. Carter a 
proposé un schéma équivalent à un filtre à cristal 
de grenat [2], ces auteurs n’envisageant d’ailleurs 
que l’équivalence en régime stationnaire idéal. 
G. Brunneret A. Loesche{[3]ont proposé un modèle 
qui rend compte du comportement d’un système de 
spins isolé et l’ont étendu du régime stationnaire à 
celui de « passage rapide ». [is n’ont pas envisagé 
la recherche d’un modèle plus complet qui puisse 
traduire l’interaction entre spins et circuit élec- 
trique de mesure. 

Nous nous proposons ici de rechercher le modèle 
électrique généralisé des transducteurs de champ 
magnétique qui utilisent la résonance nucléaire, en 


mettant l’accent sur ses conditions de validité dans 
un régime quelconque. 

L’observable utilisée dans de tels transducteurs 
est la fréquence des transitions entre sous-niveaux 
Zeeman d’un ensemble de moments nucléaires 


placés en présence du champ magnétique semi-per- 
manent à étudier et soumis à un rayonnement trans- 
versal de radiofréquence, la relat'on champ/fré- 
quence étant celle de Larmor. On peut schémat'ser 
le. montage de base sous la forme d’un échantillon 
gyromagnétique contenu dans le volume intérieur 
d’un circuit électrique générateur du rayonnement 
transversal, lui-même fermé sur une certaine charge 
électrique. L’axe u du générateur du thamp magné- 


N° 4 
tique de radiofréquence — bobine ou cavité suivant 


— 

le cas — et la direction du champ permanent, 

supposés orthogonaux pour simplifier, définissent 

une base fixe (e, e, e;) que l’on norme, et l’on a, 
ES 


par exemple : u = jule,,X, — H, e, (fig. 1). L’en- 
semble des spins nucléaires d’une part, et le réseau 
électrique linéaire d’autre part, constituent un sys- 
tème couplé dont nous étudierons successivement 
l’évolution libre et forcée. 


II. ÉVOLUTION LIBRE DU SYSTÈME COUPLÉ. 


A. Ensemble gyromagnétique isolé. — Si l’on 
considère le système de spins isolé et soumis uni- 
quement à l'influence du champ permanent, on sait 
que, sous certaines conditions, son évolution, envi- 

| sagée sous l’aspect macroscopique, est régie avec 
_ une très bonne approximation par le schème de 
Bloch [4]. Dans la base définie précédemment, l’ai- 
mantation M{(+) est reliée à l’aimantation à l’équi- 
libre m, e3 par : 


aM 1 
de —= (M VAN H;) LT (M.e;) €; 
L'iMe,)e,—L[(M.e.) e, (1) 
WT, -C2 2— 7, | C3) — Mo] €3 


où y est le rapport gyromagnétique, caractéristique 
des spins nucléaires utilisés, 7, le temps de rela- 
xation spin-milieu et 7’, le temps de relaxation spin- 
spin. Si l’on s'intéresse à la paire constituée par 
lPaimantation projetée sur l’axe de la bobine, 
m, — (M.e,), et l’aimantation projetée sur le 
champ permanent, m; — (M.e;), il est facile de 
tirer de (1) le groupe d’équations : 

2 9 

tra t (oi +73) | = 0 

Pari = 7 Mo 


Nous avons supposé implicitement que la condi- 

1447, L 
Hdi T; 
remplie, comme on le verra plus loin. L’équa- 
tion (2a) montre que le système est pendulaire et 
que la solution libre est une oscillation de fréquence 
cyclique w, — yA,, fréquence de Larmor, amortie 
suivant le décrément à — 2r/w, T On peut 
parler d’une «surtension des spins » égale à 


Qs = r[è = (1/2) Go Tr 


tion 


se trouve toujours largement 


valeur qui, en pratique atteint et dépasse souvent 
104 (*). La composante m, évolue apériodiquement 


(*) C’est en tenant compte de l’aspect pendulaire du sys- 
tème gyromagnétique isolé que Brunner et Loesche [loc. cut] 
ont utilisé expérimentalement un circuit résonant élec- 
trique qui leur a permis une très belle simulation des 
phénomènes de « wiggles » en passage rapide. 
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avec la constante de temps 7, (m, suivrait une loi 
identique à celle de m,). Il est remarquable que le 
système isolé soit un système linéaire comme le 
montre (2) et nous allons voir que c’est le couplage 
qui introduit une non-linéarité. 


B. Réseau électrique linéaire. — Nous allons 
supposer d’abord que la charge de la bobine — 
dont Z est l’auto-inductance et À la résistance de 
pertes — est passive et constituée par une simple 
capacité C. Ce réseau constitue un deuxième sys- 
tème pendulaire dont l’évolution libre est régie par 
une équation homogène du type : 


de Te + = 0 (8) 
or est la population propre : or — (LC)-V2, 

Q est la surtension : 0 = Lor/R. 

L’interaction avec le système de spins introduira 
un terme de couplage au second membre de (3). Si 
l’on tient compte de ce que la surtension du réseau, 
de l’ordre de la centaine ou même moins, est très 
fable devant la « surtension des spins » il est à 
prévoir que cette interaction imposera à (3) un 
régime forcé dont le spectre de fréquences sera 
centré au voisinage de la fréquence de Larmor, ce 
qui amène à accorder le réseau près de cette fré- 
quence : © æ y. 

À cause de l’effet de filtrage, ce spectre est néces- 
sairement étroit, sa largeur étant de l’ordre de la 
bande passante du réseau tout au plus. Il en 
résulte que le régime sous interaction peut se 
mettre sous la forme 


Alt) cos ot + Alt) sin ot avec © = &. (4) 


Les amplitudes A,(#) et A,(t) varient lentement 
dans le temps et, en particulier, leur variation 
pendant une période est très faible, d’où 


d 
2e Ai & 0 Ai. (5) 


C. Système couplé. — L’interaction entre l’en-. 
semble de spins et le réseau électrique se présente 
sous un double aspect : le courant circulant dans la 
bobine crée un champ qui s’ajoute vectoriellement 
au champ permanent et perturbe le régime libre, 
les variations de la projection de l’aimantation sur 
l’axe de la bobine se traduisent d’autre part par une 
variation de son flux d’induction et y induisent une 
force électromotrice. 


1. RELATIONS COURANT-AIMANTATION, — Consi- 
dérons d’abord la première composante du courant 
conforme à la forme (4) : 


= )icos at. (6) 


— 
Sachant que le champ magnétique & (r, t) qu'il 
produit lui est proportionnel en tout point r de 
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l’espace, Fe (r,t) = k(r).i(t), nous pouvons suivant 
I. Solomon [5], [6], définir un champ magnétique 
moyen H dans le volume + de l’échantillon : 


ba He FE {ke ) ar] it 


lequel, si l’on tient compte de la symétrie de la 
bobine, est dirigé suivant son axe : 


H — He, — KI{t) cos wt.e:. l (7) 


= il) 


Il est bien connu que ce champ oscillant peut se 
décomposer en deux vibrations circulaires tournant 
dans deux sens opposés dans le plan [001] et que, 
selon F. Bloch et A. Siegert [7], l’une de ces compo- 
santes n’a d’autre effet que de déplacer la fréquence 
de résonance des spins de la quantité 

(AH o)° (o012T) 
correction qui, dans les conditions expérimentales 
habituelles où 4 [H, < 10, est entièrement négli- 
geable. L’autre composante circulaire s’écrit dans 
la base fixe et si, par exemple, y est positif : 


H* = 5 KI() cos bte, —; Kl(tsin te, (8) 


et nous savons, d’après (5), que son module ne 
variera pas de façon appréciable pendant quelques 
cycles de rotation. Le système d’équations de Bloch 
s’écrit alors 


d 1 4 : 

& + Mi — Optez | 9 YÆI(t) sin wt.ms; (9a) 
d 4l 4 
d TT Ma = — O0 M + ; yÆXI(+) cos ot.m3 (9b) 
d 


1 1 ; 
a, M3 ose) sin of.m: 


_ ; yKI(E cos wt.m. (9c) 


Par comparaison avec les équations (2) de l’en- 
sembleïsolé, on voit que ce sont les termes de cou- 
plage qu‘ détru'sent la l'néar'té ; m; dépend de (+) 
par (9c) ce qui introduit des termes non-linéaires 
dans les seconds membres de (9a) et (9b). Cepen- 
dant, si nous nous limitons à un champ AH suffi- 
samment faible pour que sa composante circulaire 
puisse être traitée comme une perturbation, il 
apparaît [8], qu’à une approximation du second 
ordre près en ?, m, suivra sa loi d'évolution libre 
définie par (2b). Les équations (9) se trouvent ainsi 
linéarisées et nous obtenons après élimination de 
LEE 


d? 2 ol 
tant (oi + … 


—. yK (re \m3 [Lo + O&09) COS œt + ch sin ot] 
co Te 


ma 


Quelques remarques amèneront à simplifier le …| 
terme perturbateur du second membre : F. 


a) nous avons vu que w demeure voisin de wp "| 


et nous pouvons remplacer (© + &,) par 20 ; 
b)"w, étant très grand devant 1/7T,, on peut 
négliger le terme en (1/T,) sin ot ; 
c) d’après la loi d’évolution (2b) que nous admet- : 
tons dans l’approximation linéaire, FLE est de 
l’ordre de (m3 — m3) [T;, et par suite : 


t 
: My LE O3 My 


L'évolution des moments nucléaires est alors 
régie par la loi : 


d? 2 de 
dt Td c% mi 
ne fe cos ot + sh ot] (10a) 


MOT {em 


2. RELATIONS AIMANTATION-TENSION. — Les mi; 
représentent évidemment une valeur moyenne dans 
l'échantillon des composantes de l’aimantation. 
Suivant les équations de Bloch linéarisées, la valeur 
locale de la projection de l’aimantation sur le champ 
est proportionnelle en tous points à ce dernier. St À 
est ce facteur de proportionnalité à l'instant f, nous 
pouvons en tenant compte des propriétés géomé- 
triques de symétrie et de l’équation (7), mettre m, 
sous la forme 


ee xi(t) [ kr) dr = ÆXi(). 


On en déduit la valeur de l’énergie emmagasinée 
par les spins 


HER x [ dr=—{UJÆ) mil fe dr 
T UT 
et la force électromotrice induite dans la bobine 


1 2 
— x [# dr.dm,fdt. 


On a coutume [6] d'utiliser le concept de « coef- 
ficient de remplissage » défini par 


n= fraff k? dr. 
ï PJes 


Nous adopterons un système d’unités non ratio- 
nalisées. L’intégrale située au dénominateur, éten- 
due dans tout l’espace, est alors ‘égale à 4rL/u,, 
L étant l’auto-inductance de la bobine et u, la 
perméabilité, du ‘vide. Le coefficient de remplis- 
sage représente ainsi la fraction d’énergie magné- 
tique que renferme le volume de l’échantillon en 
l’absence de matière aimantée. Finalement l’équa- 


(11) 


N° 4 


tion différentielle du réseau électrique prend la 
forme 
or d 


d? : 
droit oRli=— 


kr d°m: 
DAT de 


(12) 


D. Le modèle électrique. — Nous sommes main- 
tenant en mesure de définir le modèle électrique 
dont le comportement puisse se représenter du 
point de vue formel par les relations (10) et (12). 
Il y a pour cela une condition nécessaire : 
Y. Rocard [9] a montré que, dans le cas de sys- 
tèmes couplés dont les paramètres dérivent d’un 
potentiel, la conséquence de l’existence des éner- 
gies potentielle et cinétique est que les termes de 
couplage sont nécessairement symétriques. En l’es- 
pèce, les seconds membres de (10) et (12) ne peuvent 
_ être tels que sous la condition 

Le © [or cos œot + LÉ sin et] 
e étant un scalaire indépendant du temps, et 
comme z(t) a la forme (6), on en déduit que la 
représentation par un modèle électrique n’est pos- 
sible que si p — — w, donc p # — 9, et si 


di Jdt € wo JL. 

C’est cette dernière loi de variation qui, dans le 
cas d’une charge capacitive, se trouve exprimée 
ipso-facto comme nous l’avons vu dans l’équa- 
tion (5). En introduisant la susceptibilité parama- 
gnétique statique Yo = Molto/, l'équation (10a) 
prend alors la forme (1) : 


d 24 1 
D nee 2 RES 
rage 


On obtiendrait un résultat identique avec la deu- 
xième composante du courant, en sin of, que sug- 
gère la forme (4) et la linéarité des équations en 
conserve également la validité pour le courant 
total. 

L'ensemble des relations (12) et (13) permet de 
bâtir entièrement le réseau équivalent. Il comprend 
deux mailles couplées par une mutuelle induc- 
tance (fig. 2), la première d’entre elles, parcourue 
par le courant 1, — {(t), étant formée évidemment 
par une inductance ZL, et une résistance À, iden- 
tiques à celles, L et À, de la bobine, placées en 
série avec la capacité C, = C = oœg?Lr1 Il 
existe une certaine souplesse dans la représen- 
tation de la maille 2, puisque les seuls paramètres 
déterminants sont sa fréquence propre et sa sur- 
tension, et les valeurs des éléments résultent en 


(1) On a choisi a dessein une grandeur sans dimensions, 
Mo H5 
Uo H3/8T 


car 8TYXo — 
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fait de la convention de correspondance adoptée 
pour son courant #,. Nous en profitons pour poser 
arbitrairement L, — L, — Let obtenons ainsi 


Inductance L, = L Capacité C, ? = L[oÿ + (4/T3)] 


Résistance R;, = 2L/T;, 


Courantirs _. Varn.MolXo M3 


Inductance mutuelle : A — L /4rn%o Mslmo 


ou coefficient de couplage : 


z = HJVL, Li = Vanne MslmMo 


(15) 


DDDD0DDD 
ra 
Re 
ee 


a 


Fi. 2. — Modèle équivalent en évolution libre. 


Dans la figure 2, la traduction champ/fréquence 
s’opère par les variations de la capacité C, avec la 
IE 


1103 = (Ho) = LULITH + YA. (6) 


; III. GÉNÉRALISATION. É 
ÉVOLUTION FORCÉE DU SYSTÈME COUPLÉ. 


A. Évolution forcée. — Nous avons utilisé 
jusqu’à présent une hypothèse assez contraignante, 
selon laquelle la charge de la bobine était cons- 
tituée par un unique condensateur, et nous nous 
proposons d’étendre l’étude au cas le plus général 
d’une charge quelconque, représentée par un opéra- 
teur fonctionnel linéaire Z5,, qui peut inclure en 
outre une source de tension e(t) (ou de courant) 
donnant lieu à un régime forcé. Le spectre de 
Fourier de la source peut comporter une ou plu- 
sieurs raies ainsi qu'une partie continue, et nous 
introduisons cette unique restriction que, pour que 
des transitions puissent être induites, ce spectre 
devra posséder une densité notable dans le voisi- 
nage de la fréquence de Larmor. 

Dans le réseau ainsi défini, le courant qui tra- 
verse la bobine est soumis à la relation 


(17) 


d , 
[1 Fe Zas | reel) 


(2) Bien que 1/73 << y? HE nous avons conservé la 
valeur exacte en vue d’éventuels calculs de précision. 
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et si sa densité spectrale est définie par 


Ile) = (2x) [© 


le champ moyen qu’il crée dans le volume de 
l’échantillon prend la forme 


HA 
— CO 

Dans la bande de fréquences cycliques infiniment 

étroite centrée sur « et de largeur dw,il correspond 

à la composante élémentaire dH(«) de ce champ 

rectiligne une vibration circulaire plane de sens 

convenable dont il est facile de trouver l’expres- 
sion, soit pour y positif (voir Appendice) : 


i(t) e—iot de (18) 


) eiot do.e. (19) 


dH*(o) = ; lo) ciot Le: Æ ire ea] du, avec = —1 
(20) 


et, si nous supposons remplies les conditions de 
linéarité des équations de Bloch, nous sommes 
conduits, en vertu de la règle de superposition, à 
écrire ces équations sous la forme 


Rte . Mi = Op Mao 
+ 
_- vins | it I(w) eiot dw (214) 
d 1 me 
ET Ds our +5 vEms [ en) eiot do 
(21) 
d 1 m 
Fr + _ M3 T,- (21c) 
On en tire par élimination de m 
dé el 1 L +oo : 
ie an + (oi Eye) | = VE [ET Ho) oies 
NC 1 n) 7.)| : 
[ratlol I © 0) Tr L lol (rs Re Fu 7. do (22) 


Remarquons tout d’abord que le terme imagi- 
naire situé entre les crochets du second membre 
est de l’ordre de m,/T,, donc très inférieur à w, Ma, 
ce qui permet de le négliger (*). Si, en première 
approximation, on négligeait la lente variation 
temporelle de m, et en tenant compte de ce que 
1/75 € w$, on pourrait exprimer la solution de (22) 
sous la forme du terme spectral de m, : 


[ol + ol 
@2 — ©? + 120/T; 


M,{) #5 yAm; (23) 


et la forme du dénominateur nous montre que ce 
terme n’est vraiment important que dans une 


(*) Cette Ross serait fausse dans le cas où 
ms < Moloo Le, (mis à part le cas de saturation totale qui 
ne présente aucun intérêt) ce qui correspondrait à un désali- 
gnement du moment par rapport à sa position d'équilibre 
égal à 90° à quelques minutes d’arc. Ce dernier état est peu 
probable et en tous cas de très faible durée de vie. 
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ne très étroite comprise entre &ç — (1/T2) et 


+ ({/T d). Dans ces conditions, cet effet de fil h 

a des spins nous autorise à remplacer dans (22) … 

[og + |w|] par 26?/o, ce qui n’affectera aucune-. 

ment l’argument de H,(«) et introduira une erreur n 

infime sur son module. Comme d’autre part 
+00 d2? 

Le &w? I(o) eiot do = — 0 
l'équation (22) se transforme finalement en une. 
équation identique à (13). Les résultats acquis dans 
le cadre de l’évolution libre s’appliquent donc à | 
l’évolution forcée, et, tenant compte de (17) et (21c). 
nous adoptons comme schéma équivalent au sys- 
tème formé par une bobine et un échantillon gyro- 


R, 
Z 5e 36 
FL Vie £ 
Ë Dre c 

© 
nie, 

ès NE 

zL 

F1c. 3. — Modèle équivalent dans le cas général. 


magnétique celui de la figure 3, dont les éléments 
sont déjà définis par (14) et (15). 

Ce modèle est valable aussi bien en régime per- 
manent que transitoire. Il offre ceci de particulier 
que le coefficient de couplage (15) est une variable 
susceptible de dépendre du temps par lintermé- 
diaire du rapport [m3/m0]"?. 


1. RÉGIME PERMANENT. — Le coefficient de cou- 
plage est constant d’après (21c) et lorsque le facteur . 
de saturation est faible, m,/m, # 1. Si l’on veut 
améliorer l’approximation, 1l est facile d'adopter la 
solution exacte des équations de Bloch (cf. M. Sou- 
tif [13]) soit pour {(t) — Z, cos (ot + o): 


CE SES 4 + To — w,)? 
m q| > 2 r ù 
* FEVER Ti Ts 


1 + Ti(o — wo)? 


Seuls sont touchés par cette correction le coef- 
ficient de couplage et le coefficient de proportion- 
nalité entre courant et composante axiale de l’ai- 
mantation. 


2. RÉGIME TRANSITOIRE. — Dans certains cas 
seulement, si la valeur initiale m$ de m, au temps 0 
diffère par trop de sa valeur d’équilibre m,, le 
coefficient de couplage varie avec le temps. 

La loi suivant laquelle m, tend vers m, est 


Mmalmo = 1 + [(m$[mo) — 11 exp 

1 
et comme T, > «51, le couplage demeurera prati- 
quement constant pendant un grand nombre de 
périodes, ce qui permet de traiter de façon simple 


N°& 


le régime transitoire du réseau à paramètres 


variables. 


B. Cas d’une polarisation négative. — Plusieurs 
procédés ont été proposés qui fournissent un ren- 
versement x des populations de spins nucléaires, 
soit par prépolarisation (H. Benoît et J. Henne- 
quin [10})), soit par polarisation dynamique 
(A. Abragam et all. [11]) (*). Comme ces procédés 
s’accompagnent d’une très forte augmentation de 
lamplitude de la polarisation, ils s'avèrent d’un 
grand intérêt dans la réalisation de transducteurs 
de type « MASER » en champs faibles. 

L’aimantation à l’équilibre dans un champ per- 
manent y prend une direction opposée à celle du 
champ, — mes, et comme l’évolution macro- 
scopique en présence d’une polarisation quasi-sta- 
tionnaire est toujours régie par le schème de Bloch, 
la seule modification à introduire dans les calculs 
précédents est d’y considérer la susceptibilité para- 
magnétique statique y, comme une valeur négative. 
Dans le schéma de la figure 3, le coefficient de 
couplage devient alors une grandeur imaginaire, 
ceci d’après (15) : 


g = MVL, L, = iVärnlxolmslmo. (15 bis) 
Les éléments du circuit demeurent réels, mais la 
relation de correspondance courant/aimantation, 
d’ailleurs conventionnelle, est imaginaire. 
Ces résultats paraissent logiques, si l’on consi- 
dère la nécessité de rendre compte de la fourniture 
d'énergie qu’opère dans ce cas le système de spins. 


C. Possibilités d’extension. — Il est relativement 
aisé d'étendre le schéma équivalent d’un réseau 
électrique en interaction avec un système de spins 
à des montages plus complexes mettant en œuvre 
plusieurs bobinages, dans le cas où leurs axes 
demeurent parallèles. 

Il n’en est pas de même dans le cas contraire, et 
en particulier, celui du montage de Bloch à bobines 
croisées. L'étude du groupe d’équations différen- 
tielles, auxquelles conduit un calcul analogue au 
précédent, montre alors l'existence de termes de 
couplages de signes opposés, du type de ceux que 
Y. Rocard nomme «termes gyroscopiques » [9]. On 
sait qu’il est alors impossible de rendre compte du 
comportement de l’ensemble, dans les conditions 
les plus générales, par un modèle purement élec- 
trique. Si l’on se limite aux seuls régimes station- 
naires, on peut par contre bâtir un tel modèle en y 
introduisant un simple élément déphaseur, dont la 
présence correspond au décalage temporel des com- 
posantes, suivant chaque bobine, de l’aimantation 
macroscopique lors de son mouvement de 
précession. 


(*) On trouvera une excellente étude d’ensemble de ces 
procédés dans un récent article de P. Grivet [12]. 
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IV. CONDITIONS DE VALIDITÉ. 


A. Les conditions de validité sont tout d’abord 
celles même que Bloch a introduites pour poser ses 
équations [loc. cit], à savoir : 

— On supposera que l’équilibre thermique entre 
spins et milieu est établi dans le champ perma- 
nent 1, au moment où l’on applique le champ 
transversal. Cette condition n'implique prati- 
quement qu’une restriction très limitée dans l’em- 
ploi du modèle électrique appliqué à un trans- 
ducteur. 

— Homogénéité du champ permanent : impli- 
citement, la valeur adoptée pour 7, peut tenir 
compte de l’inhomogénie, à condition qu’elle soit 
faible et contribue à un élargissement modéré de la 
raie naturelle. Autrement dit, dans le cas des 
liquides, le champ moléculaire devra constituer le 
processus de relaxation dominant, de façon à ce 
que la décroissance des composantes transversales 
de l’aimantation conserve un caractère exponentiel, 
celui qu'implique le réseau électrique dissipatif 
employé dans la représentation. Sinon le modèle ne 
pourra convenir à tous les cas d’expérience et il 
serait nécessaire de le compléter en introduisant 
une répartition de fréquences propres sur un 
ensemble de mailles de spin, de préférence continu, 
ce qui ne présente plus aucun intérêt pratique. 

B. Nous avons rencontré également au cours de 
ce calcul d’autres restrictions qu’il reste à étudier. 
Tout d’abord nous avons posé, à propos de l’équa- 
tion (2), que la stabilité à court terme du champ 
«permanent » devait être telle que 


1H 
H, dt 


4! 
HE 


Nous pouvons considérer que ce champ comporte 
une partie fixe < H, > à laquelle se superpose 
une partie fluctuante, de nature plus ou moins 
aléatoire e(t) beaucoup plus petite que la première. 


On voit alors que la condition précédente, qui 
an ee A re ét tnE 
{ e dt Ta Fe 


très largement remplie, et d’ailleurs la mesure 
de A, perdrait autrement toute signification. 

C. En outre, nous avons admis la possibilité de 
remplacer, dans l’équation (22) la quantité 

CIRE par 2&?/«$. 
0 

Cette approximation implique une erreur de repré- 
sentation si « s’écarte trop de la fréquence de 
Larmor &,. Deux cas sont à étudier, suivant que 
le spectre du courant est diseret ou continu : 


1. RÉGIME STATIONNAIRE. — Supposons que le 
spectre soit constitué par une raie unique centrée 
sur + 3, avec wi y. 


1% 
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On sait [13] que dans le cas d’une perturbation 
périodique, la solution exacte des équations de 
Bloch conduit à décrire la variation d’impédance 
de la bobine entourant l’échantillon par la relation : 


1 


= Ts PAGE CAT PeRtIEse 
NP = TX Oo Le A “Lo RESTE 


En utilisant le cireuit équivalent de la figure 3, 
cette variation d’impédance apparaît comme due 
à l’impédance couplée (zLw:)?/Z22 qui s'écrit, sui- 
vant les équations (14) et (15) : 


AZ = 27NY0 Oo RAR Lo, 
mo 
©1160 


5 me ‘EA 
1 — 1 2 [oi — ot + x] 
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On voit, par comparaison, que l’approximation 
conduit à trois sortes d'erreurs : 

a) Une erreur sur le centrage de la résonance, 
Sol &0 # (1/2) (wo T2) ?. Elle sera le plus souvent 
très inférieure à 1078, donc entièrement négligeable. 

b) Une erreur sur la forme de la raie d’absorp- 
tion, qui n’est plus exactement lorentzienne. Sa 
valeur relative, de l’ordre de [w, — w,]?/80$, est 
inappréciable pour des désaccords atteignant déjà 
plusieurs dizaines de largeurs de raie. Elle n’est 
sensible que sur la « queue » de la raie, mais il faut 
remarquer que c’est également dans ce domaine 
qu’est discutable la forme même de Lorentz à 
laquelle conduisent les équations de Bloch. On sait 
(cf. Y. Ayant [14], [15]) que, même dans le cas de 
liquides à faible viscosité, l’influence du temps de 
corrélation de Debye est susceptible de se mani- 
fester pour les grandes valeurs de |w, — «il, et 
il est plus vraisemblable que la raie d’absorption 
ait une forme générale de Lorentz avec des queues 
gaussiennes. 

Le modèle électrique peut donc être considéré 
comme satisfaisant à cet égard pour l’étude du 
fonctionnement en transducteur, la nécessité dans 
ce cas d’une énergie d’interaction appréciable entre 
le champ électromagnétique et le système de spins 
impliquant en effet une fréquence de travail telle 
que T,|o, — w,| soit de l’ordre de l’unité. 

c) Une erreur sur l’amplitude, due au facteur 
à, /w, que nous négligerons pour les mêmes raisons. 


2. RÉGIME TRANSITOIRE. — Nous pouvons dé- 
composer le spectre de fréquences du courant 7, en 
trois parties principales : 

a) 0 < [wo] < 1/T,;. L'évolution très lente qui 
correspond à cette bande conserve l’équilibre ther- 
mique entre spins et milieu et se traduit par un 
vecteur d’aimantation macroscopique parallèle et 
proportionnel au champ magnétique résultant. Il 
y a donc rotation et variation de norme de cette 
aimantation. L’effet en est négligeable si le champ 
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transversal est faible, sinon le modèle électrique, « 
qui ne peut en rendre compte, devra faire place au 
calcul direct. , à 

b) 1/7, < [wo] < «4. La perturbation intro-. 
duite par cette partie du spectre répond au critère 
de validité du théorème adiabatique. La répar- … 
tition des populations entre sous-niveaux Zeeman 
est donc sensiblement la même qu’en l’absence de 
perturbation, mais l’axe de quantification évolue » 
parallèlement au champ résultant, ce qui traduit 
par une rotation du vecteur aimantation macro- 
scopique sans variation de sa norme. Si la densité. 
spectrale est du même ordre dans cette bande 
qu’autour de w,, il correspond à Æ(«) ef” de une 
variation 

dm À Uo Xo H(o).do 


entièrement négligeable par rapport à ce que four- 
nit le voisinage de wç correspondant 


1 
dm, & 3 Vo Xo So T: H(ow).do. 


Mais si la majeure partie de l’énergie est répartie 
dans la bande adiabatique, le résultat ci-dessus 
cesse d’être valable et le phénomène de rotation . 
dominera celui de résonance. D'ailleurs, à moins de » 
conditions expérimentales défavorables, cette 
partie du spectre proviendra d’événements séparés 
dans le temps de ceux susceptibles d’induire des 
transitions. On pourra alors traiter les rotations 
adiabatiques par le calcul direct et le modèle élec- 
trique sera utilisé à l’étude des variations non- 
adiabatiques, en prenant pour conditions initiales 
celles fournies par ce calcul. 

c) [of + wget |[w| > «,. Cette troisième bande 
entre dans le cadre des « modifications soudaines » 
de l’état dynamique du système de spins. C’est la 
seule qui intervienne dans l’interaction avec le 
champ électromagnétique dont tient compte le 
schéma électrique étudié. On a vu par l’équa- 
tion (23) l’extrême étroitesse de la bande de fré- 
quences qui intéresse les transitions entre niveaux 
énergétiques. De ce fait, et plus encore que dans le 


. « (a) (a) 
cas stationnaire, le remplacement de [er par 


L . 0. 
2w2/w$ ne peut avoir qu’un effet inappréciable et 
se justifie pleinement. 


D. Linéarité. — L'emploi d’un réseau électrique 
de caractère linéaire pour représenter des phéno- 
mènes de nature physique différente ne peut se 
justifier que si ces phénomènes eux-mêmes sont 
régis par des lois linéaires, du moins avec une 
bonne approximation. 

Or, nous avons vu que la linéarité des équations 
de Bloch est détruite par les termes de couplage, 
dès lors que la composante de l’aimantation sui-. 
vant le champ permanent dépend de la valeur du 
champ transversal. 


| No 4 


4. Dans le cas de régimes stationnaires, le critère 
de linéarité est bien connu [4]. Nous l’écrirons : 


(HT, T, <1 


H — K.1, étant l'amplitude du champ rectiligne 
sinusoidal créé par le courant Z, cos «t qui cireule 
dans la bobine. 

2. La définition d’un tel critère est beaucoup 
plus délicate lorsqu'il s’agit de perturbations tran- 
sitoires, et nous allons essayer de l’aborder ici : 

Considérons donc une modification soudaine de 
Phamiltonien, mise en route au temps & — 0 et 
provoquée par un champ transversal issu d’un 
courant 1,(t) dont nous supposons la densité spec- 
trale notable autour de &,, de telle façon que les 
probabilités de transitions soient importantes. 

Sauf dans un cas unique (oscillation très peu 
amortie qui, d’ailleurs, si sa largeur spectrale est 
inférieure à 1/T,, pourra se ranger dans le cadre 
_ d’une perturbation périodique) le courant z, ainsi 
_ défini est tel que son spectre est nécessairement 
très étalé relativement à la largeur de raie. 

L'aspect de l’équation (23) nous a montré que 
l’effet de filtrage des spins est tel que les compo- 
santes spectrales situées en dehors de la raie 
_ d’absorption ne jouent aucun rôle. On va profiter 
de cette propriété pour substituer au spectre du 
courant 1,({) un autre spectre ayant même allure 
autour de &w,, mais dont les parties inopérantes 
seront adaptées de façon adéquate. On s’arrangera 
en effet pour faire correspondre à ce spectre de 
remplacement la transformée d’une fonction du 
temps de forme simple. Nous calquons ce procédé 
sur celui qu’on utilise fréquemment dans la Théorie 
des Réseaux, et qui permet de substituer à une 
impulsion brève une « impulsion réduite » équi- 
valente du point de vue de ses effets [16]. Dans 
notre cas cependant, le courant £,(t) n’a certai- 
nement pas la forme idéale d’une impulsion de 
Dirac et la recherche de son impulsion réduite doit 
tenir compte dans une certaine mesure de l’évolu- 
tion du spectre à l’intérieur de la bande 
d'absorption. 

Si K(w)=H(w) et? est la densité spectrale, à la 
fréquence cyclique w, du champ transversal A(t), 
le fait que ce dernier sera, dans la plupart des cas, 
créé par un courant issu d’un réseau physiquement 
réalisable implique, d’après les relations de 
Bode [17], que la forme du module H(w) n’est pas 
indépendante de celle du facteur de phase ®(o), 
qu’il suffit alors de considérer seul. 

Si nous représentons l’évolution de ® autour 
de «w, par un développement en série 


Do) = Dos) + (wo — wo) (AD/do)s 
+3 (0 — o)° (d° DJde?)s + ... 


nous pouvons admettre que l’épanouissement du 
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spectre en regard de la bande d’absorption permet 
de n’en retenir que les deux premiers termes, ce qui 
constitue la relation d'équivalence de forme : 


pour © æ + wo. Do) # Dos) + (o— wo) (dP do) = 0 
(24) 


Notons qu’un étalement satisfaisant, eu égard à 
notre hypothèse, se traduit par : 


1 D [d «| com 609 & 27 T2. (25) 


Du module, il suffit de ne retenir que sa valeur à 
la fréquence cyclique w, 


pour © æ + &5 H{(o) # H{(w5) (26) 


ce qui permet de conserver l’équivalence énergé- 
tique. 

Dans le cas où le courant est créé d’une autre 
manière, nous admettons qu’il est tout de même 
possible de conserver ces relations d'équivalence en 
considérant entre autres que les variations de 
l’argument, plus rapides que celles du module, ont 
une influence dominante sur la forme de la fonction 
image. 

Nous disposons ainsi, par (24) et (26), de trois 
paramètres pour définir la fonction que nous allons 
substituer à A(t) (les équations correspondant à 
& = — w, ne sont pas indépendantes des précé- 
dentes à cause des propriétés de parité de H et ®), 
et cette latitude permet de centrer sur «, le 
spectre de substitution, ce qui lui fait correspondre 
une perturbation à la fréquence de Larmor dont 
l'amplitude est variable avec le temps. 

Choisissons par exemple : 


het) = Hr sin (ost— ©) pour 0 << O 
Rr(t} =" 0 partout ailleurs. 


Son spectre, par ailleurs bien connu, répondra 
aux conditions (24) et (26) si l’on écrit : 


(27) 


O = — 2{(4®/do)s=%0 (28a) 
? ——[D{w) + r/2] (28b) 
Hr = — 2x mod K{ws).(dP/do)oL we (28c) 


ce qui définit la durée, la phase et l’amplitude de 
l'impulsion réduite de pulsation w,. 


À mod GC («w) 


Sir 


Spectre initial 


F1G. 4. — Découpage du spectre. 
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Fic. 5. — Allure de l’impulsion réduite. 


La largeur du spectre de A,(t) est de l’ordre 

de 4r/O, soit d’après la relation (28a), 
2r|[d®/dofs lo, 

et en tenant compte de l'hypothèse (25), on voit 

que cette valeur est grande devant la largeur 2/T, 

de la raie d’absorption, ce qui justifie l’approxi- 

mation effectuée. 

Traitée dans le référentiel mobile de Rabi [18] 
tournant à la vitesse angulaire — ©, autour de e:, 
il correspond à l’impulsion réduite un champ trans- 
versal ayant la forme d’une impulsion rectangulaire 
de hauteur H,/2 et de durée O. Au temps t — 0 
l'axe 1 de ce référentiel tournant, choisi suivant la 
direction du champ transversal, forme avec l’axe 1 
correspondant du référentiel fixe l’angle — + défini 
par (28c). 

Au temps { — © le champ transversal s’annulle, 
le système de spins est alors abandonné à son évolu- 
tion libre et retrouve son équilibre avec la loi 
exp (— {/T;). C’est donc l’état acquis au temps 
t — © qui nous renseignera sur la linéarité de la 
réponse à la perturbation. 

La présence, pour £{ compris entre 0 et ©, de 
coefficients constants, permet, comme l'a fait 
Torrey [19], de soumettre à la transformation de 
Laplace les équations de Bloch rapportées au repère 
tournant, ce qui fournit le groupe : 


(p + 1/73) U1 = 0 (29a) 
(P + 1/Te) Be + (1 L2) YHr Ba = 0 (29b) 
P + 11Ti) Us — (1/2) YHr Be = Uso + UolPTi, (29c) 


Les u représentent les transformées de Laplace, 
. +co 
a = [07 est pale) de 
—Co 


des composantes ui(t) de l’aimantation suivant les 
axes du repère tournant, 9 est la valeur initiale 


1 1 Lo 
Uo P°? + (7 ts + 7 bo) P FT T 


CES (30) 


1 


Ty 


PRE 
JP+ivx eo 


: 1 
P [r % Fe dE 

D’après la forme du dénominateur, 1l apparait 
que u(t) contiendra des facteurs exponentiels avec 
la constante de temps 2{(1/7;,) + (1/72) 1. 

Or, d’après (28a) et (25), ® est petit devant 7, 
et a fortiori devant T,, ce qui fait que l’évolution 
de u3{t) pendant la durée de l’impulsion réduite 
peut être représentée correctement par un dévelop- 
pement en série autour de { — 0. La propriété bien 
connue : 


limeolti(t)o] = limp>co[(p)] 
nous permet de tirer ce dernier du développement 
asymptotique de (30) : 


ere 1 _n Ho == Hso “a 
EE p Ë 7 p? 
Fee = U30 
(T1)? 


15 1 
+ E V2 HE bso] 23 + ... 


ce qui, par transformation inverse, donne immé- 
diatement pour { — © : 


O 
Us(©) = 30 + (Ho — U3o) T. 
1 


1 


2 [tes Uso) + (31) 


dl O\?2 
Re HET? uso] (7) ter 


L'influence de A, apparaît dans (31) à partir du 
terme quadratique, ce qui permet de prendre 
comme critère de linéarité : 

puae<i+ (mea). (32) 
Dans la grande majorité des cas, m30 diffèrera 


suffisamment peu de m, pour que le second membre 
reste très voisin de 1. Nous en déduisons alors un 


résultat très simple : la condition nécessaire pour. 


conserver leur linéarité aux équations de Bloch en 
régime transitoire est indépendante de l'allure 
exacte de A(t) et se résume, d’après (32), (28a) et 
(28c) à : 


2r|y| mod H{(ows) < 1 (33a) 


soit : 


mod [r Fra e—®ot h(t) a] HA (33b) 
— CO 

La comparaison entre cette condition et le cri- 

tère de non-saturation en régime stationnaire, 

montre que T, et T, n’interviennent plus. Ce résul- 

tat paraît justifié si l’on considère que la durée de 


?. 


N° & 


la perturbation transitoire est nécessairement brève 
devant l’une ou l’autre de ces quantités. 

N. B. — Si, au contraire, la valeur initiale de m3 
était très faible devant m,.0/T,, il faudrait tenir 
compte indirectement de l’allure de A(t) par l’inter- 
médiaire de la vitesse de variation de l’argument de 
son spectre autour de w,. La condition de linéarité 
deviendrait : 

Er 


2 Rlo)]? < 225 
[27y mod R{w0)]? < crrnel 


aD: 
do 


O = Oo 


_ Il semble d’ailleurs que ce dernier cas ne puisse 
correspondre qu’à des perturbations dont l'allure 
se rapproche davantage de celle des perturbations 
quasi-stationnaires (O et donc d® /do importants). 


V. Fluctuations thermiques. — Le bruit de 
fluctuations thermiques dans le circuit électrique 
utilisé pour la mise en évidence des phénomènes 
gyromagnétiques peut être attribué à la résis- 
tance À, de ce circuit, avec une densité spectrale 
énergétique à la température T du circuit donnée 
par la relation de Nyquist : E?{(v) — 4XTR, (*). 

La présence de la résistance À, dans le modèle 
équivalent de la figure 3 suggère de lui attribuer 
un bruit propre dont l’origine réelle ne peut être 
évidemment due qu’à des fluctuations dans la répar- 
tition des populations des niveaux Zeeman. Il est 
alors nécessaire de vérifier si les valeurs calculées au 
moyen du modèle équivalent s'accordent avec la 
réalité physique. 

Sous un aspect macroscopique, considérons le 
moment magnétique résultant M en présence du 
champ permanent H,, en équilibre thermique à la 
température 7” du milieu et soumis en outre au 
champ de rayonnement isotherme incohérent. Il 
en résulte des fluctuations autour de sa position 
d'équilibre, qui, en dehors de la bande de fré- 
_quences de 0 à 1/T,, laissent invariante la norme 
de M : ; 

1 
2 (mx)? = (m0)? = u8 x H3. 

L'énergie potentielle du système par unité de 

volume 


(m:,)? + (m)? | 
- : Aro et 
2(M5) * 


aurait pour valeur à l’état « figé » :W, = — m5 Hs. 
Comme le système macroscopique possède deux 
degrés de liberté il en résulte une énergie moyenne 
de fluctuation dans le volume + de l’échantillon 
qui vaut : 
r.mod (< W, >— Wi) = AT’ 


W, =—(MH;)#— mAH [1 


et en vertu de l’équipartition, on en déduit la 
variance des composantes : 


kT'’ 
nee (rm) — Here [Xol- 


(*) Æ représente ici la constante de Boltzmann. 
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Ce résultat sera conservé lors du mouvement de 
précession causé par la présence supplémentaire 
d’un rayonnement cohérent, tant que la valeur 
moyenne de l’angle entre M et H, reste faible. 
Alors, d’une part, la température des spins demeure 
très voisine de celle du milieu et, d’autre part, les 
équations de Bloch demeurent linéaires et auto- 
risent ainsi l’emploi du modèle électrique équi- 
valent. 

Ce modèle fournit justement, en appliquant la 
formule de Nyquist à la maille 2 considérée comme 
isolée (1, — 0) la relation 


ire KUILS 
soit, d’après le tableau de correspondance (14) : 
€ (nm)? = (u8 K?JarnL).kT/|{0]. 


Or il est facile de voir, d’après les définitions 
de X et de n du $ C, que u$ K?/4rnL #7, ce qui 
conduit pour < (m,)? >. aux mêmes résultats que 
précédemment. 

Ainsi le modèle électrique, qui fournit la même 
variance et le même caractère « microscopique » de 
la fonction de corrélation, s’avère-t-il parfaitement 
apte à traduire quantitativement et jusque dans 
leurs ultimes conséquences tous les phénomènes de 
fluctuations. 


VI. CONCLUSION. 


Nous avons pu voir qu’une équivalence analy- 
tique entre un montage à résonance magnétique 
nucléaire et un réseau électrique de structure parti- 
culièrement simple paraît se justifier pour des con- 
ditions d’emploi qui sont celles d’un appareil des- 
tiné à la seule mesure d’un champ magnétique : 
stabilité et homogénéité du champ, en relation 
même avec la nécessité d’une définition précise de 
Pobjet de la mesure, et travail dans une bande de 
fréquences évidemment voisine de la fréquence de 
Larmor, ou dans le cas transitoire, respect du critère 
de non-adiabaticité. 

Bien qu’il s’agisse en fait d’une fiction mathé- 
matique qui ne peut prétendre se substituer en tous 
points au traitement physique, le fait que l’emploi 
du modèle équivalent puisse s’étendre au régime 
transitoire est d’un grand intérêt, en ce sens que les 
problèmes délicats d’adaptation ou d’interpré- 
tation qu’on y rencontre peuvent trouver une solu- 
tion simple au moyen des méthodes de calcul 
connues de la Théorie des Réseaux. Citons par 
exemple leffet complexe d’amortissement par 
rayonnement cohérent [1] qui, par l'intermédiaire 
du circuit équivalent, peut être ramené à l’étude 
classique de l’évolution libre de deux mailles 
couplées. 

De même, en régime stationnaire, les relations 
d'équilibre d’un pont, les conditions d’oscillation 
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d’un « Maser » en champ faible [6], [10], ou l'effet 
d’entraînement de fréquence [20], etc. peuvent 
être traités quantativement de manière très simple. 


APPENDICE 


Champ tournant déduit d’un spectre de bande. 
— Considérons, dans le repère de la figure 1, un 
champ magnétique transversal rectiligne, dont 
l’expression en intégrale de Fourier s’écrit : 


+oco 


R{t) — H(o) eiot do.e.. 


—00 
Si l’on choisit une bande de fréquences cycliques 
infiniment étroite, centrée sur « et de largeur do, 
la composante élémentaire 


dh — H{w) etot do.e, 
représente un champ oscillant rectiligne cisoïdal 
qu’on écrira : 


dh = H(«) [cos est + rrsin lle] dorer (sa) 
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A ce champ rectiligne on peut faire correspondre 
deux vibrations à polarisation circulaire, et si l’on 
ne retient que celle dont le sens est convenable, en 
relation avec le signe du rapport gyromagnétique 
des noyaux étudiés, on obtiendra un vecteur élé- 
mentaire dh* ayant une vitesse de rotation 


(a) 
RER) EL Ege 


Dans cete correspondance, la contribution du 
terme de (34) en cos ot est (on suppose y positif) : 


= H(w) [cos wt.e, — sin [o|t.e,] do. 


Celle du terme en sin ol! : 
5 H(o) (sin [olt.e, + cos ot.e] y do. 

Soit au total : 

2 


lol 


Manuscrit reçu le 23 novembre 1960, 


dh* — 5 Ho) ciot [e: De e| do. {54 
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DÉVELOPPEMENT DU MOMENT DIPOLAIRE DES MOLÉCULES 
EN FONCTION DES COORDONNÉES DE SYMÉTRIE. I. 


Par GEORGES GRANER, 


Laboratoire d’Infrarouge, 12, rue Guvier, Paris (5°) 
(Laboratoires de Chimie Physique). 


… Résumé. — On rappelle l’état actuel des calculs sur les intensités des bandes d’absorption 
infrarouge, pour les molécules di- et polyatomiques, en mettant numériquement en évidence 
l’importance relative des différents termes du développement du moment dipolaire. 


Abstract. — The present state of computations on intensities of infrared absorption bands is 
reviewed, for di-and polyatomic molecules, showing the relative importance of diverse terms of 


development of the dipole moment. 


I. GÉNÉRALITÉS. 


1-1. Introduction. — Les mesures d’intensité des 
bandes d’absorption dans le domaine infra-rouge 
présentent un intérêt croissant dans différents do- 
maines : astrophysique, cinétique chimique, ana- 
lyse de la structure moléculaire. 

En effet, les intensités font intervenir direc- 
tement les populations des divers niveaux énergé- 
tiques, ce qui permet des évaluations de tempé- 
rature, le moment dipolaire, qui donne des rensei- 
gnements importants sur la liaison chimique, les 
densités optiques enfin, ce qui permet d’estimer 
l'abondance de certaines molécules ou radicaux 
libres dans des flammes ou des atmosphères plané- 
taires. 

La technique des mesures d’intensité a fait de 
sérieux progrès. On sait maintenant, non seulement 
faire des mesures pour toute une bande de vibra- 
tion-rotation, mais, dans des mesures en phase 
gazeuse, pour chacune des raies de rotation de cette 
bande : ceci est maintenant possible grâce aux spec- 
trographes très résolutifs en usage et aux tech- 
niques de correction des largeurs de fente. Nous 
n’insisterons pas sur cet aspect de la question, si 
ce n’est pour dire que la précision des mesures 
d'intensité, qui atteint quelques pour cent pour 
des raies larges dans les cas les meilleurs, reste 
très médiocre dans les autres cas. 

D'autre part, les mesures de dispersion d’indice 
fournissent des informations supplémentaires, voi- 
sines de celles que donnent les intensités [1]. 

Le calcul théorique des intensités permet d’ex- 
ploiter les informations apportées par les mesures. 
Compte tenu de sa difficulté, on ne l’a développé 
jusqu’à ces dernières années que pour les molécules 
diatomiques. Cependant, le problème des inten- 
sités des molécules polyatomiques commence à 


être abordé, en particulier, celui des molécules 
linéaires. 

Nous voudrions exposer rapidement les données 
du problème en mettant en évidence le rôle des 
développements en série du moment dipolaire. 


1-2. Théorie des intensités. — L’intensité d’ab- 
sorption théorique correspondant à une transition 
d’un niveau énergétique +” à un état +’ est donnée 
par : 

Ton OO hr) < «(Mile © 0°) 
3hCuz 
où h est la constante de Planck ; c, la vitesse de la 
lumière ; v, la fréquence en em-t de la transition 
Tr" — +’; N”, la population du niveau initial ; g”, 
son degré de dégénérescence ; X la constante de 
Boltzmann ; T, la température absolue du corps 
absorbant. 

Considérons une bande de vibration corres- 
pondant à un saut de nr unités pour le nombre 
quantique de vibration v. On est amené à calculer 


2e 
< 0 + r[Mlv Pi 


— 
Soit A7, les projections de 47 sur les axes Ox, Oy, 
Oz liés à la molécule (x = x, y, z). On a 


nine a Co Panne 2 
œ 


On est donc ramené au calcul de trois moments 
de transition < v + n|M,]v > c’est-à-dire : 


[rt MA dr 


1-3. Intensité des molécules diatomiques. — 
Pour les molécules diatomiques, il suffit de consi- 
+ 


dérer la composante W, de M sur la droite joignant 
les noyaux. L’intensité est donc proportionnelle à 
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< v + n|M,lv >? pour la transition v — v + nou 
à 5 <v+n|M.\v >? pour l’ensemble de la 


v=0 
bande. | 
C’est à ces calculs d’intensité, toujours relatifs 


aux intensités intégrées des bandes de vibration, 
que se rapportent les travaux les plus anciens [2] 
à [5]. 

Er. des mesures et des calculs ont été faits 
sur les raies de rotation individuelles dans une 
bande de vibration-rotation [1], [6] à [12]. Les 
calculs ont montré qu’il fallait tenir compte de 
l'interaction vibration-rotation, même pour l’inten- 
sité totale des bandes de vibration : cette intensité 
fait en réalité intervenir, non pas < v + n|M|v >? 
mais > <v + n, J'|Mlv, J" >2,J' et J" dési- 

J’,J" 
enant les nombres quantiques de rotation. Dans ce 
qui suit, nous supposerons cette correction 
effectuée. 


1-3-1. CHoOIx DES FONCTIONS POTENTIELLES. — 
Le calcul de < v + n|Mlv = exige la connais- 
sance des fonctions d’onde Ÿ,,, et W,, donc le 
choix d’une fonction potentielle. Dès le début deux 
types de fonction potentielle ont été utilisés : 

— La fonction de Morse [1], [2], [3], [7]. [9] 
[12], [13], [14], [15] 

Ul(r) = def1 re))l?. (3) 

— La fonction de l’oscillateur harmonique per- 
turbé, contenant des termes cubiques, quar- 
tiques.., qui rend bien compte des phénomènes près 


de la position d’équilibre de la molécule [1], [3], [5], 
[6], [7], [9], 107, [13], [14], [15] 


U(r) = aolr 


exp (— a(r 


re)? + 


où r représente la distance internucléaire et re sa 
valeur à l’équilibre. 

Un certain nombre d’autres fonctions poten- 
tielles, ont été proposées, qu’on trouvera récapi- 
tulées rapidement dans [1] et d’une manière 
exhaustive dans [16]. 

Il convient de rappeler que c’est en général pour 
obtenir une bonne vérification des fréquences ex pé- 
rimentales qu’on choisit l’une ou l’autre de ces 
fonctions potentielles. Les mesures d'intensité qui 
sont considérablement moins précises que les 
mesures de fréquence n’interviennent pratiquement 
pas dans ce choix. 

L'utilisation de la fonction de Morse permet un 
calcul rigoureux d'expressions du type 


JF @.q".#"(g).dg. 


Pour l’oscllateur anharmonique, le calcul de telles 
expressions se fait par les méthodes de pertur- 
bation. Ces méthodes ont été rendues plus systé- 
matiques par la méthode de transformation de 
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N° 4 
contact, exposée dans [10] et exploitée dans [9] 
et[11]. 


1-3-2. CHoix DE LA FONCTION M(r). — La forme 
de la fonction M(r) donnant le moment dipolaire en 
fonction de la distance internucléaire n’est connue 
qu’à partir des mesures d’intensité et de quelques 
considérations de chimie théorique [4], [17]. Ceci 
explique que la quasi-totalité des auteurs aient 
utilisé des développements limités de M(r) au voisi- 
nage de re : 


M = M, + Mifr— re) + Milr — re)? + Malr — re)? 
(5) 

ou encore 
M = do + d19 + di19* + din 9° + (6) 


en désignant par g la coordonnée normale sans di- 
mension correspondante 


se — re) 4 


we et Be désignant les constantes spectroscopiques 
de la molécule. 

Un tel développement limité ne se justifie que si 
les termes successifs vont en décroissant, c’est-à- 
dire en fait, pour de petites valeurs de (r — re)fre. 

Nous montrerons plus loin par quelques exemples 
numériques que l’on ne peut pas plus négliger les 
termes M,, M; ... qui définissent l’anharmonicité 
électrique qu’on ne peut négliger l’anharmonicité 
mécanique. 

Récemment des auteurs [12], [15], [18] ont uti- 
lisé une nouvelle méthode de développement de 
M(r) : on suppose connus, à partir des mesures 
d'intensité, n éléments de matrice A", c’est-à-dire 


F0 = € v|M(r)l0 >. (8) 


Dans ce cas, M(r) sera approché par l'expression 


Mr) = SRE). (9) 


v=0 


Ce développement, rigoureux pour n infini, pos- 
tule lorsqu'on l’arrête au n°" terme, que les 
coefficients R°0 suivants tendent rapidement vers 
zéro, ce qui semble justifié par l’expérience. Le 
choix d’une fonction potentielle correcte devient 
crucial par cette méthode. 


1-3-3. EXEMPLES DE DÉVELOPPEMENT DE MO- 
MENT DIPOLAIRE. 

Cas DE HCI. 

Nous emprunterons les formules qui suivent et 
les valeurs numériques à Legay [1]. D’après ses 
mesures, les coefficients 47; de la formule (5) sont : 


M; — 0,993 DJÀ 
Ms = —1,093 DJ(AŸ: 


M, = 1,138 D 
M3 = — 0,309 DJ(À)? 


N° 4 


Les éléments de matrice R°0 du moment dipolaire 
(ef. formule (9)) s'expriment en fonction des HW; par 


ñn 
(Ro)? — [È ea HP | (10) 

i=1 
ou x représente une constante de la molécule 
x — AIT? cu«./h et où R?° est l’élément de matrice 
< v|g\ 0 >, q étant la coordonnée sans dimension 
donnée par (8). 

Ces éléments R?° sont calculés dans [1] en utili- 
sant les fonctions d'onde de l’oscillateur anharmo- 
nique. 

Pour la bande fondamentale, (10) s’écrit : 


QRIO) = 7 (M, + AN M, + à F1 M.)° 
= —1(0,70750 — 0,010317 — 0,026667)2. (11) 


Dans cette expression, c’est le terme , qui joue 
le rôle principal. Mais si, comme il est classique de 
le faire, on se contente de ce terme, on commet une 
erreur de 10 % sur l'intensité, ce qui est loin d’être 
négligeable. 

Pour la bande première harmonique, (10) devient : 
(Aro) = o (RU M, + oo VI M, to RM) 
— g—1(— 0,04265 — 0,02188 + 6)2. (12) 


Le terme en M, et le terme en #7, jouent donc un 
rôle comparable et ni l’un ni l’autre ne peuvent 
être négligés. 

Pour la bande deuxième harmonique, on a : 
(R30}2 — En ge M, + RTE LES M, + Fra us M:)° 
= x—1(0,011817 + 0,0038775 — 0,010888)2 (13) 


Donc les trois termes jouent des rôles compa- 
rables. 

Les considérations précédentes montrent qu’au 
moins dans le cas de HCI, il est nécessaire d’utiliser 
au moins 4 termes dans le développement de la for- 
mule (5) si on désire calculer les intensités de la 
deuxième harmonique, de la première, et même de 


la bande fondamentale. 

Ces conclusions restent valables si on utilise la 
fonction potentielle de Morse au lieu de celle de 
l’oscillateur anharmonique, si on applique les for- 
mules de Dunham [3] ou de Crawford et Dins- 
more [5], ou encore les valeurs expérimentales 
d’autres auteurs. 


Cas DE CO. 


Nous emprunterons les formules qui suivent et 
les valeurs numériques à Herman et Shuler [14]. 
D’après leurs mesures, les valeurs les plus pro- 


bables des coefficients de la formule (6), sont : 
d) = 1,470 X 10—1. d,, = — 1,307 x 105 
dire 2,040 107. 


Les éléments de matrice R°” du moment dipo- 
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laire s'expriment en fonction de ces coefficients par 


(REG RTS ET EE Le 0) Rte 


où A?% est l’élément de matrice < v’Iglv" >. 


Nous avons recalculé un certain nombre de 
ces A?% et nous obtenons les résultats suivants : 
Pour la bande fondamentale, (14) devient : 


(R10)2 = (0,10430 — 0,00025 — 0,00075)2. 


Le terme en d, est de loin le plus important. 
Négliger d,,etd,,, revient à commettre sur (R10)2 
une erreur de 0,2 %, ce qui est négligeable. 

Pour la bande première harmonique, (14) donne : 


(R20)2 = (— 5,203 X 10—8 — 0,91 x 10—3 + 0}. 


Le terme en d, est plus important, mais négliger . 
le terme d,, reviendrait à commettre une erreur 
de 25 %, sur l'intensité. 

Pour la bande deuxième harmonique 
devient : 


(14) 


(R30)2 = (5,664 X 10—4 + 1,167 x 10—4— 3,066 x 10—4)2. 


Les trois termes sont du même ordre de grandeur. 

Les conclusions sont presque les mêmes que 
pour HCI : le calcul de l’intensité de la bande pre- 
mière harmonique exige l’utilisation de d, et de d,,, 
celui de la bande deuxième harmonique l’utilisation 
de d;, dis, dinse 

Nous n’avons pas trouvé de travaux concernant 
d’autres molécules diatomiques que HCI et CO, et 
qui permettent les mêmes discussions sur l’impor- 


tance relative des termes du moment dipolaire. 


1-4.-Cas des molécules polyatomiques. — La 
résolution du problème des intensités des molécules 
polyatomiques est bien plus difficile et bien moins 
avancée que celle des molécules diatomiques. Les 
mesures d'intensité ne sont guère plus difficiles et, 
au moins pour des molécules simples, un bon 
nombre de résultats ont été publiés : tout au plus 
peut-on signaler que la plupart des mesures con- 
cernent les intensités intégrées sur toute une bande 
de vibration et non, comme cela se fait de plus en 
plus pour les molécules diatomiques, les intensités 
des raies individuelles de rotation. 

La partie la plus difficile est l’interprétation de 
ces mesures en termes de moment dipolaire et de 
fonction potentielle. En effet, pour une molécule à 
n atomes, il y a 3n — 6 (ou 3n — 5 pour une molé- 
cule linéaire) coordonnées normales g;. Ceci veut 
dire que le développement du moment dipolaire, 
qui a cette fois 3 composantes, comprend, a priori, 
3 termes constants, 3(3n — 6) termes du premier 
degré par rapport aux g, 3(3n — 6)(3n — 5)/2 
termes du deuxième degré, au lieu de 3 termes en 
tout pour les molécules diatomiques. 

De même si la fonction potentielle est celle d’un 
oscillateur anharmonique, il y a (3n — 6) termes 
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du deuxième degré, (3n — 6) (3n — 5) (3n — 4)/6 
termes du troisième degré... 

Cette complexité explique qu’on puisse classer 
les calculs théoriques sur les intensités de molécules 
polyatomiques en deux types : 1° ceux où on se 
préoccupe d’interpréter les résultats de manière 
simple et quasi-intuitive, de manière à ce que le 
nombre de paramètres inconnus à déterminer soit 
inférieur ou égal au nombre de mesures effectuées ; 
29 ceux où la rigueur théorique est grande mais 
qui, pour l'instant, sont invérifiables faute de 
données expérimentales. 


4-4-4. PREMIER TYPE DE CALCULS. — Îls uti- 
Bsent un certain nombre d’hypothèses, de validité 
souvent discutable : 


Hypothèse 1 : L’intensité d’une bande fondamen- 
tale v; est proportionnelle au carré de la dérivée 
première de 47 par rapport à la coordonnée nor- 
male g; corres pondante 


—\ 2 
dM 
(2%) 


Cette hypothèse sous-entend l’harmonicité élec- 
trique et mécanique. Mais, nous avons vu dans 
l'exemple de HCI et de CO, qu’elle approche d’assez 
près la réalité, même si ces deux harmonicités ne 
sont pas réalisées. Une telle hypothèse ne peut être 
faite pour des bandes harmoniques. 

— 


(15) 


est déterminé par les 


Dans ces conditions, | — 
i 
mesures d'intensité. Si une composante seulement 


— 
de M, M; par exemple intervient, on a alors 
dM,/0dq; au signe près : 


d M 
da 


dMx 
— "+ È 16 


ee 
dM 
Page les 


Hypothèse 2 : La fonction potentielle est harmo- 
nique. 

Compte tenu de cette hypothèse et des fré- 
quences des bandes observées dans le spectre infra- 
rouge de la molécule, on peut exprimer numéri- 
quement la fonction potentielle en fonction des 
coordonnées de déplacement : allongements Ar des 
liaisons, variations A des angles de valence, etc... 
Il est alors possible de trouver les relations qui 
relient les coordonnées normales g; aux coordonnées 
symétriques 153 (1) 

Sy — > Luigi où qi = > (L— 4e 154. (17) 

Les coordonnées symétriques peuvent elles- 
mêmes être exprimées uniquement à partir de con- 
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sidérations géométriques (1) en fonction des coor- 
données internes Ëy 


Si —=Dua 4 ou 4 = 2 (ur—1)y S%. (18) 


Les relations (3) permettent de déterminer les 
quantités 2, [084 car : 


dMe f dMe d dMy 
= — = Leaf, . 
DS À ou 38e — à LU So 


Il peut y avoir une incertitude sur les signes 


(19) 


* Steph 


des dW,[dq ce qui entraîne plusieurs valeurs pos- 


sibles pour 2,108: 
Ensuite les relations (18) donnent les dérivées 
par rapport aux coordonnées internes 


Me LS y Me 
CEE dS4 


Hypothèse 3 : Additivité des moments de liaison. 
Cette hypothèse a été introduite par Mecke [19], 


(20) 


— 
[20]. Le moment dipolaire M de la molécule est la 


somme de moments dipolaires partiels attribués à 


chaque liaison et ayant la direction de la liaison 
correspondante. 

Soit LL Le ... y Ces moments de liaison. 

Il existe des relations faciles à trouver entre les 


EL 
composantes 7, de M et les moments de liaison uy. 
Soit : 

M4 = 2 Dex. Uk. (21) 
Les ®,4 sont les cosinus des angles des liaisons 
avec les axes Ox, Oy, Oz. 
D’où 


ù ex) 


dM, d® 
= pa ( CRE LiDE: 
k dE 


dE dE 

Les dérivées d0,:/dE sont calculables à partir de 

considérations géométriques. Les dérivées 2M,/dE 

sont données par les équations (19) et (20). Il en 

résulte qu’on peut, à partir de (22) calculer les 

moments de liaison u, et les dérivées du,/0Ë, de ces 

moments de liaison (?) par rapport aux coor- 
données internes de la molécule. 

Cette méthode est très couramment utilisée [21] 

à [32] mais sa validité a été mise en doute [33], [34] 


parce qu’elle amène à des résultats inexplicables. » 


D’une part, les diverses vibrations de type de 


symétrie différents, donnent des valeurs très diffé- 


rentes pour la charge effective d’une même liaison 
[23], [26]. D'autre part, les moments de liaison puy 
déduits des mesures de bandes de déformation ont 
des valeurs différant nettement de celles que don- 
nent d’autres méthodes. 


(1) Cf. Wilson [43], p. 115. 

(?) Lorsque la coordonnée interne & 
gement Arz de la liaison dont le moment est ux, la quan- 
se dux [06 est souvent appelée « charge effective » de la 
iaison, 


(22). 


est l’allon- 


N°4 


C’est pourquoi Sverdlov [34] à [39] a jugé utile 
de modifier quelque peu les approximations utili- 


sées sans pour cela rejeter les trois hypothèses 
précitées. 


1-4-2. DEUXIÈME TYPE DE CALGULS. — Un deu- 
xième groupe d'auteurs rejette les trois hypothèses 
prsdentes et considère le moment dipolaire glo- 


bal M comme une fonction des coordonnées nor- 
males dont les seules propriétés connues sont celles 
que donne la théorie des groupes. Par commodité, 
on considère ses composantes M, comme dévelop- 
pables en polynômes par rapport aux coordonnées 
normales. 

Ensuite le calcul des éléments de matrice du 
type < v'[M,|v" > se ramène à des calculs de 
terme du type < v’[(g}"|v” >, c’est-à-dire à des 
problèmes de détermination de fonctions d’onde. 

Ce genre de calculs a été fait par la méthode 
« classique » des perturbations par Eggers et 
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Crawford pour CO, et N,0 [33]. Par la suite, la 
méthode de la transformation de contact a été 
appliquée aux calculs d'intensité des molécules di- 
et polyatomiques [10] puis au cas général des 
molécules polyatomiques [40], des molécules liné- 
aires [40], [41], des toupies symétriques [42]. 

Ces méthodes permettent d’obtenir des relations 
entre les constantes d’anharmonicité du potentiel 
et les coefficients du moment dipolaire, mais les 
applications directement pratiques sont quasi- 
inexistantes. Un des écueils est que, contrairement 
à ce qui se passe pour les molécules diatomiques, on 
manque des données suffisantes pour connaître 
avec beaucoup de précision les fonctions poten- 
tielles, du moins dans leur partie anharmonique. 

Le but du présent article est d’éclaircir un peu le 
problème en montrant le rôle que jouent les coef- 
ficients du développement du moment dipolare 
dans les calculs d’intensité et, en utilisant la théorie 
des groupes, de trouver des relations importantes 
entre ces coefficients. 


BIBLIOGRAPHIE 


[1] Lecay (F.) 
[2] Dunxam (J. 
[4] Dunxam (J. 
[4] MULLIKEN | 


, Thèse, Paris, 1957. 

), Phys. Rev. , 1929,84, 438. 

), Phys. Reo., 1930, 35, 1347. 

84, J. Chem. Physics, 1934, 2,400 et 712. 

[5] CrawrorDp (B. L.) et Divsmore (H. Se J. Chem. 
Physics, 1950, 18, 983 et 1682. 

[6] HERMAN (R.) et WALLIS (RL), 


iB£ 
L. 
IR 


J. Chem. Physics, 


1955, 23, 637. 
[7] Herman (R.) et Rugin (R. J.), Astrophys. J., 1955, 
121, 533. 


8] Lecay (F.), Cahiers Physique, 1956, 75, 11. 
9] Herman (R.) et Rugin (R. J.), Canad. J. Physics, 
1956, 34, 1119. 

[10] Hawson (H.), NiezseN (H. H.), Suarrer (W. H.) et 
WaAGconer (J.), J. Chem. Physics, 1957, 27, 40. 

[11] Lecay (F.), Cahiers Physique, 1958, 99, 416. | 

[12] Herman (R.) et Rugin (R. d.), J. Chem. Physics, 
1960, 32, 1393. 

[13] Herman (R.) et Saurer (K. E.), J. Chem Physics, 


1953, 21, 373. 
[14] Herman (R.) et Saucer (K. E.), J. Chem. Physics, 
1954, 22, 481. 
[15] Triscuka (J.) et SALWEN (H.), J. Chem. Physics, 1959, 
31, 218. 
[16] Varsani (Y. P.), Rev. Mod. Physics, 1957, 29, 664. 
[17] Gray (B. F.) et PrircnarD (H. O.), J. Mol. Spect., 


1958, 2, 137. 

[18] Bzinper (S. M.), J. Chem. Physics, 1960,82, 582. 

(ol Tim (B.) et MEcxe (R.), Z. Physik, 1935, 98, 373. 

[20] Mecxe (R.), Z. Physik, 1937, 107, 595. 

[21] ScxaATz (P. N.) et Hornic (D. F.), J. Chem. Physics, 
1953, 21, 1516. 

[22] WILLIAMS (D. G.), Person (W. B.) et CrAwronD (B.), 
J. Chem. Physics, 1955, 28, 179. 


[23] MacKEaAN (D. C.) et ScnarTz (P. N.), J. Chem. Physics, 
1956, 24, 317. 
24] MAcKEaAN (D. C.), J. Chem. Physics, 1956, 24, 1002. 
25] GoziKke (R. C.), Mizzs (I. M.), Person (W. B.) et 
CRAWFORD (B.), J. Chem. Physics, 1957, 26, 1266. 
] Mayuoop (J. E.), Canad. J. Phys., 1957,85, 954. 
] Nyquisr (I. N.), Mizzs (I. M.) et Person (W. B.), 
J. Chem. Physics, 1957, 26, 552. 
] Dicxson (A. D.), Mizzs (I. M.) et Crawrorp (B.), 
J. Chem. Physics, 1957, 27, 445. 
] Seuuzrz (J. W.) et Eccers (D. F.), J. Mol. Spect., 
1958, 2, 113. 
[30] Mizzs (I. M.), Person (W. B.), Scuerer (J. R.) et 
CrAwWFORD (B.), J. Chem. Physics, 1958, 28, 851. 
[31] ScHarTz (P. N.) et Levin (I. W.), J. Chem. Physics, 
1958, 29, 475. 
[32] Barton (D. R.)et THompson (H. W.), Proc. Roy. Soc., À, 
1959, 251, 475 et 458. 
[33] Eccers (D. F.) et CrAwrorD (B.), J. Chem. Physics, 


LOAMOMISS 
[34] SVERDLOV (L. M.) Opt. 1. Spekt., 1959, 7, 600. 
[85] SverpLov (L. M.), Opt. 1. Spekt., 1958, 4, 697. 
[86] SverpLov (L. M.), Opt. 1. Spekt., 1959, 6, 729. 
[37] SverpLov (L. M.), Opt. i. Spekt., 1959, 7, 21. 
[38] SverpLov (L. M.), Opt. 1. Spekt., 4959, 7, 452. 
[39] SverpLov (L. M.), Opt. : Spekt., 1960, 8, 96. 
[40] LeGay (F.), Cahiers Physique, 1958, 99, 416. 
[41] WAGGonER (J.) et SHAFFER (W. H.), J. Mol. Spect., 


1960, 4, 224. 
] Hanson (H. M.) et Nrezsen (H. H.), J. Mol. Spect., 
1950, 4, 468. 
[43] Wizsow (E. B.), Decrus (J. C.) et Cross (P. C.), Mole- 
cular Vibrations, McGraw-Hill Book Company, 
New York, 1955, 


[42 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


FLUORESCENCE, PHOSPHORESCENCE ET « RÉDUCTION » 
D'UN ÉCHANTILLON DE WILLEMITE ACTIVÉ AU MANGANÈSE 


Par Hector À. ARDUINO, Juan C. GUERCI, Anprea LEVIALODI et Norgerro MAJLIS, 


Departemento de Fisica, Facultat de Ciencias Exactas y Naturales, Buenos-Aires, Argentine. 


Résumé. — L'utilisation d’un phosphoroscope à oscilloscope nous a permis de photographier 
les courbes d’excitation et d'extinction d’un échantillon de willemite Zn,SiO, (Mn 1 %) entre 15 °C 
et 200 oC. 

Sur les clichés agrandis, nous avons mesuré la somme de lumière et nous l’avons comparée à 
la lumière totale émise à l’extinction. L’analyse des courbes nous a permis de définir deux vies 
moyennes fluorescentes et une phosphorescente. Nous avons déterminé la « réduction » (variation 
relative du rendement) qui traduit une variation dans le transfert d’énergie interne avec la tempé- 
rature. Cette note n’est qu’une introduction à la mesure du rendement absolu au moyen d’un 
radiomètre différentiel absolu qui va être mis au point chez nous. 


Abstract. — The use of an oscillograph phosphoroscope gave us the possibility of photographing 
the excitation and decay curves between 15 °C and 200 °C of a willemite sample Zn,SiO, (Mn1 %). 

By means of an enlarged negative, we measured the sum of light and compared it with the 
total extinction light. The analysis of these curves allowed us to define two different fluorescent 
life-times and one phosphorescent one. We measured the so-called ‘‘ reduction ” (relative varia- 
tion of efficiency) which corresponds to a variation with temperature of the internal energy transfer. 
This paper is an introduction to a research on the measurement of the absolute efficiency by means 
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of an absolute differential radiometer. 


Introduction. — Nous avons observé une diffé- 
rence entre la loi suivie par la courbe d’émission 
sous excitation et la courbe d’extinction. Cette 
différence dépend beaucoup de la température dans 
le domaine observé (entre 15 et 200 0C). En outre, 
le rendement lumineux, comme on le sait, dépend 
de la température. 

Un travail qui a déjà fait l’objet d’une publi- 
cation antérieure [1] a été poursuivi dans le but 
de mesurer la variation de la « réduction » (quen- 
ching) en fonction de 7. Comme on le sait, toute 
émission lumineuse où les centres excités ne sont 
pas soumis à interaction avec d’autres atomes 
pendant le temps qui sépare l’absorption de l’émis- 
sion, donne une courbe d’émission, sous excitation, 
égale à la courbe complémentaire de la courbe 
d'extinction lorsque les centres excités sont un 
faible pourcentage de tous les centres présents (1). 
Ceci est en particulier certain pour des courbes 
exponentielles. 

Dans ces conditions, le rendement sous excitation 
et le rendement à l’extinction sont complémen- 
taires ; le phosphore restituant pendant l’extinction 
toute l’énergie lumineuse emmagasinée pendant 
l'excitation. 

Nous disons que la surface A4’ est complémen- 
taire de :B quand À — 4’ (fig. 1). 

Nous appelons rendement pour une longueur 
d’onde d’émission em excitée par une longueur 
d’excitation Xex, le rapport entre le nombre de 


(1) Loin de la saturation. 


L) 


photons émis et le nombre de photons d’excitation 
absorbés. 

Si, de plus, la courbe d’émission sous excitation 
est exponentielle avec une vie moyenne Tem et Si 
la vie moyenne à l’extinction est +, 1l arrive que : 


Tem — Tex. 


Nous admettons, en plus, en accord avec presque 
tous les auteurs, que même lorsqu'il y a des inter- 
actions, qui enlèvent toute symétrie entre courbe 
d’excitation et courbe d'émission, le coefficient 


PTE 1 


d'absorption du phosphore à la radiation excita- 
trice n’est pas modifié par la température. Dans 
ces conditions, toute dissymétrie doit être attribuée 
exclusivement à linteraction qui modifie le trans- 
fert interne d'énergie. Il en résulte une variation du 
rendement. C’est cette variation relative de rende- 
ment que nous appellerons ci-dessous «réduction ». 

Nous pouvons définir une « réduction instan- 
tanée », une « réduction totale », une « courbe de 


— "7 ee 
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réduction » en fonction du temps et une «courbe de 
réduction » en fonction de la température. Les 
définitions en seront les suivantes : 

Réduction instantanée : Rapport entre ordonnées 
à l’instant considéré ; soit rapport entre ordonnée 
de la courbe d’extinction et ordonnée de la courbe 
complémentaire de la courbe d’excitation choisies 
pour des temps égaux, en prenant pour origine le 
début des deux processus différents. 

Réduction totale: Rapport entre la surface totale 
limitée par la courbe d’extinction et la surface 
limitée par la courbe complémentaice à celle d’émis- 
sion. 

Courbe de réduction en fonction du temps : Définie 
par les ordonnées qui donnent les rapports entre 
surfaces partielles successives limitées par les 
courbes d’extinction et la complémentaire de la 
courbe d'émission entre {, — 0 et £ variable. 

Courbe de réduction en fonction de la température : 
Réduction totale en fonction de la température. 
Si l’excitation est assez intense pour amener le 
phosphore très près de la saturation, une variation 
moindre de lexcitation peut produire une grande 
variation de la courbe d’excitation sans presque 
changer la courbe d’extinction. 

Dans nos expériences, nous sommes certainement 
très loin de la saturation. En effet, bien que nous 
ne connaissions pas encore le rendement absolu, 


nous pouvons admettre que la radiation excitatrice 


(A = 2537 À) a une intensité constante et son 
absorption par le phosphore ne varie pas entre 15 
et 200 oC. 


PHOSPHOROSCOPE UTILISÉ. — Dans les figures 2 
et 2 bis on montre le dispositif utilisé, légèrement 
modifié par rapport. au schéma donné précé- 
demment [1]. 

A : Lamelle qui permet de supprimer les radia- 
tions ultraviolettes en laissant seulement le visible. 

B : Lamelle déplaçable qui permet de couper 
toute la radiation (y compris l’ultraviolet). On 
VPutilise pour obtenir une onde d’excitation carrée. 

C : Porte échantillon de métal pour y loger le 
phosphore. 

D : Système de chauffage électrique à résistance 
(en courant continu). 

E : Thermopile. 

F : Photomultiplicateur IP 21 R. C. A. 

G: Tube germicide (germicide tube) servant de 
source d’ultraviolets 8 watts Sylvania (en courant 
continu). 

H : Transformateur pour allumer le tube par 
impulsions de courant alternatif. 

I : Oscilloscope Tektronix No. 545. 

J : Relais Siemens 3402/5 modifié. 

K-L : Commutateur mécanique actionné par un 
moteur synchrone. 

Le fonctionnement du dispositif expérimental est 
le suivant : Le faisceau ultraviolet est arrêté pério- 
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Ca 


Il 
Tektronix 
545 


IG 2 bis. 


diquement par la lamelle transparente au visi- 
ble (B). La lamelle B est actionnée par le relais (J) 
qui est commandé par le commutateur K-L qui, à 
son tour est actionné par le moteur synchrone. Le 
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faisceau ultraviolet dont la largeur est égale à 
celle de la cuvette porte-échantillon, soit 3 mm, 
est intercepté en un temps { à vitesse constante. 

Dans une série d'observations le temps d’exci- 
tation utilisé était : , — 1,125 seconde et le 
temps d’extinction : 4 — 3,375 secondes, donc #; 
et t, très grands par rapport au temps de coupure 
(é — 0,002 seconde). La lamelle B est presque par- 
faitement transparente au visible de façon que nous 
puissions obtenir une impulsion ultraviolette carrée 
superposée à une irradiation continue. 

Le faisceau, ainsi modulé, tombe sur l’échantillon 
de phosphore dont l’émission excite le phototube 
EM PM ROC. À 

L’impulsion électrique créée dans le photomul- 
tiplicateur est appliquée au préamplificateur du 
Tektronix et simultanément déclenche le «trigger » 
au moyen d’un commutateur (J-L). Le déclen- 
cnement s’effectüe 0,01 seconde avant que le fais- 
ceau ne soit coupé. Le chauffage de l’échantillon 
se fait par une résistance alimentée en courant 
continu pour éliminer tout bruit de fond ; sa tempé- 
rature est mesurée par une thermopile. La lamelle 
semi-transparente qui intercepte l’ultraviolet (A) ; 
est disposée de telle sorte que l’on obtienne une 
impulsion carrée visible servant de repère quand(B) 
fonctionne comme précédemment. 

Afin de permettre l’observation du phénomène 
pendant un temps plus prolongé, nous avons appli- 
qué un balayage exponentiel à l’axe des X ; de 
cette façon, les modifications en ordonnées devien- 
nent plus évidentes. Les temps ont été mesurés en 
superposant des impulsions électriques séparées par 
des intervalles de temps fixes. 


Résultats. — Les résultats de nos expériences 
sont représentés sur les divers diagrammes ci-Joints. 


Fc. 3a. — 20 0C. 


Nous avons d’abord déterminé la variation de 
Pémission sous excitation et de l'extinction en 
fonction du temps (paramètre : la température). 
Les photographies de la figure 3 en donnent 
2 exemples, respectivement à 20 et à 100 degrés. 
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Nous en avons déduit les courbes normalisées en 
échelle semi-logarithmique représentées par les dia- 


FIG: 35.100" °C: 


grammes des figures [4:, 4», 43, 44, et 45], qui 
représentent à diverses températures comprises 
entre 25 et 200 ©C, la variation de l’intensité pour 
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FIG. 4e. 


cent normalisé de la décroissance (extinction sans 
excitation) : courbe D, et de l'excitation : courbe E. 

L’intensité de l’émission sous excitation et celle 
de extinction en fonction de la température (para- 
mètre : le temps) sont représentées par le dia- 


N° 4 


10° 1% normalisé 
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Fic. 5. — Émission sous excitation 
en fonction de la température (paramètre : le temps). 
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gramme de la figure 5 pour t — 0 (intensité de 
SALULALION) MAN S7 TO PE NES 10! 
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FiG. 6. — Réduction en fonction de la température. 


La figure 6 donne la «réduction » en fonction de 
la température en prenant comme paramètre le 
temps. 

La figure 7 donne la «réduction » en fonction du 


F1G. 7. — Réduction en fonction du temps. 


tempsen prenant comme paramètre la température. 
L’intensité d'émission sous excitation (courbe 

supérieure) et d’extinction (courbe inférieure) en 

fonction du temps est rapportée à l’équation : 


æ = Al1 —e—vt). 


Le problème est un cas compliqué de transfert 
d’énergie interne. Nous n’avons pu étudier jusqu'ici 
que la variation relative du rendement. Nous 
mettons actuellement au point un radiomètre 
absolu pour parfaire ces études. 

Néanmoins nous pouvons dès maintenant tirer 
quelques conclusions. 

L’intensité de saturation relative varie linéai- 
rement avec la température (fig. 4) ce qui prouve 
que les centres d’émission sont réduits linéairement. 
En effet, on connaît bien la relation linéaire entre 
émission et excitation à température constante et 
loin de la saturation. 

La réduction en fonction du temps suit entre 60 
et 200 0C la loi suivante : 

19 Pendant un temps inférieur à 10 millisecondes 
la variation (de réduction) est linéaire et rapide et 
ensuite presque constante (variation linéaire très 
lente). 

29 La réduction en fonction de la température 
suit une loi très semblable. Si la température avait 
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été modifiée linéairement avec le temps, avec un 
gradient particulier, la courbe de réduction ther- 
mique et celle de réduction avec le temps pourrait 
se comparer. Cette comparaison nous parait être 
assez proche de celle entre « thermal glow » et 
«after-glow ». 

30 Toutes les courbes normalisées d'extinction et 
d’excitation montrent que leur décalage dépend 
beaucoup de la température. 

Il y a deux cycles en évidence : un cycle rapide 
et un cycle lent. La réduction avec le temps varie 
beaucoup plus pendant le eyele rapide que pendant 
le cycle lent. 
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Nous avons ensuite appliqué la méthode de dé- 
composition en exponentielles [2] à nos courbes, en 
considérant toujours des temps totaux tels que 
l'émission finale observée soit inférieure à 1 % à la 
valeur initiale. Nous avons ainsi obtenu : 

1) Deux émissions fluorescentes de vie moyenne : 
ty 2, O2 108$ et mr = 5508 AUS 
indépendantes de la température. 

2\ Une émission phosphorescente due à une 
famille de pièges autour de la valeur : 


E — 0,41 + 0,02 eV (S = 108 s—1). 


Manuscrit reçu le 15 octobre 1960. 
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CrAwrorD (A. E.), Technique des ultra-sons, traduit par 
J. Palmé. (1 vol. 13,5 X 21,5 cm, 456 pages, Dunod, 
Paris, 1959, prix : 48 NF.) 

Les ultra-sons ont rapidement quitté le domaine du labo- 
ratoire pour faire l’objet de nombreuses applications indus- 
trielles. Aussi l’ouvrage de B. Crawford groupant toutes ces 
applications est-il particulièrement intéressant. 

Cet ouvrage comprend trois grandes parties : 

Un rappel descriptif sur la propagation des ondes ultra- 
sonores et le phénomène de cavitation, ces deux points 
étant à la base de la majorité des applications. Ce rappel 
a été précisé de manière plus théorique par J. Palmé qui 
cite dans ce même chapitre les formules principales régis- 
sant les ondes ultra-sonores. | 

— Un chapitre sur les émetteurs : phénomène de piézo- 
électricité, magnétostriction, générateurs à veines fluides. 

— Un chapitre d'applications qui constitue la partie prin- 
cipale de l’ouvrage (applications métallurgiques, chimiques, 
biologiques). Cette partie de l’ouvrage est autrement riche 
et parfaitement documentée. On regrettera cependant qu’un 


effort de synthèse sur le rôle des ultra-sons pour un phéno- 
mène physique donné n’ait pas été suffisamment fait. Les 
paragraphes supplémentaires introduits par J. Palmé sont 
précieux à cet égard car ils résument et concluent. 

Ch. MorerTrTi. 


Sync (J. L.), The relativistie gas. (1 vol. 16 X 22,5 em, 
120 p., 15 fig., North-Holland publishing company, 
Amsterdam, 1957, prix : 15 f1.). 

Cet ouvrage fait suite à « Relativity ; the special Theory» 
du même auteur, et en utilise les notations (en particulier 
le formalisme à quatre dimensions de Minkowsky) et les 
résultats. 

On y étudie un gaz relativiste composé de particules (ou 
éventuellement de photons) sans dimension interagissant 
uniquement par collision directe, les phénomènes relevant 
de la théorie générale de la relativité ou de la mécanique 
quantique n’étant pas pris en considération. On suppose 
que les photons sont des particules d’énergie äv, de mo- 
ment äv/e dont le quadri-vecteur d’énergie-impulsion se 
trouve sur le cône de lumière et obéissant aux lois de la 
statistique de Maxwell-Boltzmann. 

Après quelques considérations géométriques sur l’espace 
des quadri-vecteurs d’énergie-impulsion d’une particule de 
masse »2 (ou d’un photon), l’auteur définit une fonction de 
distribution invariante et un tenseur d’énergie-impulsion du 
gaz. 

Puis par des considérations statistiques, il montre qu’en 
l’absence de forces extérieures la distribution la plus pro- 
bable est celle de Maxwell-Jüttner. Il en déduit les équa- 
tions du mouvement du gaz considéré, la distribution 
d’énergie, définit sa température, et la vitesse la plus pro- 
bable d’une particule du gaz. 

Après en avoir tiré quelques conséquences simples, l’au- 
teur consacre son dernier chapitre à l’étude des ondes de 
choc dans les gaz relativistes et, en appendice, discute cer- 
taines formules sur les champs de rayonnement dans les 
gaz dues à Li. H. Thomas (1930). 

Bien que n’abordant qu’un petit nombre de questions, cet 
ouvrage, d’une clarté remarquable, sera utile à celui qui 
désire étudier le comportement du gaz relativiste. 

À. MANGENEY. 
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ÉLARGISSEMENT D’UNE RAIFE PAR UN SÉLECTEUR D'AMPLITUDE 


Par P. MOATTI, 


Résumé. — On étudie l'élargissement d’un spectre d’énergie par sommation dans les canaux 


d’un sélecteur d'amplitude. 


On précise le cas de la déformation d’une courbe de Gauss lorsqu'on la construit par paliers 
rectangulaires. On donne les caractéristiques de la courbe de Gauss d’origine à partir de la courbe 


expérimentale. 


Application à l’étude des fluctuations de gain d’un photomultiplicateur. 


Abstract. — The author studies the broadening of an energy spectrum by summation in the 


channels of an amplitude selector. 


The author emphasizes the case of the deformation of a Gaussian curve when it is built with 


rectangular steps. 
one. 


He gives characteristics of the original Gaussian curve from the experimental 


Application to study of gain fluctuations of a photomultiplier. 


I. Position du problème. — Si l’on considère un 
spectre d’énergie, les impulsions d’une raie de 
forme quelconque, à leur entrée dans un sélecteur 
d'amplitude, ont une distribution en hauteur qui 
est définie par la densité de probabilité : 


où », est le nombre total d’impulsions par seconde 
relative à la raie, et v leur hauteur. L’histogramme 


MD) 
(D) 


0 LP D 


lite dl 


expérimental fourni par le sélecteur représente le 
nombre de coups par seconde : 


An(v, s) = n@[F{v + s) — F(v —s)] (1) 


F(v) = ['it) dv est la fonction de répartition en 


hauteur des impulsions ; 2s est la largeur d’un canal 
du sélecteur. 

La densité de probabilité expérimentale à la 
sortie du sélecteur devient donc : 
1 AV. ‘À 
g(v, s) rer [Fo +) = Fo —s)] 
au lieu de f(v). C’est l'effet d'intégration dans un 
canal par sélection d'amplitude, que nous nous pro- 
posons d’étudier ici (voir fig. 1). 


II. Élargissement d’une raie dans le cas général. 
— 19 CALCUL DE LA FONCTION DE RÉPARTITION 
VRAIE A L'ENTRÉE DU SÉLECTEUR. — Supposons les 
canaux juxtaposés ; on a, en sommant : 


(2) 


Il est donc toujours possible de connaître exac- 
tement la’ fonction de répartition à l’entrée du 
sélecteur. 

Pour avoir la densité vraie f(v), il faudrait diffé- 
rentier la courbe de discrimination, sans oublier le 
décalage des abscisses, le milieu du g°°e canal 
ayant pour abscisse v, — s et non pas Us. 


29 EFFETS D’EMPILEMENT ET DE TEMPS-MORT. — 
Ces effets ont déjà été étudiés par l’auteur [1], [2]. 
Par suite de l'effet d'intégration, qui a lieu simulta- 
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nément, le spectre d'énergie vraie {(v) à l’entrée du sélecteur est donné par : 


T8 Ho) do + (1 — en) il É h(v) dv (3) 


avec 


DANSE SANS ñ 
AE D diem (4) 


h(v) est la densité générale des impulsions empilées 5 T est le 
temps actif où le sélecteur fonctionne à la suite de l’arrivée 
d’une impulsion ; nr est le nombre d’impulsions déclen- 
chant le temps-mort 7 pendant la mesure de An. 


30 LARGEUR A MI-HAUTEUR DE LA RAIE. — Il est 
commode, et usuel, de mesurer la largeur d’une raie 
par sa largeur 2] à mi-hauteur ; cette grandeur est 
définie par l’équation : 


FD ts) Mots) 5 LF (vo + 5) — F(vo—s)] 


où v, est l’abscisse du sommet du pic. 


La largeur à mi-hauteur n’est donc fonction que 
de la largeur du canal. 


CDI 0,5 1 
Fix. 2. 


La courbe [(s) a la forme générale indiquée par la 
figure 2. Le calcul de /(s) et de ses dérivées est facile 
à partir de l’équation précédente. Nous donnerons 
simplement : 


(0) — Lo, demi-largeur de la courbe à l’entrée du sélec- 
teur 
l’(0) = 0 
” AUD) 
l’(0 EE ——— : 
REA 


Cette dernière relation montre que l’élargis- 
sement est minimum lorsque f/(v,) est petit et 
(vo + lo) grand, c’est-à-dire lorsque la courbe se 
rapproche d’une forme carrée. Le sélecteur tend 
done à «équarrir » un pic délié, et il ne modifie pas 
un pie carré. Celà est dû à l’intégration inévitable 
dans un canal de largeur non nulle. 


40 DÉGALAGE DU SOMMET. — En dérivant l’équa- 
tion (1) par rapport à v, on voit que le maximum de 
l’histogramme a lieu pour : 


AU nt à 
ce qui n’est vrai pour v — v, que si f(v) est symé- 
trique. 
Dans le cas général, le sélecteur crée donc un dépla- 
cement du sommet du pic, indépendant du dispositif 
de mesure, par le seul effet de la largeur des canaux. 


59 POSITION DES CANAUx. — L’équation (1) 
montre assez clairement qu’étant donné un spectre 
de densité f(v), l’histogramme An(v, s) ne dépend 
que de la largeur des canaux ; en particulier, leur 
juxtaposition exacte, et leur position par rapport au 
sommet du pic, n'ont aucune influence. Ce qui 
importe, c’est la connaissance de leur abscisse v en 
valeur absolue. 


ITT. Cas de la raie gaussienne. — C’est en parti- 
culier le cas du pic photoélectrique d’un rayon- 
nement y, et de certaines courbes de coïncidences. 
Il a déjà été examiné en partie [°] ; nous nous pro- 
posons d'apporter quelques précisions et appli- 
cations pratiques. 

La densité de probabilité à l’entrée du sélecteur 


est : 
a 1 __ @—0) 
1 Vie ml 


Elle est, à la sortie du sélecteur : 


(0, 9) = pe LED — m9 + 5) — Go — vo — s)] 


où Gest la fonction de Galton |” {(v) de, et co 
l’écart-type de la courbe de Gauss. 


19 LARGEUR A MI-HAUTEUR. — Cette largeur 2! 
est définie par l’équation : 


GR Pan ere an 7) 


No 4 


en prenant maintenant pour variable l'écart réduit : 
æ = (0 — v)/60 
La courbe 100(7 — 7,)/1, en fonction du rapport 


sl représente l'élargissement relatif à mi-hauteur 
de la courbe (voir fig. 2). 


29 LARGEUR AU QUART DE LA HAUTEUR. — Cette 
largeur 22 est définie par l’équation : 


De set 
On trouve que la loi : est très bien 


suivie quel que soit s, en particulier pour s — 0. 

Cela montre que l’histogramme expérimental ne 
s’écarte que très peu d’une forme gaussienne sur les 
flancs. C’est bien ce que l’on constate expérimen- 
talement. 


30 ÉGARTS-TYPES A LA SORTIE 6 ET A L'ENTRÉE Oy 
DU SÉLECTEUR. — L’écart-type de sortie 6 est 
donné par : 


ce qui conduit à : 


0 = 02 + 52/3 |. (5) 


Cette relation, à la fois rigoureuse et d’une grande 
simplicité, prouve qu’au moins au 2° ordre, on peut 
assimiler l'effet d'intégration par le sélecteur à l’addi- 
ton d'une variable aléatoire, de valeur moyenne 
nulle et d’écart-type s/V/3, indépendante de l’écart- 
vrai v — v,et s’ajoutant à lui. 

La courbe 100(6 — 6,)/o,en fonction du rapport 
sl (voir fig. 2) représente l’élargissement relatif 
de la raie gaussienne, exprimé en termes d’écarts- 
types. 

La courbe s/c en fonction de s/l (fig. 3) montre 
que l’écart-type © de la courbe de sortie s’écarte 
très peu de la loi linéaire : 


6 — 0,85 1 


c’est-à-dire d’une loi gaussienne. 

Par contre, la courbe s/6, en fonction de s/l 
(fig. 3) indique que o s’écarte beaucoup de 6,,. Celà 
signifie que la courbe expérimentale conserve son 
allure pseudo-gaussienne, même pour des largeurs 
de canal très grandes, mais qu’elle s'éloigne beaucoup 
de la raie gaussienne d’origine. 

En pratique, dans les domaines les plus divers, 
on rencontre souvent de telles courbes « pseudo- 
gaussiennes ». On risque de commettre une erreur 
importante en ne tenant pas compte de l’effet 
d'intégration dans les canaux. 
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297 
—_———— = OS = (os 
rare 
C} La droite représente 
1,5 la loi de Gauss ni 


0,85 —---- Y (So) 
 y{6) Fe 
» Î 
La droite représente ; 
Des Ja loi de Gauss 1 


0,75 


0,65 


49 ORDONNÉES RELATIVES AUX ÉCARTS-TYPES © 
ET 69. -— L’ordonnée relative au maximum de l’his- 


togramme est : 
S S 
G(x +2) à afs—è) 


AE) 


Pour æ = 1 et c/o, respectivement les courbes 
Y(60) et y(o), représentées figure 4, conduisent à la 
même remarque que précédemment. 


y(v, s) = 


50 CALCUL DE L'ACTIVITÉ TOTALE. — Le nombre 
de coups comptés dans le canal central est : 


Anmax = 270 G(s][60). 


228 


La courbe représentant Anñmax/r9 en fonction 
de s/l (fig. 5) est rectiligne dans tout le domaine 
pratique, de sorte que l’on a : 


| | Animax fo © 0,94 sh | (6) 
A TL max 
ES OS= OS 


0,5 


0,5 
FIG: 


et pour densité de probabilité maximum à la sortie 
du sélecteur : 


(0) & 0,4 Lt |. (7) 


La densité maxima est inversement proportionnelle 
à la largeur à mi-hauteur. C’est dire que la courbe 
de sortie «s’aplatit » assez vite lorsqu’on augmente 
la largeur des canaux. 

On peut tirer de la relation (6) un calcul très 
simple de l’activité totale n, relative à une raie 
gaussienne. Celà présente en outre l’avantage 
d'éviter une addition fastidieuse et moins exacte 
(par suite de la moins bonne précision statistique 
sur les flancs, de l’effet Compton sur l’aile gauche, 
de l'effet d’empilement sur l’aile droite, de l’eftet 
de temps-mort, etc...). 


E RS 1,06(1/5s) Anmax |. 


Cette formule a été très bien vérifiée expérimen- 
talement par l’auteur. Elle ne fait intervenir que le 
rapport //s, qui est bien caractéristique du phéno- 
mène étudié 1ci, et est facile à mesurer. 


(6 bis) 


60 INFLUENCE DU SÉLECTEUR DANS-LE CAS OU 
S — 69 — Supposons la demi-largeur du canal 
égale à l’écart-type 0, de la courbe gaussienne 
d'entrée. C’est très sensiblement le cas pressenti par 
R. Foucher, et autres [1]. En fait, on prendra 
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s — 0,85 [, au lieu de 0,8 Z,. Cela semble présenter 
les avantages suivants : la largeur de la bande 
centrale est exactement la largeur théorique de la 
raie à l’entrée du sélecteur, soit 20,, que l’on a 
ainsi sans calculs ; le réglage est très facile : le 
point d’écart s, c’est-à-dire situé sur la frontière 
du canal central, doit avoir pour ordonnée rela- 
tive 0,702. Or cette ordonnée y(s) est une fonction 
décroissante de la largeur du canal. 

On peut aussi utiliser la largeur à mi-hauteur, 
qui est une fonction croissante de s et de s/l. 
Pour s — 09%, on doit avoir: s/l —0,722. On 
compte 68 % des impulsions de la raie dans le 
canal central, nombre bien connu et élevé. L’élargis- 
sement relatif à mi-hauteur est de 17,5 %, c’est-à- 


F (x) se 
g (D) — 


HG 20: 


dire pas trop élevé. La figure 6 représente les den- 
sités de probabilité gaussienne et de sortie dans le 
cas où $ — 69: On peut y constater une autre pro- 
priété intéressante : pour l'écart s = 6, les 2 den- 
sités sont égales ; autrement dit, pour ce point privi- 
légié, l'intégration par le sélecteur n’a aucun effet. 


IV. Comparaison de l’élargissement d’un pie 
photoélectrique par un photomultiplicateur et par 
un sélecteur d’amplitude. -— 10 PosirIoN pu PRo- 
BLÈME. — Si le nombre moyen de photoélectrons 
par impulsion entrant dans le photomultiplicateur 
est 0, on sait que la dispersion au second ordre à 
l'entrée du P. M. est 1 /V/0. 

L'effet des fluctuations de l’émission secondaire 
est de multiplier cette dispersion par le facteur 


1/2 
1 Ho | 
| ji Bit 1] ? 


N°4. 


A 


N° & 


où u, et u, sont les 2 premiers semi-invariants du 
gain par étage [*]. Dans le cas d’un pie photo- 
électrique, la courbe de sortie du photomultipli- 
cateur est gaussienne, d’abscisse v, et d’écart-type 
Co et l’on a : 


A la sortie du sélecteur, l’histogramme a pour 
écart-type 6, et son abscisse centrale est toujours v, 
(voir IT, 40). Les courbes précédentes permettent de 
déterminer immédiatement l’élargissement rela- 
üf (o — 69)/09 à partir du rapport s/l (voir fig. 2). 

On peut alors comparer l’élargissement relatif de 
la raie dû au P. M. seul, soit 


Ho Le 
{ | 1, 
| F U(ts — 1) 


à l'élargissement relatif dû au sélecteur seul, 
soit (6 — 69)/50- 


20 RÉSULTATS EXPÉRIMENTAUX. — Le spectre 
d'énergie de lSAu a été tracé par l’auteur avec un 
tube E. M. I. 6097-B, un scintillateur INa, et un 
sélecteur rapide à 10 canaux dont la largeur était 
de 2 volts (C. E. N., Saclay). Les résultats (*) sont 
les suivants : 


(*) Les dispersions sont exprimées en termes théoriques, 


(o] L4 e [2 . 
de la forme = ; les expérimentateurs ont coutume de définir 
0 


; ; 21 : 
la résolution d’une raie par le rapport D» Ce qui donnerait 
0 


un ordre de grandeur double. 
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Énergie des y (MeV) 0,411 0,67 1,09 
Abscisse v, (volts) 33 5220 82 
sl 0,67 0,5 0,54 
Dispersion 
à l’entrée du P. M. (%) 2,7 24 1,7 
Dispersion 
à la sortie du P. M. (%) 3,4 3,0 1,8 
Dispersion 
à la sortie du sélecteur %) 3,8 3,2 1,9 
Élargissement 
relatif par le P. M. (%) 2? 22 22 
Élargissement (°/) 
relatif par le sélecteur 12 6 7 
30 Conclusions. — On voit facilement que le 


photomultiplicateur seul élargit le pic de 2 à 4 fois 
plus que le sélecteur seul. 

Si l’on tient compte : 

1) que l'effet du sélecteur peut être rendu aussi 
petit qu’on le veut (il suffit de réduire la largeur 
des canaux), et d’ailleurs facilement calculé quelle 
que soit cette largeur ; 

2) que, pour un scintillateur et une photocathode 
donnés, la dispersion à l’entrée du P. M. est une 
constante, que les progrès techniques rendent de 
plus en plus petite ; 


3) que le facteur de bruit [1 + He 


| est 
s (li — 1) . 
lié au earactère statistique de l’émission secondaire 
dans les photomultiplicateurs [4] ; il varie beau- 
coup avec la haute tension appliquée au P. M., le 
modèle de couche multiplicatrice, ete. ; 

4) que l’élargissement relatif par le photomul- 
tiplicateur seul est couramment de l’ordre de 
25 % ; on conçoit tout l’intérêt et la nécessité de 
l’étude des fluctuations de gain d’un photomultipli- 
cateur, liées à l’émission secondaire [°]. 


Manuscrit reçu le 22 octobre 1960. 
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ERRATUM 


Dans l’article de messieurs GRuN (J. B.), Sieskinp (M.) 
et NIKITINE (S.), J. Physique Rad., mars 1961, 22, p. 176, 
sur la figure 1, en haut à gauche, intervertir les mots rouge 


et jaune, et lire : 


fond continu jaune calculé 
fond continu rouge calculé. 
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utile d’effectuer une synthèse de cet ensemble assez disparate, afin d'établir les conditions pra- 
tiques d’emploi de ces sources. 

Après avoir décrit les modèles usuels de sources H. F., nous procédons à leur classement selon 
le mécanisme de la décharge (annulaire, linéaire, avec ou sans un champ magnétique auxiliaire). 
Le problème de l’extraction des ions est ensuite abordé ; nous passons en revue les principaux 
systèmes d'extraction, ainsi que les paramètres qui conditionnent l’intensité du faisceau extrait. 
Nous indiquons enfin quelle est la valeur de l’énergie moyenne des ions et la dispersion autour 
de cette valeur moyenne. 


Abstract. — Sources excited by high frequency oscillators (or H. F. sources) give high currents 
of positive ions, they are of current use in particle accelerators, ion microscopes, thermonuclear 
apparatus. Producing nearly monoenergetic ions, they are of value for mass spectroscopy. The 
papers about them, rather numerous, are related to very different matters, so that the data on 
these sources are found often in a more general study. We have considered that a synthesis of 
this rather varied material would be useful, in order to state the conditions of practical use of 
the sources. 

After describing the usual models of H. F. sources, we class them according to the mechanism 
of the discharge (annular, linear, with or without addition -of a constant magnetic field). Then, 
we describe the problem of the extraction of ions ; we describe the principal extraction systems, 
and we review the various parameters influencing the intensity of theextractedstream. Finally, 


TOME 22, AVRIL 1961, PAGE 230. 


we give the average ion energy, and the dispersion around this average value. 


INTRODUCTION 


Cockroft et Walton ont construit le premier accélé- 
rateur de protons en 1932. Depuis lors, les sources 
d'ions (protons, deutons, hélions, etc...) sont les auxi- 
liaires indispensables de toutes les machines accélé- 
ratrices, quels que soient leur importance et le prin- 
cipe de leur fonctionnement. L’arsenal des sources 
utilisables est important : tubes à rayons canaux 
(décharge sous pression élevée), sources à arc (basse 
pression ou capillaire), systèmes magnétiques à 
cathode chaude ou à cathode froide, sources à électrons 
oscillants, sources à émission secondaire, etc... Dans 
les meilleures conditions, de telles sources fournissent 
des faisceaux d’ions positifs dont l'intensité ne dépasse 
pas quelques centaines de microampères. 

Or, de nombreuses recherches nécessitent l'emploi de 
faisceaux ioniques très denses et très stables. Citons, 
par exemple, les études des interactions nucléaires des 


particules chargées, la production de neutrons rapides, 
Pinjection d’un faisceau dans un accélérateur de 
grande puissance, la création de plasmas dans les appa- 
reils thermonucléaires, etc... Les sources à excitation 
électrique de haute fréquence, que nous désignerons 
par «sources H. F. » dans ce qui suit, constituent la 
meilleure solution actuelle : elles produisent un cou- 
rant ionique très intense par une consommation de 
gaz réduite, le pourcentage des ions monoatomiques 
est très élevé, et ces ions sont presque mono-cinétiques ; 
très robustes, stables et d’un emploi facile, elles ne 
nécessitent pas un apport d'énergie élevé. 

Pour que ces qualités apparaissent nettement, les 
sources H. F. doivent être placées dans des conditions 
de fonctionnement bien définies et vérifiées avec beau- 
coup de soin, ce qui explique pourquoi les utilisateurs 
ne tirent pas toujours le maximum des sources qu’ils 
emploient. Nous croyons utile de grouper dans cette 
mise au point les conditions usuelles d’emploi (voir le 
tableau 1). 


ps nn - ape 
| CE D, OU 


LES SOURCES D'IONS A EXCITATION ÉLECTRIQUE 
TABLEAU I 


CONDITIONS D'EMPLOI DES DIVERSES SOURCES ÉTUDIÉES. LEURS CARACTÉRISTIQUES ÉLECTRIQUES. 


Les références bibliographiques sont classées dans l’ordre chronologique. 
(1) À = annulaire. ZL — linéaire. 
(2) H3 — champ longitudinal (œrsteds). 44 — champ transversal (œrsteds). 


} 
(3) Dimensions : diamètre (mm) ; longueur (mm). 
M Nature. Dimensions (mm). 
) 


(kilovolts) . = 


Rapporté aux conditions normales de température et de pression. 


RÉFÉRENCES [Ru-1] [Ba-3] [Ha-1] [Th-3] 
[Ha-2] [Th-4] 
Nature de la décharge (1) L A L ni 
Gaz utilisé Hydrogène Hydrogène Hydrogène Hydrogène 
Deutérium Deutérium Deutérium 
Intensité Z (milliampères) 0,4 0,5 0,4 05 
Ions atomiques (%) 60 57 (D+) 60 92 (H+) 
DAS ETES) 89 (D+) 
Dans la décharge (w.) 30 — 60 50 
Puissance 
Dans le générateur (w.) — 400 — — 
Fréquence (mégahertz) 180 SE 450 20 
Débit gazeux (cm*/heure) (5) 15 15 30 415 
Champ magnétique auxiliaire (2) Hy= 130 — He" 1000 — 
Pression dans la source (microns 
de mercure) 15 10 100 28 
Canal d'extraction (3) 3,175-12,70 2-12 DE 2,38-19 
Nature de la sonde Aluminium Aluminium Dural Dural 
Dural 
Chambre à décharge (4) Pyrex Pyrex Pyrex ou Pyrex 
g = 31,75 9 = 75 mm quartz 5 — 40 
MS mme ==" 11 0/rnm D — 1000 
1="10 
Tension d'extraction (kilovolts) 2,8 12 0 3 
Première tension d’accélération 
(kilovolts) 14 is 40 30 
RÉFÉRENCES [Mo-1] [Bu-1] [Sw-1] [Ba-1] 
[Sw-2] 
Nature de la décharge (1) L A Â A 
Gaz utilisé Hydrogène Hydrogène Hydrogène Hydrogène 
Deutérium 
Intensité Z (milliampères) 41,25 0,058 0,14 0,175 
Ions atomiques (%)) 92 69 80 61,5 
{ Dans la décharge (w.) 60 100 —- — 
Puissance 
Dans le générateur (w.) — — 100 à 200 200 
Fréquence (mégahertz) 100 25 25 22 
Débit gazeux (cm*/heure) (5) 6 — 1 — 
Champ magnétique auxiliaire (2) H, Hy (faible) (Hi 2 
Pression dans la source (microns 
de mercure) 15 — — 10 
Canal d’extraction (3) 1,6-6,35 4-12 0,5-3 180124 
Nature de la sonde Aluminium  Magnésium Acier Aluminium 
Aluminium 
Chambre à décharge (4) Pyrex Pyrex Verre Pyrex 
D 020 RNA) D 50 D —k( 
0=222,95 1E=0120 =" 110 LE 50 
Tension d'extraction (kilovolts) 5 5 5 29 
Première tension d’accélération 
100 50 14,9 


[Ne-1] 


A 
Hydrogène 


0,6 sans At 


2,0 avec H+t 
40 


150 


4-10 
Aluminium 


Pyrex 
Sphère de 2 L. 
0 à 10 
10 


[Be-1] 


L 
Deutérium 


0,25 


150 
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[AI-3] 


A 


Hydrogène 
Tritium 
0,2 


70 


25 


1,78-19,05 


Dural 
Aluminium 


Pyrex 

A — 040 

LS 
3,9 
30 
[Gi-1] 


L 
Hydrogène 


50 
D — 2 
Aluminium 


Pyrex 
g = 30 


1 


0 
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RÉFÉRENCES 
Nature de la décharge (1) 
Gaz utilisé 


Intensité Z (milliampères) 
Ions atomiques (%) 


Dans la décharge (w) 
Puissance : 
Dans le générateur (w) 


Fréquence (mégahertz) 

Débit gazeux (em°/heure) (5) 

Champ magnétique auxiliaire (2) 

Pression dans la source (microns 
de mercure) 


Canal d'extraction (3) 
Nature de la sonde 


Chambre à décharge (4) 


Tension d’extraction (kilovolts) 


Première tension d'accélération 
(kilovolts) 


RÉFÉRENCES 
Nature de la décharge (1) 
Gaz utilisé 


Intensité I (milliampères) 
Ions atomiques (°) 


Dans la décharge (w) 
Puissance 
Dans le générateur (w) 


Fréquence (mégahertz) 
Débit gazeux (cm°/heure) (5) 
Champ magnétique auxiliaire (2) 


Pression dans la source (microns 
de mercure) 


Canal d'extraction (3) 
Nature de la sonde 


Chambre à décharge (4) 


Tension d'extraction (kilovolts) 


Première tension d’accélération 
(kilovolts) 
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[Go-1] 


— 


A 


Hydrogène 
Deutérium 


0,63 


92 


45 
2,38-14 
Aluminium 


Pyrex 


250 


[Ci-2] 


L 
Hydrogène 


0,12 
70 


30 
Ht — 180 


37 
2,5-12 


Aluminium 


Pyrex 
0 = 25 


1,8 


k 


TABLEAU I (suite) 


[Ni-1] 


A 
Hydrogène 


2 


80 


10 
3-15 
Aluminium 


Pyrex 
9 = 40 


5 


50 


[Er-1] 


A 
Hydrogène 


0,7-0,8 


FI =—=150 


10 
2-15 
Aluminium 


[Eu-1] 


À 
Hydrogène 


15 


90 


55 
3,175-15,24 


[BI-5] 
L 
Hydrogène 


0,46 
75 


Hy = 500 


1 
1,6-1,2 
Aluminium 


Pyrex 
6 — 20 
l = 50 


2,9 


[La-1] 


L 
Hydrogène 


1,0 


75 


[Ci-3] 


A 
Hydrogène 


0,139 
85 à 91 


H}$—="60 


26 
2,5-12 
Aluminium 


Pyrex 
OÙ —=M2) 
P=M100 


2 


k 


[Mo-3] 


A 
Hydrogène 


0,1 
85 


100 


300 
50 

1,5 

H$ 


10 
1,1-15 
Aluminium 
Quartz 


Di 0 
— 90 


5) 


10 


[Co-4] 


A 
Hydrogène 


3-4 
Acier inox 


Quartz 

D —10150 

DER 
45 ï 


. Nou 


[AT-4] 


A 


Hydrogène 
Hélium 


0,3 (H+) 
0,2 (Het) 


0,03 Het + 


60 


60 
0,8-2,75 
Aluminium 


[Do-1] 


A 
Deutérium 


0,5 à1,4 


D — 3 
Aluminium 


N°4 


RÉFÉRENCES 
Nature de la décharge (1) 
Gaz utilisé 
Intensité Z (milliampères) 


Jons atomiques (%) 


Dans la décharge (w.) 


Puissance 


Dans le générateur (w.) 


Fréquence (mégahertz) 
Débit gazeux (cm°/heure) (5) 


Champ magnétique auxiliaire (2) 


Pression dans la source (microns 


de mercure) 
Canal d’extraction (3) 
Nature de la sonde 


Chambre à décharge (4) 


Tension d'extraction (kilovolts) 


Première tension d’accélération 


(kilovolts) 
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TABLEAU I (suite) 
[Sk-1] [Go-2] [De-1] [Ga-2] [Ho-1] [Se-1] 
[Se-2] 
A À ou L A A A 
Hydrogène Azote Hydrogène Hydrogène H, He Hydrogène 
Oxygène Deutérium Ne, A 
0,2 1 3 (H) 5-6 Æ 0,95 
1 (D) 
70 20 à 40 = æ _ 80 
de N+ ou O+ 
- æÆ 180 (H) = ce 130 
30 à 200 (D) 
80 — — 300 300 — 
27 25 30 40 29 42 
10 36 — 30 à 60 — 8,9 
— H; H+ ou H, Hi = "30 Fi =00 Ht 
H=20550 
45 0,8 20 — Best ANT 15405 
2,6-15,8 D — 10 De 3-10 0,9-2,5 1,5-7,5 
Aluminium — — Aluminium — — 
Pyrex — — Sphérique volume Pyrex 
D QE G) = 85 cmè QG —=2) 
DNS =n50 l = 90 
3 3 6 5 1,6 5 
== — 2 — 29,2 60 
RÉFÉRENCES [B1-3] [Ko-2] 
Nature de la décharge (1) A A 
Gaz utilisé Argon Hydrogène 
Deutérium 
Intensité Z (milliampères) (DE TEME ENT #) 2 
Ions atomiques (%) — — 
Dans la décharge (w.) — 300 
Puissance 
Dans le générateur (w.) 40 — 
Fréquence (mégahertz) 195 110 
Débit gazeux (cm°/heure) (5) 0,2 à 25 rues 
Champ magnétique auxiliaire (6) — Hy = 600 
Pression dans la source (microns de 
mercure) DTA 10 
Canal d’extraction (3) 3-6 avec 2,5-16 
diaphragmes 
Nature de la sonde Aluminium == 
Chambre à décharge (4) Pyrex Pyrex 
PE M'A) 
Ov 
Tension d'extraction (kilovolts) [A y 
Première tension d’accélération 
(kilovolts) 15 10 


I. MODÈLES DIVERS DE SOURCES H,. F. 


Les premières réalisations. — Getting [Ge-1] a été 
le premier à utiliser une décharge continue, obtenue 
dans un gaz par une excitation électrique de haute 
fréquence [12,8 méga-hertz (MHz), 100 watts], pour 


en extraire un faisceau d'ions positifs ; un tube de 
verre, de forme torique, contenait le gaz ; il était 
entouré par une bobine, traversée par le courant H. F. 
Les résultats furent peu encourageants : l’extraction 
des ions, faite par diffusion, conduisait à une intensité 
ionique très faible (130 microampères) pour un débit 
de gaz important. 
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La source d’Alfvën et Cohn-Peters [A1-1] (fig. 1) est 
plus proche des modèles actuels : une décharge H. F. 
(100 MHz) est créée entre les deux électrodes planes E; 
et E, placées dans un cylindre de verre, rempli d'hydro- 
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Fi1G. 2. — La première source H. F. de Thonemann [Th-2]. 


gène. L’ionisation est augmentée par les électrons 
secondaires venant des surfaces métalliques ; par 
contre, la recombinaison des ions sur ces surfaces est 
très importante, L’électrode négative E; a pour rôle 
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d'extraire les ions positifs. Cette source a l’avantage de 
fonctionner à basse pression (0,25 à 0,4 micron de 
mercure) mais elle manque de stabilité. 

La première réponse au problème de l’extraction des 
ions a été donnée en 1946 par Thonemann [Th-2] : un 
ballon sphérique de pyrex (2 litres) est entouré de 
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F1G. 3. — Source cylindrique de Rutherglen et Cole [Ru-1]. 


deux anneaux externes À et B, reliés à un oscil- 
lateur H. F. (60 MHz, 200 watts) (fig. 2) ; l’anode E, 
placée à un potentiel positif compris entre 0 et 20 kilo- 
volts et l’électrode d'extraction E,, à la masse, créent 
dans la source un champ électrique constant et 
assurent l’extraction des ions. 

Les sources H. F. ont pris en 1948 leurs formes 
actuelles : 

1. Rutherglen et Cole [Ru-1|, utilisant deux élec- 
trodes externes B (fig. 3) et un oscillateur de 180 MHz 
(300 watts), réduisent les surfaces de métal exposées à 
la décharge pour diminuer la recombinaison des ions. 
Ils ajoutent un champ magnétique axial, uniforme, 
(0 à 300 œrsteds). Pour l'hydrogène, le courant ionique 
est de 400 microampères, la proportion des protons 
de 60 %,. 

2. Le modèle de source H. F. de Thonemann, 
Moffatt, Roaf et Sanders [Th-3] est excellent (voir les 
fig. 4 et 15): un étranglement placé au sommet de 
l'enveloppe de la source (pyrex) permet de réduire au 
strict minimum la surface métallique exposée à la 
décharge ; de même, l’électrode d’extraction est en- 
tourée d’un cylindre de verre. Une bobine entourant la 
chambre à décharge est couplée inductivement à un 
générateur H. F. (20 MHz, 50 watts). Pour l'hydrogène, 
la proportion des protons est comprise entre 84 et 92 0, 
la dispersion en énergie étant inférieure à 50 eV. 

3. L'originalité de la source de Bayly et Ward [Ba-3] 
réside dans son système d’extraction : l’anode et l’élec- 
trode d’extraction, distantes de 1 cm environ, sont 
placées à la partie inférieure de l’ampoule de verre 
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Fic. 4.— Source H. F. cylindrique de Thonemann, Moffatt, 
Roaf et Sanders [Th-3]. Voir aussi la figure 15. 
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Fic. 5. — La source de Bayly et Ward [Ba-3]. 
Voir aussi la figure 14. 


(fig. 5). Hélas, la proportion des ions monoatomiques 
ne dépasse pas 50 %. 

4. La source de Hall [Ha-1}, [Ha-2], est nettement 
différente ; elle est formée par une chambre de faible 
volume (0, 5 em?) en pyrex (fig. 6 et 13), où la pression 
est beaucoup plus grande que dans les sources précé- 
dentes ; un oscillateur de 450 MHz (50 watts) y main- 
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tient une décharge sans électrodes. La formation des 
ions est facilitée par un champ magnétique axial. L’ex- 
traction s'effectue par diffusion à travers une ouverture 
disposée à la partie inférieure de la source. Une source 
analogue, mais de volume beaucoup plus grand, a été 
réalisée récemment par Bloom, Rytz et Staub [B1-5] : 
grâce à un système d'extraction par électrode, on a 
fortement diminué le débit de gaz, tout en obtenant 
un débit ionique voisin de 0,5 mA. 
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F1G. 6. — Source de faible volume réalisée 
par Hall [Ha-1], [Ha-2]. Voir aussi la figure 413. 


On trouvera dans les articles de Ward (Wa-1] et de 
Hoyaux et Dujardin [Ho-2] une synthèse détaillée des 
travaux que nous venons d’énumérer. 


Modifications apportées à ces modèles de sources. — 
Ces prototypes, et surtout la source de Thonemann 
(fig. 4) ont inspiré de nombreux modèles, assez voisins. 
Nous citerons les principaux, par ordre chronologique. 

1. Le système d’extraction de la source de Thone- 
mann a été modifié par Moak, Reese et Good [Mo-1]: 
le manchon isolant 7 (voir les fig. 7 et 16), en silice 
fondue, couvre le canal de l’électrode d'extraction 4 et 
isole de la décharge la cathode d’aluminium. Lorsque 
l’on applique le potentiel d'extraction, c’est-à-dire la 
tension constante entre l’anode et l’électrode d’extrac- 
tion, ce manchon fonctionne comme une anode vir- 
tuelle. Enfin, le disque de pyrex 2, à la partie supé- 
rieure de la chambre, empêche les électrons secon- 
daires venant de l’électrode d’extraction, de frapper 
l’anode. Allison et Norbeck [Al-4] ont simplifié ce type 
de sources en supprimant la sonde et son système de 
protection. 

2. Budde et Huber [Bu-1] refroidissent leur source, 
du type Thonemann, par un courant de pétrole, passant 
entre la source et une double paroi de verre. Lorsque la 
puissance dépensée est importante, 1l devient indis- 


. pensable de refroidir l'enveloppe de la source, par de 
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Pair comprimé, par exemple [La-1]. De plus, le refroi- 
dissement augmente, semble-t-il, [Ni-1], le pourcen- 
tage des ions monoatomiques contenus dans le faisceau 
extrait. 

3. La source de Swingle et Swann [Sw-1], [Sw-2] est 
une chambre de verre dont le fond plat (fig. 8) est 
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Fic. 7. — Modifications apportées par Moak, Reese et 
Good [Mo-1] à la source de Thonemann. Voir aussi la 
figure 16. 


4. Pyrex (5 — 25 mm). — 2. Bouclier de pyrex. — 
3. Anode. — 4. Sonde. — 5. Aiïimant permanent. — 
6. Disque d’acier. — 7. Silice. — 8. Entrée de gaz. 
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Fic. 8. — Source à fond plat de Swingle et Swann [Sw-1], 
[Sw-2]. 

1. Vers l’oscillateur 25 MHz.— 2. Lavite. —3. Cathode 
d'aluminium. #. Isolant. — 5. Champ magnétique 
transversal, — 6. Verre, — 7. Entrée de gaz. — 8. Porce- 

laine, 
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formé par un disque amovible. L’électrode d’extraction 
est un disque d’acier ou d'aluminium, percé axialement. 

4. Une source à fonctionnement réglable a été réa- 
lisée par Goodwin [Go-1], (fig. 9 et 17). Elle est ana- 
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\ F1G. 9. — Source de Goodwin [Go-1] 
à électrode d’extraction mobile. Voir aussi la figure 17. 
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F1G, 10, — Source à chambre mobile. Morozov [Mo-3]. 
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logue au modèle de Moak (fig. 7), mais l’électrode 
d'extraction est mobile selon son axe. De même, 
Morozov [Mo-3] a construit une source modifiable 
(fig. 10), dans laquelle le cylindre isolant de quartz (2), 
solidaire de la chambre (1), également en quartz, peut 
se déplacer par rapport à la cathode (3), fixe. 

5. Reifensweiler [| Re-3] a souligné l'importance du 
rôle que joue un diaphragme placé dans le canal de 
l’éléctrode d'extraction. Des diaphragmes de ce genre 
nous ont permis [B1-3], [Bl-4], d'obtenir des débits 
ioniques élevés (2 mA) pour une puissance H. F. 
dépensée faible (40 watts) et des débits gazeux réduits 
(0,2 à 25 cm par heure) [*], sans qu’il soit nécessaire 
d'employer un champ magnétique auxihaire (voir plus 
loin). 

6. En vue de leur montage sur les accélérateurs du 
Commissariat à l'Energie Atomique, Moreau et Vienet 
[Mo-2] ont construit des sources associées à des puis- 
sances H. F. plus grandes (400 à 500 watts). Le modèle 
définitif du C. E. A. [Co-4] peut fonctionner avec des 
oscillateurs de 1 200 à 1 500 watts. 

Cette énumération technologique de sources assez 
voisines est peut-être fastidieuse ; elle est cependant 
indispensable pour montrer que les sources H. F., 
variables dans leur géométrie et leur mode de couplage, 
forment à première vue un ensemble assez disparate. 
Cependant le fonctionnement de toutes ces sources 

- obéit à des lois précises, qu’il faut dégager. C’est ce 
que nous allons maintenant analyser. 


II. CARACTÉRISTIQUES DE LA DÉCHARGE H. F. 


Modes de couplage. — Pour exciter par courant H.F. 
un gaz sous basse pression (0,1 à 0,01 mm de mercure) 
— ce qui est le cas dans les sources du type Thonemann 
— on peut opérer de deux façons [Se-1], [Sk-1] : 


1. PAR COUPLAGE cAPACITIF à l’oscillateur H. F. 
(fig. 11a). L’ampoule qui contient le gaz est placée 


a) b) 


F1. 11. — Modes de couplage des sources H. F. 
(a) Capacitif. — (b) Inductif. 


entre deux anneaux métalliques formant les armatures 
- du condensateur d’un cireuit oscillant : c’est le champ 
électrique alternatif créé entre ces deux armatures qui 
provoque et entretient la décharge. Cette décharge 
est dite « linéaire », d’après la direction des forces 
électriques. 


2. PAR COUPLAGE INDUCTIF à l’oscillateur H. F. 
(fig. 11). L’ampoule est placée à l’intérieur de la self 
du circuit oscillant. Au champ magnétique alternatif 
créé par la bobine s’ajoute un champ électrique axial 


[*] Les débits sont toujours rapportés aux conditions 
normales de température et de pression (0 °C, 76 cm de 
mercure).' 
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créé par la différence de potentiel aux bornes de la 
self. La décharge est « annulaire ». 

3. Dans l’un comme l’autre de ces deux cas, la 
superposition d’un champ magnétique constant aug- 
mente considérablement l'intensité de la décharge 
[Ga-2], [Ko-1], [Ne-1], [Ne-2], [Sk-11. Dans ce qui suit, 
nous appellerons « champ transversal » un champ 
magnétique dont la direction est normale à l’axe de 
révolution des anneaux capacitifs ou à celui de la 
bobine de couplage. Un « champ longitudinal » est 
parallèle à cet axe. 


Structure de la décharge H. F. — L'émission d’élec- 
trons secondaires depuis les parois de la source joue un 
rôle fondamental dans l’amorçage et l’entretien de la 
décharge [Gi-1], [Sa-2], sans que l’on ait encore pu 
donner l'importance quantitative de ce phénomène. 
Par contre, il y a une recombinaison des ions sur les 
parois, surtout sur les électrodes, métalliques : le coef- 
ficient de recombinaison, qui est de 2.105 pour le 
pyrex et de 7.104 pour le quartz, atteint l’unité pour 
les métaux. C’est pourquoi il faut que les surfaces 
métalliques exposées à la décharge soient aussi faibles 
que possible. 

Cette décharge, dite «sans électrodes », est constituée 
par une région centrale fortement ionisée, identique- 
ment neutre, où le champ électrique global est nul ; les 
densités des porteurs de charges positives et négatives 
étant les mêmes, il s’agit d’une sorte de plasma. Ce 
plasma est entouré par une gaine ionique, d'épaisseur 
inférieure à 0,1 mm, qui lisole des champs électriques 
extérieurs (produits par les différences de potentiel 
entre les électrodes). Le potentiel du plasma est presque 
égal à celui de l’anode ; la forme du plasma est déter- 
minée par sa densité ionique, la géométrie des élec- 
trodes et les potentiels qui leur sont appliqués ; les 
systèmes d'extraction exploitent cette structure. 


Décharge linéaire. Son emploi est peu courant. 
Nous nous contenterons de citer quelques réalisations. 


1. SANS CHAMP MAGNÉTIQUE [Ci-1], [La-1]. — Peu 
de renseignements ont été publiés. Pour l'hydrogène, 
la source de Lacoste, Salmon et Watjsbrun/[La-1]donne 
une intensité inférieure à 1 mA pour une puissance H.F. 
de 660 watts (96 MHz). 


2. CHAMP TRANSVERSAL. — On peut citer [Gi-2] et 
[Th-2]. 


3. CHAMP LONGITUDINAL [Be-1], [Bl-5], [De-1], 
[Mo-1}, [Ru-1]. Voir la figure 3. La fréquence de l’oscil- 
lateur est généralement comprise entre 100 et 200 MHz ; 
Hall [Ha-1], cependant, utilise une fréquence de 
450 MHz. Le champ le plus intense (750 œrsteds) a été 
employé par Gabovich [Ga-3]. 


Décharge annulaire. — Très employée, elle a fait 
Pobjet de nombreuses études. 


1. SANS CHAMP MAGNÉTIQUE (voir les figures 4 et 5). 
Dans la plupart des cas, ces sources consomment une 
puissance comprise entre 100 et 200 watts [AI-3], 
[Ba-1], [Ga-1], [Ne-1], [Sk-1]. Certains auteurs em- 
ploient cependant des puissances plus grandes : Bayiy 
et Ward (400 watts) [Ba-3], Coutant, Prévot et Viénet 
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(700 watts) [Co-4], [Co-5], Van Dorsten et Spaa 
(500 watts) [Do-1], Eubank (300 watts) [Eu-1]. Pour 
des débits ioniques identiques, et grâce à l'emploi de 
diaphragme, ïl est possible d'utiliser une puis- 
sance H. F. beaucoup plus faible (40 watts) [BI-1]. 


2. CHAMP TRANSVERSAL (fig. 7 et 8). La puis- 
sance H. F. est également comprise entre 100 et 
200 watts, en général, [Bu-1], [Gi-3], [De-1], [Ne-1}, 
[Se-2], [Sw-1], mais parfois plus grande : Gabovich 
(300 watts) [Ga-1], Morozov (300 watts) [Mo-3], Niko- 
le et Crnilovic (300 watts) [Ni-1], Osetinski]j (300 watts) 
[Os-1]. Gabovich [Ga-1] propose de diriger le champ 
magnétique selon l’axe du canal d'extraction (fig. 25). 


3. CHAMP LONGITUDINAL. — Ces sources ont fait 
leur apparition en 1956 [Al-4]}, [Er-1]. Elles ont fait 
l’objet d’un certain nombre d’études, dont [Co-2] et 
* [Ko-2]. 


Conclusion. — D'une façon générale, dans des condi- 
tions identiques (puissance H. F., répartition des 
potentiels, pression, etc.….), les sources à décharge 
linéaire débitent un faisceau ionique moins intense que 
ceux des sources à décharge annulaire : ceci explique 
pourquoi on ne les utilise que peu. De plus, pour les 
décharges annulaires, les fréquences d’excitation (15 à 
70 MHz) sont beaucoup plus faibles que celles associées 
aux décharges linéaires (60 à 450 MHz). 

Notons que la comparaison des performances des 
divers modèles de sources doit être faite avec beaucoup 
de précautions, car le mode d’extraction des ions joue 
un rôle essentiel, que nous allons maintenant étudier. 


III. LE PROBLÈME DE L’EXTRACTION DES IONS. 


Généralités. — Cette question est liée à deux quan- 
tités : 

1. LE DÉBIT GAZEUX. — Il faut éviter un flux de 
gaz trop important dans la région où sont accélérés 
les ions. De plus, certains gaz sont coûteux : leur 
débit doit être réduit dans toute la mesure du possible. 


2. LA CHARGE D’ESPACE. — Elle limite l'intensité du 
faisceau ionique extrait. D’après Langmuir et Campton 
(Ree. Mod. Physics, 1931, 8, 237), cette intensité J du 
faisceau, en ampères par cm? de la surface émissive, 
est donnée par : 


J — Ke VVTM. 


Avec : d — distance séparant la surface du plasma 
de l’électrode d’extraction (cm). 
V — différence de potentiels (volts). 
M = masse atomique des ions (u. m. a.). 


Du côté de l’électrode négative (ou « sonde ») de la 
source, la limite du plasma est une surface équipo- 
tentielle, que les ions quittent normalement pour aller 
converger vers l’orifice d'extraction. 

La forme de cette surface équipotentielle conditionne 
la qualité de l’extraction. D’après Thonemann [Th-5] 
(fig. 12), cette surface varie en fonction de la différence 
de potentiel entre les électrodes (1) et (2) : pour une 
valeur convenable de V, la concavité de la surface 
équipotentielle favorise la focalisation des ions vers 
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l’ouverture de la sonde. Quand l’électrode d’extraction 
est un cylindre, la limite du plasma est rigoureusement 
sphérique [Sk-1], mais il est difficile de définir exac- 
tement sa position, les résultats expérimentaux restant 


V = OPTIMUM 


V=0 


F1. 12. — Formes possibles de la surface équipotentielle 
limitant le plasma. D’après Thonemann [Th-5]. 


incomplets. On peut seulement affirmer que la limite 
du plasma s’éloigne de la sonde quand le potentiel 
d'extraction augmente, ou que la puissance H. F. 
décroît [Wa-11]. 


Procédés d’extraction. — 1. EXTRACTION PAR DIFFU- 
SION des ions vers l’espace d’accélération, extérieur à 
la source. Cette diffusion s’effectue à travers un dia- 
phragme d’épaisseur aussi faible que possible ; les ions 
quittent la surface émissive avec l’énergie moyenne 


6 ar), k étant la constante de Boltzmann et T la 


température cinétique ; cette énergie est faible, et la 
charge d’espace a tendance à envoyer les ions sur les 
parois de l’orifice d’extraction, où ils se neutralisent ; 
c’est pourquoi le diaphragme doit avoir une épaisseur 
très faible. La pression dans la source est de l’ordre 
de 0,1 mm de mercure, et la consommation de gaz est 
élevé : il faut effectuer un pompage supplémentaire, 
au-dessous de la source et avant le système d’accélé- 
ration et de focalisation des ions. 


, 
7 


LIZZLL 


F1G. 13. — Détail du système d’extraction par diffusion 
de la source de Hall [Ha-1]. Voir aussi la figure 6. 


À titre d’exemple, on trouvera sur la figure 13, tirée 
de [Ha-1], le système d’extraction de la source de 
Hall (fig. 6). 


2. EXTRACTION PAR SONDE. — Cette méthode est 
presque universellement utilisée. La forme de cette 
sonde (électrode d'extraction) définit les performances 
de la source ; selon cette forme l'intensité du faisceau 
ionique peut varier dans le rapport de 1 à 40 [Th-3]. 
Dans la source, la pression est comprise entre 0,1 et 
0,01 mm de mercure, et il faut utiliser une ouverture de 
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diamètre assez faible pour réduire les fuites de gaz 
vers le système d'extraction, où règne un vide élevé. 
Les diamètres sont toujours inférieurs à 4 mm, ce qui 
rend inutile tout système de pompage supplémentaire. 
Selon la position de l’anode, on distingue deux sys- 
tèmes possibles d'extraction : 

— Système de Thonemann (fig. 4 et 14a) : l’anode 
est à la partie supérieure de la chambre d’ionisation ; 
la tension d'extraction est appliquée dans une région où 
la pression est relativement élevée (0,1 à 0,001 mm de 
mercure). C’est ce système qui conduit aux meilleurs 
résultats pratiques. 
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Fic. 14. — Les deux modes possibles d'extraction ionique 
par sonde : (a) Thonemann, (b) Ward. 


— Système de Ward (fig. 5 et 14b). Ce système a été 
déjà décrit plus haut : les deux électrodes, rigoureu- 
sement coaxiales, sont à la base de la chambre. L’anode 
constitue le support de la source et se trouve en contact 
avec le plasma. Ici, la tension d’extraction est appli- 
quée dans une région où la pression est faible (10—4mm 
de mercure ou moins) [Wa-1], [Mo-2], [Re-4]. 


Description de quelques systèmes d’extraction. — 
L’extraction par diffusion, pratiquement abandonnée, 
ne sera pas étudiée. 


1. DISPOSITIFS CLASSIQUES. — Le premier système 
de Thonemann (fig. 15), celui de la source de la figure 4 
[Th-3], se présente sous la forme d’un canal long de 
19 mm, dont le diamètre intérieur est 2,38 mm. La 
cathode (1) est en dural, ce qui évite une émission 
importante d'électrons secondaires. Le cylindre de 
verre (2), isolant, fait saillie dans la chambre ; il 
sépare la surface métallique du plasma, dont la limite 
se fixe à l’entrée de (2). 

Nous avons indiqué plus haut que Moak, Reese et 
Good [Mo-1] ont modifié ce dispositif en introduisant 
un manchon isolant de silice, coaxialement au canal 
de la cathode. On voit sur la figure 16 que ce manchon 
recouvre entièrement le canal. Il est essentiel que les 
axes de ces deux pièces coïncident exactement [Sk-1|, 
[Th-5]. D’après Moak, le courant extrait est réduit 
de 80 % quand le sommet de la cathode est de forme 
conique (ouverture : 1400) et de 95 % lorsqu'il est 
hémisphérique. Il faut donc que le bord supérieur du 
canal soit plat [Mo-1]. 

Le manchon de silice joue un rôle essentiel parce 
que son potentiel est approximativement égal à celui 
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du plasma : la limite du plasma se maintient à l’entrée 
du manchon [Mo-1], [Th-5], qui définit la forme de la 
surface émettrice d'ions et gouverne ainsi la topo- 
graphie du faisceau extrait [Re-3]. Servant d’anode 
virtuel [Sk-1], ‘ce manchon isole le plasma de la 
cathode [Er-1] : on évite ainsi les claquages, qui sont 


Fi. 15. — Détail du système d’extraction de la source 
de Thonemann [Th-3]. Voir la figure 4. 
1 Cathode: 9 — 2,33 mm, {—19mm.— 2""Partie 
isolante (verre ou quartz). 


TOLÉRANCE SILICE 


FiG. 16. — Détail du système d’extraction 
de la source de Moak [Mo-1]. Voir la figure 7. 


dus, en général, à la souillure de la surface de cette 
électrode par une condensation de vapeur d’huile rési- 
duelle du pompage, aux endroits les plus exposés au 
bombardement ionique [Er-1]. 
Dans le système de Ward, l’anode forme un cylindre 
devant la sonde et remplit le rôle du manchon isolant, 
La position du bord supérieur du canal d’extraction 
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par rapport au bouclier isolant (disque de verre ou 
manchon) constitue un paramètre critique du fonc- 
tionnement de la source [A1-3], [Eu-1], [Ga-3], [Re-1], 
[Re-3]. D’après Goodwin [Go-1], la profondeur hk de 
l’électrode d’extraction dans le manchon isolant doit 
être égale au diamètre externe de cette électrode : on 
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Fi. 17. — Détail du système d’extraction 
de la source de Goodwin (4 — 4 mm) [Go-1]. Voir la figure 9 


trouvera sur la figure 17 le système d’extraction de la 
source de Goodwin. 

Un choix convenable de diamètre D du manchon 
(fig. 19) améliore considérablement la focalisation 
[Mo-1] ; c’est ce qu'ont montré Gabovich et Pjankov 
[Ga-3] : d’après la figure 18, obtenue pour le système 
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F1G. 18. — Variation de l’intensité Z du faisceau ionique 
extrait en fonction du diamètre d du canal d’extraction 
[Ga-3] pour la géométrie de la figure 19. 


ATEN 


de la figure 19, une valeur de 5 mm correspond aux 
conditions optimales. 

Pour la disposition géométrique de la figure 20, 
Serbinov [Se-2] a fait une étude expérimentale très 
détaillée de l'influence des dimensions du système 
d'extraction sur les performances de la source. La foca- 
sation optimale du faisceau d’ions dans le canal de la 
cathode correspond à ZL/d = 5et Dfd = 2, 
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(0,6 D < I < 0,8 D). Le rapport D/d influe sur la 
valeur de la tension d'extraction pour laquelle le 
rapport du courant extrait au courant ionique tombant , 
sur la cathode est maximal ; cette tension optimale 
croît avec (D/d) TAI-3], [Eu- 1], [Se-2]. Plusieurs 
auteurs ont mis en évidence l’existence de cette tension 
optimale [B1-31, [Mo-3], [Th-5], [Th-6], [Wa-1]. Enfin, 


Fic. 19. — Système d’extraction 
de Gabovich et P’Jankov [Ga-3]. 


F1G. 20. -— Système d'extraction 
étudié par Serbinov [Se-2]. 


quand on utilise un champ magnétique auxiliaire, la 
position de la sonde dans le manchon isolant devient: 
fonction de la distribution du champ [Ni-1]. 

Le nombre élevé des paramètres fait que, dans l’état 
actuel des choses, les dimensions à adopter pour tout 
nouveau système d'extraction doivent être déterminées 
expérimentalement en tenant compte des diamètres du 
canal et du bouclier isolant ; Reifenschweiler [Re-3], 
[Re-4], a bien énoncé des lois de similitude concernant 
les dimensions du canal, mais elles ne sont pas encore 
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étayées par des vérifications expérimentales satis- 
faisantes. 

Pour compliquer encore les choses, la focalisation du 
faisceau dépend également de l’état superficiel de la 
paroi interne du col de l’isolant : la pulvérisation sur 
cette paroi d’une couche de métal venant de la cathode 
détériore les qualités de la lentille ionique que forment 
la limite du plasma et l’électrode d’extraction. Ce dépôt 
modifie en effet la distribution du potentiel électrique 
sur l’isolant [Mo-3], [Sk-1]. 


2. EMPLOI D'UN DIAPHRAGME DANS LE CANAL D’EX- 
TRACTION. — Les études de la propagation des ions 
dans le canal d’extraction, effectuées par Reifens- 
chweiïler [Re-1] constituent le point de départ de 
l'utilisation d’un diaphragme. Cet auteur utilise une 
cathode (d — 3 mm, Z — 19 mm) divisée en trois 


sections électriquement isolées (fig. 21), ce qui permet 
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Fic. 21. — Électrode d’extraction coupée en trois sections. 
Reiïfenschweiler [Re-1]. 


de mesurer le courant ionique capté par chacune de ces 
sections. Reifenschweiler constata ainsi qu’il existe une 
région de resserrement du faisceau dans la partie supé- 
rieure du canal. Plaçant un diaphragme dans le canal 
(fig. 22a), Reifenschweïler [Re-3] trouva que le débit 
ionique est augmenté dans un rapport voisin de 1,5, 
et que la consommation de gaz est divisée par deux. 
L'introduction d’un diaphragme dans les dispositifs 
d'extraction du type Ward [Re-4] conduit à des résul- 
tats analogues (fig. 22b). Bien qu’encore assez peu 
utilisés, [Do-1], [Ho-1], les diaphragmes permettent 
de réaliser des sources très stables et à fort débit 
ionique [B1-3], [BI-4]. 
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FiG. 22. — Utilisation d’un diaphragme : 


a) Dans un système d’extraction type Thonemann. 
b) Dans un système d’extraction type Ward-Reifens- 
chweiler [Re-3]. 


3. SOURCES A GRAND DÉBIT IONIQUE. — Il faut que 
la focalisation des ions soit excellente pour que l’on 
puisse parvenir à un débit élevé : la surface limite du 
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Fic. 23. — Système d'extraction 
de Thonemann et Harrison [Th-6]. 


GRAPHITE ACIER INOXYDABLE ALUMINIUM QUARTZ 


F16. 24. — Système d'extraction de Coutant, Prévot et 
Viénet pour source à débit élevé [Co-5]. 


plasma doit donc être hémisphérique. Pour réaliser 
cette condition, Thonemann et Harrison [Th-6] uti- 
lisent des électrodes planes (fig. 23) : S, est non conduc- 
teur, S, est la cathode. Une telle disposition donne 
des sources très stables, dans lequelles on n’observe pas 
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de phénomènes disruptifs entre la cathode et le plasma 
[AÏ-4], [B1-1], [Co-4], [Ga-3]. Moreau et Viénet [Mo-2] 
ont adapté ce dispositif à un système du type Ward, 
ce qui assure un isolement excellent entre le plasma 
et l’anode. Coutant, Prévot et Viénet [Co-5] ont réalisé 
un système analogue (fig. 24) pour des potentiels 
d'extraction très élevés (10 à 15 kilovolts) : la source, 
de puissance élevée (700 watts) fournit 7 mA pour une 
consommation de gaz de 75 em® par heure (conditions 
normales de pression et de température). 


4. ROLE DE LA CHARGE D’ESPACE. — L'existence 
d’un canal d’extraction a pour conséquence que la 
charge d’espace joue un rôle important. D’après Span- 
genberg, cette charge d’espace limite le courant à 
travers le canal à une valeur maximale proportionnelle 
à (d/L)?, d'étant son diamètre L sa longueur (voir la 
fig. 20). Pour diminuer les effets de cette charge d’es- 
pace, on pourrait augmenter d, mais la formule de 
Knudsen montre que le débit de gaz passant dans le 
canal est proportionnel à (d%/L) ; il est donc préférable 
de diminuer la longueur ZL et de placer un diaphragme 
de diamètre convenable [B1-4] dans le canal d’extrac- 
tion. 

A la sortie de la source, la charge d’espace produit 
une divergence du faisceau et il faut réaliser une foca- 
lisation ultérieure. On place en général au-dessous du 
canal une lentille électrostatique à champ intense, 


IV. PARAMÈTRES GOUVERNANT L’INTENSITÉ 
DU FAISCEAU EXTRAIT. 


Nous venons de faire une analyse globale de la 
décharge H. F. et du processus d'extraction. Nous 
basant sur les résultats publiés, nous allons maintenant 
étudier l'influence individuelle des divers paramètres 
de la décharge. 


Potentiel d’extraction — Nous avons vu que ce 
potentiel définit la surface limite du plasma, donc la 
focalisation des ions dans le canal d’extraction : sa 
valeur doit être choisie avec soin. Cette valeur varie 
avec la pression [Fo-1],[La-1]et avec la puissance H.F. 
dissipée dans la source [Co-4], [Go-f]. Lorsque l’on 
place un diaphragme dans le canal d’extraction, la 
valeur de ce potentiel joue un rôle très important [B1-2]. 
Un champ magnétique constant, transversal [Ga-2] ou 
longitudinal, [Ko-2], produit les-mêmes effets qu’une 
augmentation du potentiel d'extraction. 


Champ magnétique auxiliaire. — Des chercheurs 
soviétiques [Ga-2], [Ko-2], ont mesuré l’augmentation 
du courant ionique résultant de l’adjonction d’un 
champ auxiliaire : pour une puissance H. K, de 
300 watts, consommée dans la décharge, le débit peut 
atteindre 5 mA ; l'intensité de la décharge est considé- 
rablement augmentée [Be-1], [Ne-1}, [Ne-2]. 


1. CHAMP TRANSVERSAL. — Pour une fréquence 
H. F. donnée, l’intensité du faisceau extrait passe par 
un maximum très marqué en fonction de l’intensité du 
champ [Ko-1]. Il se produit un effet de réso- 
nance [Ga-2], [Mo-3], [Se-2], car la fréquence H. F. 
correspond alors à la fréquence f du cyclotron [Ha-1] : 


NA) Pntioe 
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L'intensité du champ correspondant est généra- 
lement faible et comprise entre 10 et 60 œærsteds. 
L'augmentation de l'intensité du champ transversal 
produit une concentration de la luminosité de la dé- 
charge vers l’axe de la source [Ci-2], [Ni-1]. Les cher- 
cheurs soviétiques [Mo-3], [Os-1] effectuent la régu- 
lation de l'intensité de la décharge en modifiant l’in- 
tensité du champ transversal. 


2. CHAMP LONGITUDINAL. — La densité ionique de 
la décharge diminue à mesure que l’on se rapproche de 
l’axe de la source [Ga-2]. On peut établir un champ 


uniforme dans la totalité du volume de la chambre . 


[Ho-1], [Ru-{], ou se contenter d’un champ local à la 
base de la source : le champ agit en effet sur la position 
de la limite du plasma près de la cathode et réduit, 
dans le canal, la dispersion due à la charge d’espace 
[AI-4], [Mo-1]. L’augmentation de l'intensité du fais- 


-ceau extrait ne vient pas d’un effet de résonance et les 


champs employés sont beaucoup plus élevés (300 à 
1 000 œrsteds) que dans le cas d’un champ transversal. 

Gabovich [Ga-2] a étudié avec beaucoup de soin 
linfluence du champ magnétique. Il a montré que, 
dans le cas d’un couplage capacitif, le champ doit être 
longitudinal mais que, pour un couplage inductif, il 
faut utiliser un champ transversal, c’est-à-dire normal 
à l’axe de la bobine d’induction. Il a été ainsi conduit à 
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F16. 25. — Source à décharge annulaire 


et champ transversal. Gabovich [Ga-2]. 


réaliser une source (fig. 25) à décharge annulaire, dans 
laquelle le champ transversal est parallèle de l’élec- 
trode d’extraction. 


Puissance H. F. dissipée dans la source. — L’aug- | 
mentation de la puissance H. F. intensifie la décharge : | 


on le remarque à l’accroissement de la luminosité du 
gaz [Co-4], [Do-1], [Eu-1]. Cependant, en particulier à 
cause de la charge d’espace, le courant extrait passe 
par un maximum en fonction de cette puissance [Co-4], 
[Ho-1]. Pour des performances pratiques données, un 
champ magnétique constant permet — nous l’avons 
vu — de diminuer la puissance H. F. [Ho-1], [Sw-1]. 


Pression dans la source. — Liée au débit de gaz, elle 
constitue un paramètre fondamental du fonction- 
nement : le courant ionique extrait passe en effet par 


; 
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un maximum en fonction de la pression [B1-4], [Ci-2], : 


[Eu-1], [Ho-1]. 
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Pour des diaphragmes de différents diamètres intro- 


duits dans le canal d’extraction, nous avons prouvé . 


[BI-4] qu’il existe une pression optimale, et une seule 
(fig. 26). Dans une source remplie d’argon sous cette 
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Fic. 26. — Variation du courant ionique extrait avec la 
pression, pour des diaphragmes de différents diamètres 
placés dans le canal d’extraction [B1-4]. 


pression optime (21,5 microns de mercure) la consom- 
mation de gaz varie de 25 cm par heure pour un dia- 
phragme de 3 mm de diamètre, à { em$ par heure pour 
un diaphragme de 0,8 mm de diamètre. L’intensité du 
courant ionique correspondant varie de 1,7 mA à 
1,05 mA. 


Fréquence de l’oscillateur. — Les utilisateurs ont 
souvent le tort de ne pas prêter d’attention à l’in- 
fluence de ce paramètre. Nous ne pourrons done donner 
que quelques indications assez fractionnaires. 

— Dans le cas d’un couplage inductif, les fréquences 
sont généralement de l’ordre de 100 MHz [Ba-1}, 
[Go-1], [Mo-1]. Ces valeurs assez peu élevées ont 
l'avantage de réduire le chauffage diélectrique des 
parois de la chambre [Sk-1] ; par contre on est conduit 
à employer des sources de grand volume [Gi-1], [Mo-1]. 

— Il est important de noter que, pour un couplage 
inductif, l'absorption d'énergie dans le plasma est fa- 
vorisée par une augmentation de la fréquence [Ni-1}, 
[Se-2]. Nous pensons donc que, chaque fois que lon 
veut employer un oscillateur de faible puissance, il 
faut, par contre, choisir une fréquence élevée [BI-21, 
[Ba-3]. 

— En présence d’un champ magnétique transver- 
sal, nous avons vu plus haut que le courant ionique 
extrait passe par un maximum quand la fréquence 
correspond à celle du cyelotron [ Ha-1]. 

— Enfin, on a remarqué qu’il faut augmenter avec 
la fréquence l’intensité du champ magnétique, quand 
il est longitudinal [Ba-1], [Co-4]. 
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V. CONSTITUTION DU FAISCEAU EXTRAIT. 


Cette question est d’une très grande importance 
pour les accélérateurs de particules ét pour la spectro- 
scopie de masses. 


Proportion des ions monoatomiques. — On a surtout 
étudié l'hydrogène et le deutérium, par analyse magné- 
tique. 


1. RÉSULTATS OBTENUS. — Accélérant des ions d’hy- 
drogène sous 15 kilovolts et faisant l’analyse par un 
champ électrique et un champ magnétique (400 œrs 
teds) croisés, Rutherglen et Cole [Ru-1] ont obtenu 
60 % de protons dans le faisceau extrait de la source 
représentée sur la figure 3. Avec la source de la 
figure 5, Bayly et Ward [Ba-3] ont obtenu 51 % de 
protons. Par contre, avec la source de la figure 4, 
Thonemann, Moffatt, Roaf et Sanders [Th-3] ont 
obtenu 89 % de deutérons pour une pression de deuté- 
rium de 25 microns de mercure dans la source, et 92 % 
de protons avec une pression d'hydrogène de 26 mi- 
crons de mercure dans la source : la réduction des 
surfaces métalliques exposées à la décharge explique 
ce pourcentage élevé [Ba-1], [Mo-11, [Ni-1]. 

Pour les sources actuelles, le pourcentage des ions 
monoatomiques est compris entre 70 et 92 %, selon 
les modèles (voir le tableau 1). 


2. PARAMÈTRES DIVERS. — Ils sont nombreux et 
ont des effets assez complexes. 

a) Recombinaison des ions dans la chambre à décharge. 
— Nous avons eu l’occasion de souligner l'importance 
de la nature des parois [Ha-2], mais leur état de pro- 
preté superficielle constitue un élément prépondérant 
[Gi-1], [Mo-1]. Cilensek, Cvelbar et Ramsak ont cons- 
taté [Ci-2] que le nettoyage de la chambre par une 
solution diluée d'acide fluorhydrique (5 à 10 %) aug- 
mente fortement la proportion des protons. Skurmik 
[Sk-1] obtient 90 % de protons après un nettoyage par 
une solution d’acide fluorhydrique, alors que, après 
nettoyage par de l’éthanol puis de l’eau distllée, cette 
proportion ne dépasse pas 50 %. La structure super- 
ficielle des parois se modifie lentement lors du fonc- 
tionnement de la source, et le pourcentage des ions 
monoatomiques diminue après 70 heures de fonction- 
nement, environ | Kh-1], [Sk-1]. 

b) Le refroidissement régulier des parois de la source 
assure une bonne stabilité à la proportion d'ions mono- 
atomiques extraits [Ni-1}, [Sw-1]. 

c) Gette proportion augmente avec la puissance H. F 
[Ha-2], [Re-3], [Sk-11. 

d) Elle augmente aussi si l’on ajoute un champ 
magnétique constant, transversal [Ne-1], [Ni-1], ou lon- 
gitudinal [Sk-1]. l 

e) La pression, enfin, joue un grand rôle, [Ba-3], 
[Ha-2], mais on manque de données quantitatives con- 
cernant son action. Notons que Cilensek, Cvelbar et 
Ramsak [Ci-2] considèrent que l'influence d’un champ 
magnétique ou de la pression serait nulle. 

Î) L’emploi d’un mélange de gaz permet d'augmenter 
le pourcentage des ions monoatomiques. Bayly et 
Ward ont expérimenté des mélange [Ba-3] d'hydrogène 
et d’hélium, d’argon, de vapeur d’eau ou d'oxygène ; 
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le débit ionique est plus faible, mais le pourcentage des 
protons est augmenté. Par contre, d’après Cilensek, 
Cvelbar et Ramsak [Ci-3], la présence de vapeur d’eau 
diminuerait le pourcentage des protons ; tout semble 
dépendre, en fait, des proportions relatives des deux 
constituants gazeux. Arnold [Ar-1] a constaté que, 
pour un mélange d’hélium et d'hydrogène (12,5 % en 
proportions atomiques), le pourcentage des protons est 
le même que pour de l'hydrogène pur ; un mélange 
d’hélium et de tritium conduit au même résultat en ce 
qui concerne les tritons. 


Énergie des ions extraits. — On s’est rendu compte 
très vite [Ha-2] que les énergies des ions doivent être 
assez homogènes : d’après Bayly et Ward [Ba-3], la 
dispersion ne dépasserait pas 100 eV. Les recherches 
concernant cette question sont encore peu nombreuses, 
et nous ne pouvons présenter que des résultats frag- 
mentaires. 


1. VALEUR MOYENNE DE L'ÉNERGIE. — La presque 
totalité des ions est créée à l’intérieur du plasma. 
D’après la théorie de Langmuir, leur énergie corres- 
pond donc au potentiel du plasma, supérieur de 
quelques volts à celui de l’anode. Ce point a été vérifié 
expérimentalement [La-1]. Les électrons secondaires 
venant de la cathode et les ions émergeant du plasma 
produisent une ionisation dans la couche obscure, voi- 
sine de la cathode ; l’énergie des ions ainsi créés corres- 
pond au potentiel du point où ils sont formés. Ils ne 
constituent heureusement qu’une fraction très faible 
({ à 2 %) du nombre total des ions (Er-1). 

La première mesure directe du spectre d’énergie du 
faisceau a été faite par Erô en 1956 [Er-1], avec un 
analyseur électrostatique placé directement à la sortie 
d’une source à décharge annulaire : l’énergie moyenne 
est supérieure de 50 à 100 volts à la tension continue 
appliquée à la décharge. D’après Erô, cette énergie 
moyenne pourrait dépasser la tension anodique d’une 
quantité égale à la valeur maximale de la tension H. F. 
appliquée à la bobine de couplage. Sur l'hydrogène, 
Forst a effectué en 1958 [Fo-2] une séparation magné- 
tique : aux ions dont l’énergie correspond au potentiel 
de l’anode s’ajoutent des ions d’énergie moyenne plus 
faible, qui proviendrait de la variation du potentiel 
superficiel sur les paroïs de la chambre. Holz et L6b 
[Ho-1] enfin, ont montré que des ions se forment au 
voisinage du canal d'extraction. 


2. DISPERSION AUTOUR DE LA VALEUR MOYENNE. — 
Les résultats sont assez variables, car la dispersion 
dépend des conditions d'emploi de la source : 40 à 
50 eV [Th-4], 200 eV [Ci-1], 60 à 70 eV [Mo-1], ete. 

Ero a obtenu [Er-1] une dispersion inférieure à 
50 eV pour des ions ayant une énergie de 4 000 eV, 
correspondant à une tension d’extraction de 4 000 volts. 
La dispersion diminue avec le champ H.F., le courant 
anodique et le potentiel d'extraction [Er-1], [Sk-1]. Ce 
même auteur, utilisant un analyseur électrostatique 
(1279) à grande vitesse, a enregistré en 20 micro- 
secondes la totalité du spectre [Er-2], avec une grande 
précision : la dispersion obtenue est voisine de 50 eV. 
L'emploi d’un champ H. F. de démodulation permet 
de parvenir à une dispersion inférieure à 10 eV. 

Forst [Fo-2], avec de l'hydrogène, obtient une dis- 
persion de 70 eV pour les ions dont l’énergie moyenne 
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correspond au potentiel de l’anode, et une dispersion 
voisine de 1 000 eV pour le groupe d'énergies plus 
faibles (voir plus haut). Il a noté que, si la dispersion 
dépend de la puissance H. F., de la pre et du 
potentiel anodique, elle est indépendante de l’état 
superficiel des parois et de la géométrie de l’électrode 
d'extraction. Enfin, la dispersion est d’autant plus 
faible que la masse des ions est plus grande [Fo-2]. 


VI. APPLICATIONS DES SOURCES H. F. 


Ces sources, dont la durée de vie est comprise entre 
100 et 200 heures de marche [Co-4], [Go-1], [Kh-1], 
[Re-4], trouvent de nombreuses applications. 


Accélérateurs de particules. — Fournissant une pro- 
portion élevée d’ions monoatomiques, elles sont très 
utilisées dans les accélérateurs de particules, où on les 
préfère en général aux sources du type Penning. 


1. PRODUCTION DE PARTICULES DE FAIBLE ÉNERGIE. 
— Tous les modèles de sources décrits plus haut sont 
employés dans les appareils de ce genre, généralement 
du type Cockroft-Walton. Les publications abondent : 
[AI-4], [Bi-1], [Da-1], [Eu-1], [Go-11, [Ha-1], [Ko-21, 
[Ma-2], [Mo-3], [Se-1], ete. 

La source de Moak, par exemple, équipe l’accélé- 


rateur de 0,3 MeV du Centre d’Oak-Ridge ; le faisceau 


de protons, d'ions (,Heë#)+ ou (,He4)+, est focalisé 
sous la forme d’un spot de 3 mm de diamètre, à 
2 mètres de la base de la source [Mo-1]. 

Dans les générateurs de neutrons rapides utilisant 
les réactions (D-D) ou (D-T), ces sources permettent 
d'obtenir un faisceau très intense [Ba-1], [Be-1], [Co-11, 
[De-1], [Ko-2], [Ma-3], ete. A titre d'exemple, l’accé- 
lérateur de Dorsten et Spaa [Do-1], de 300 à 1 000 kilo- 
volts, produit plus de 10 milliards de neutrons par 
centimètre et par seconde, en utilisant une source 
de 30 MHz et 500 watts, débitant 0,5 à 1,4 mA de 
deutérons ; ce flux neutronique est identique à celui 
d’une source de (Ra-Be) de 1 000 curies. 


2. ACCÉLÉRATEURS DE VAN DE GRAAFF. — Les 
sources H. F, y sont couramment employées [Go-1], 
[Gr-1], [Sk-11, [Sw-11]. Magnac-Valette et Liess [Ma-1] 
ont conçu une source consommant 5 cm par heure 
d’un mélange d'hydrogène et de tritium (0,5 cm® de 
tritium pour 25 cm° d'hydrogène) spécialement des- 
tinée à un accélérateur de Van de Graaff, 


3. DOMAINE DES HAUTES ÉNERGIES. — Des sources 
fonctionnant en impulsions permettent d’injecter des 
bouffées de particules dans le système d’accélération ; 
la source de Malkiel et Sukhanov [Ma-5] produit, pen- 
dant une microseconde, 100 à 200 mA de protons ou 
de deutérons, avec une fréquence de répétition de 4 000 
par seconde, Cette source (30 MHz) consomme 
30 watts en moyenne ; son rendement en protons ou 
en deutérons est compris entre 80 et 90 % ; elle débite 
30 em® par heure, environ. 


Spectroscopie de masses. — Ces sources sont appli- 
cables à la spectroscopie de masses [B1-1], où il est 
important de disposer d’un faisceau ionique d’inten- 
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sité très stable ayant une très faible dispersion d’éner- 
gie [Ha-2], [Se-2]. 


Microscopie ionique. — Des sources H. F. ont équipé 
le microscope protonique par transmission de Magnan 
et Chanson [Ma-4]. 


Production d'ions négatifs. — En inversant les pola- 
rités, Khirnij [Kh-1] a produit des ions (H)— (20 micro- 
ampères) et (D) (15 microampères) dans le capillaire 
d’une source H. F.;le débit de gaz est important 
(28 em$ par heure). Il en est de même pour la source de 
Phillips et Tuck [Ph-21, qui produit des ions (H)— 
avec un débit gazeux de 20 à 25 cm$ par heure. 


Produetion d'ions lourds. — Cornides, Roosz et 
Siegler [Co-3] ont construit une source d'ions (Li)+ du 
type Thonemann ; l’échantillon est chauffé extérieu- 
rement dans un tube de nickel. Utilisant de l'hydrogène 
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comme support gazeux, ils obtiennent une intensité 
de 0,5 microampère en ions lithium. Un tel dispositif 
pourrait rendre de grands services,.en spectroscopie de 
masses, pour l’analyse de certains échantillons solides. 


Nous sommes heureux d'exprimer notre très respec- 
tueuse gratitude à M. le Directeur Général G. Dupouy, 
pour l'intérêt constant qu'il a pris à cette étude, et 
pour lappui qu'il nous a prodigué. 


Manuscrit reçu le 17 août 1960. 


Note ajoutée aux épreuves, — Signalons également la 
source H. F. de Guidini [Gu-1] fonctionnant en impulsions 
successives. et qui a fait l’objet d’une note intérieure au 
CG. E. A. Cette source, avec une fréquence d’excitation de 
20 MHz, peut fournir jusqu’à 33 mA de protons, pour une 
puissance instantanée de 10 kW et un débit gazeux de 
52 cm® T. N. P. par heure. 
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LETTRE A LA RÉDACTION 


DISPOSITIF DE DÉTERMINATION SOUS VIDE 
DE LA CONDUCTIBILITÉ ÉLECTRIQUE 
ET DE LA CONSTANTE DE HALL 
DES LAMES MINCES MÉTALLIQUES 
ET SEMI-CONDUCTRICES 


Par S. MARTINUZZI, 


Laboratoire de Physique P. C. B., 
Faculté des Sciences, Alger. 


Les lames minces dont on étudie les propriétés élec- 
triques sont préparées par évaporation thermique, 
dans l'enceinte d’un évaporateur permettant 
d'atteindre et de maintenir un vide de l’ordre de 
10—$ mm de mercure. 


1. Appareïllage. — Il est représenté par la figure 1. 
Sur le support, où s’effectue la condensation de l’élé- 
ment ou du composé étudié, sont préalablement dépo- 


Fi. 4. — Le dispositif, 
sorti de l’enceinte après une évaporation. 


sées par évaporation sous vide, des électrodes en 
argent [1]. Leur disposition géométrique est donnée 
par la figure 2. Elles permettent la détermination de la 
conductibilité électrique et de la constante de Hall, 
sous le vide de la préparation. A cet effet, un aimant 
permanent placé à l’intérieur de l’enceinte peut être 
assujetti de l'extérieur à un mouvement hélicoïdal dis- 
continu alternant avec un mouvement de translation. 
Ce mécanisme sert à éloigner l’aimant de la lame 
pendant l’évaporation, ou à renverser le sens du champ 
magnétique. 


Deux fours permettent, outre une variation de la 
température interne, le dégazage du support et 
de l’enceinte, ainsi que des recuits prolongés sous vide 
jusqu’à 200 0C. 

Le premier, amovible, n’est utilisé que pour le déga- 
zage du support : le second, constitué par deux cy- 
lindres de laiton contenant les résistances chauffantes, 
peut porter, à 200 0C, toute l’enceinte (dont les parois 
intérieures sont argentées, et les joints en néoprène ou 
du type SIL). Les températures sont évaluées par des 
couples chromel-alumel. Enfin la détermination de 
l'épaisseur des lames minces obtenues se fait par pesée 
à la microbalance d’une lamelle de verre témoin avant 
et après la projection. 


2. Technique des mesures. — Le champ magnétique 
et le courant électrique utilisé sont continus. Les di- 
mensions de la lame sont telles que le rapport longueur 
sur largeur soit égal à 3, ce qui simplifie les calculs de 
la constante de Hall (on se place pratiquement dans 
le cas du conducteur infini) [2]. 

La connaissance de l’intensité du courant traversant 
la lame, et de la tension aux bornes des électrodes 3 et 4 
de la figure 2 donne la résistance électrique d’une 


F1a. 2, — Disposition géométrique des électrodes et mon- 
tage correctif de la tension de Hall proposé par Lambert. 


longueur connue (1 cm) de la lame mince, Il est néces- 
saire d'effectuer chaque fois deux mesures en inversant 
le sens du courant. 

Pour la constante de Hall, il suffit de mesurer l’in- 
tensité du courant, et la tension aux bornes des élec- 
trodes 1 et 4 de la figure 2. Mais il faut pouvoir éliminer 
les effets parasites TR-drop, Nernst, Ettinghausen, et 
Righi-Ledue. L'emploi de la méthode proposée par 
O. Lindberg [3], qui consiste en quatre mesures consé- 
cutives faites en inversant chaque fois, soit le sens du 
champ, soit le sens du courant, donne de bons résultats. 
Toutefois, comme les lames minces étudiées ont tou- 
jours des résistances élevées, l'effet TR-drop est très 
marqué et il convient de l’éliminer directement. La 
méthode proposée par Lambert[4],schématisée figure 2 
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est préférable à celle d’Aigrain [5], car elle permet cette 
élimination quelle que soit la résistance /? 2, 2° com- 
prise entre les deux électrodes 2 et 2’ de la figure 2. 
Ce dispositif est actuellement utilisé pour l'étude des 
lames minces d’arsenic, de gallium et d’arséniure de 
gallium. 

Lettre reçue le 8 décembre 1960. 
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REVUE DES LIVRES 


Brun (E. A.), MARTINOT-LAGARDE (A.), MAT4IEU (J.), 
Mécanique des fluides. Tome 1 : Généralités — Statique 
— Écoulements à une dimension — Conduites. (1 vol. 
25 /16 cm, 581 pages, Dunod, Paris, 1960, 58 N. F.) 

Ce cours de mécanique des fluides s’adresse spécialement 
aux élèves des facultés des Sciences et des Grandes Écoles 
d'Ingénieur : il contient en particulier les éléments des 
certificats de Mécanique des fluides I et IT institués par la 
nouvelle licence. 

L’avertissement nous apprend qu’il est prévu 3 tomes qui 
traiteront plus spécialement : 

Tome I. — Rappel des fondements de la mécanique, de 
la thermodynamique, de l’acoustique et de l’optique, envi- 
sagé du point de vue des applications à la mécanique des 
fluides. — Statique. — Écoulement à une dimension. — 
Pertes de charges dans les canalisations. 

Tome II. — Écoulement à 2 dimensions et à potentiel de 
vitesse des fluides incompressibles. — L’aérodynamique 
des grandes vitesses dans les domaines subsoniques, super- 
soniques et hypersoniques. — Phénomènes de transport 
de quantité de mouvement, de chaleur et de masse, princi- 
palement ceux qui ont pour siège, la couche limite. 

Tome III. — Méthodes de mesure et domaine d’appli- 
cation de la mécanique des fluides : machines hydrauliques, 
lubréfaction, entrainement des solides, aérosols, rhéologie, 
aéroélasticité.…. 

Le tome I qui nous est présenté ici traite des généralités 
(chap. 1), de la statique (ch. 2), et des écoulements à une 
dimension (ch. 3). 3 

Dans le chapitre 1, une Jongue partie est consacrée à la 
particule fluide, l'hypothèse de continuité et ses limitations 
(qui prennent une importance considérable en superaéro- 
dynamique). Toutes les grandeurs caractérisant les pro- 
priétés d’un fluide : compressibilité, dilatabilité, viscosité, 
conductibilité, diffusibilité, sont ensuite définies dans le 
cadre de leur utilisation future. Un bref rappel de la théorie 
cinétique des gaz permet de justifier au moins qualitative- 
ment les relations entre ces grandeurs caractéristiques et les 
grandeurs de base de la théorie cinétique. Après les rappels 
de mathématiques, de thermodynamique, des lois de simi- 
litudes et de leurs applications, l’auteur termine les géné- 
ralités par la description et la justification sommaire des 
néthodes optiques classiques et récentes d’études expéri- 
mentales des écoulements. A la fin de chaque sous-chapitre 
sont rassemblées des séries d'exercices très simples qui allè- 
gent l’exposé principal et en facilitent la compréhension. 

Le 2e chapitre du livre reprend l’exposé des notions clas- 
siques de statique des fluides. Il faut noter dans cette partie 
une série d’excellents exercices sur l’hydrostatique. 

Le chapitre 3 reprend d’abord les équations générales de 
la dynamique, puis précise les cas où il est permis de ne 
considérer schématiquement qu’une variable d’espace, et 
avec quelle approximation. 

Les écoulements unidimensionnels sont ensuite systéma- 
tiquement étudiés : 

— équation de Bernouilli et ses divers aspects. Une tren- 
taine d'exercices suivent l’exposé de cette partie et per- 
mettront à l’élève de manier avec virtuosité l’équation de 
Bernouilli dans tous les cas possibles. 
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— équations de Saint Venant à masse volumique varia- 

le ; 

— équations de Hugoniot et les écoulements subsoniques 
soniques et supersoniques dans les tuyères ; 

— les écoulements avec frottements et transmission de 
chaleur (ondes de choc et de combustion) ; 34 

— les écoulements dans les canaux ouverts qui trouvent M 
ici logiquement leur place si l’on considère les analogies ” 
entre cet écoulement et celui des fluides compressibles dans 
les tuyaux. 

Le dernier chapitre est consacré aux écoulements vis- 
queux et aux différents types de régime de ces écoulements. 

Pour l’ensemble de ce livre, plus de 200 exercices, tous 
extrêmement simples, mais test excellent pour l’élève de. 
son assimilation active au cours. Cette présentation et cet 
effort porté par les auteurs sur le choix des exercices confère. 
à cet ouvrage un intérêt pédagogique extrême. 

S. PAUTHIER. 
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La théorie des gaz neutres et ionisés. (Cours de l’école ” 
d’été de physique théorique des Houches) 1 vol. relié de. 
469 pages, 17 X 24 cm, prix : 60 N. F. Hermann, 
Paris, 1960. 

Il s’agit de recueil de cours, en langues anglaise et fran- 
çaise, à l’usage des chercheurs, rendant compte des der- 
nières et plus récentes recherches. Il ne s’agit pas d’un 
ouvrage pour étudiant, ni d’un ouvrage de culture générale. 
Voici la liste des articles : 

MonNTROLL (G. W.). — Questions de mécanique statis- 
tique des particules en interaction. I. Ensembles de Gibbs 
et fonction de partition. IL. Théorie des intégrales d’amas 
(cluster). LIT. Théorie des phénomènes de transport. IV. 
Déplacements aléatoires dans les espaces discrets. V. Tran- 
sitions de phase, problème d’Ising. 

VAN HOVE (L.). — Mécanique statistique des phénomènes 
de déséquilibre. (Équations de Boltzmann et de Pauli). 

. Dercroix (J. L.). — Théorie microscopique des gaz 

ionisés. z 
I. Fonctions de distribution simples et multiples. IL.Ani- « 

sotropie et moments de distribution simple. III. Équation 

d'évolution microscopique. IV. Étude des corrélations dans 
un plasma. V. Rayonnement des plasmas. ; 

KruskAL (M. A.). — L’hydromagnétique et la théorie # 
des plasmas dans un champ magnétique fort, et les prin- « 
cipes d'énergie pour l’équilibre et la stabilité. 

B) Théorie asymptotique de systèmes d'équations difté- « 
rentielles ordinaires avec toutes leurs solutions presque 
périodiques. 

C) Amortissement de Landau. | 

KAUFMAN (A.). — Théorie du transport pour les plasmas 
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